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FOREWORD 


This  report  was  prepared  fay  Aerotherm  Corporation  under  USAF 
Contract  F04611-68-C-0086  and  presents  the  development  and  demon¬ 
stration  of  total  heat  flux  and  radiative  heat  flux  measurement 
systems  for  use  in  ablative  nozzles.  The  work  was  administered 
under  the  direction  of  the  Air  Force  Rocket  Propulsion  Laboratory 
with  Lt.  John  R.  Ellison  and  Mr.  Harold  I.  Binder  as  the  Project 
Officers . 


Mr.  John  w.  Schaefer  was  the  Aerotherm  Corporation  overall 
program  manager  and  principal  investigator  in  the  area  of  test 
evaluation.  Mr.  Duane  L.  Baker  was  the  principal  investigator 
in  the  areas  of  accuracy  analysis,  measurement  systems  design, 
and  nozzle  design.  Mr.  Mitchell  R.  Wool  was  the  principal  inves¬ 
tigator  in  the  areas  of  data  reduction  techniques  and  analysis  of 
test  results.  Mr.  John  J.  Reese,  Jr.  provided  support  in  the  test 
evaluation  program. 


The  authors  wish  to  acknowledge  the  contribution  of  Mr.  D. 
Larry  Jones,  President,  Heat  Technology  Laboratory,  Inc.,  for 
his  support  in  the  design  of  the  radiation  heat  flux  systems  and 
the  contribution  of  the  AFRPL  test  crew  who  performed  the  motor 
firings  under  this  contract. 


This  technical  report  has  been  reviewed  and  is  approved 


Charles  R.  Cooke 

Chief,  Solid  Rocket  Division 

Air  Force  Rocket  Propulsion  Laboratory 


ABSTRACT 


Total  heat  flux  and  radiative  heat  flux  measurement  systems 
vero  developed  for  use  in  ablative  solid  rocket  nozzles.  The 
total  heat  flux  system  is  a  thermocouple  instrumented  plug  of  the 
same  material  and  layup  angle  as  the  ablative  part  in  which  it  is 
installed.  The  total  heat  flux  through  a  firing  is  determined 
from  an  analysis  of  the  measured  in-depth  temperature  histories. 

The  radiative  heat  flux  system  is  a  narrow  view  angle,  Gardon-type 
radiometer  installed  in  a  cavity  below  the  surface  of  the  ablative 
part.  The  total  radiative  heat  flux  through  a  firing  is  determined 
from  a  calibration  curve  relating  millivolt  output  and  total  radia¬ 
tion  flux  of  a  uniform  calibration  source.  The  ablative  materials 
considered  under  the  program  were  MX4926  carbon  cloth  phenolic  and 
MX2600-96  silica  cloth  phenolic.  Detailed  accuracy,  design,  fab¬ 
rication,  and  installation  studies  were  performed  to  define  the 
system  configurations.  A  preliminary  evaluation  program  was  per¬ 
formed  in  which  the  measurement  systems  were  installed  in  two- 
dimensional  nozzles  and  successfully  tested  in  simulated  rocket 
environments  using  an  arc  plasma  generator.  Final  systems  evalua¬ 
tion  and  demonstration  was  performed  through  two  nozzle  firings  on 
the  AFRPL  40-inch  char  motor.  A  total  of  six  total  flux  systems 
and  three  radiation  flux  systems  were  tested  in  these  firings.  A 
high  flame  temperature  aluminized  propellant,  APG  112D,  was  used 
and  peak  pressures  in  both  firings  were  approximately  1000  psia 
with  firing  durations  of  approximately  60  seconds.  Successful 
performance  of  the  two  measurement  systems  was  demonstrated. 
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SECTION  1 


INTRODUCTION 


A  knowledge  of  the  convective  and  radiative  heat  flux  incurred 
by  ablative  materials  in  solid  rocket  nozzles  is  important  in  the 
design  of  these  nozzles  and  in  the  analysis  of  material  response. 

No  practical  and  effective  methods  for  measuring  these  heat  fluxes 
have  been  available,  however.  The  objective  of  the  program  de¬ 
scribed  herein  was  the  development  of  measurement  systems  for  use 
in  ablative  nozzle  materials.  Total  and  radiation  heat  flux  inci¬ 
dent  to  the  material  surface  are  measured  by  these  systems,  convec¬ 
tive  heat  flux  is  defined  by  difference. 


The  total  heat  flux  measurement  system  is  a  thermocouple 
instrument  plug  of  the  same  ablative  material  and  reinforcement 
tape  layup  angle  as  the  nozzle  part  in  which  it  is  installed.  The 
primary  instrumentation  is  several  fine-wire  thermocouples  which  are 
installed  in  the  plug  to  accurately  measure  the  thermal  response  of 
the  ablative  material.  The  total  surface  heat  flux  is  defined  from 
an  analysis  of  the  in-depth  temperature  measurements. 


The  radiation  flux  measurement  system  is  a  narrow  view  angle 
Gardon-type  radiometer  which  is  installed  in  a  cavity  in  the  abla¬ 
tive  nozzle.  The  radiometer  is  gas  purged  and  employs  a  sapphire 
window  to  eliminate  convective  heating  and  particle  deposition. 

The  millivolt  output  from  the  narrow  view  angle  radiometer  is  con¬ 
verted  directly  to  a  hemispherical  radiative  heat  flux  utilizing 
a  calibration  based  on  measurements  of  a  known  uniform  radiative 
energy  source.  This  simplified  technique  for  obtaining  the  hemi¬ 
spherical  radiative  heat  flux  was  used  because  of  the  difficulty 
in  directly  measuring  the  total  hemispherical  radiative  heat  flux 
in  an  ablative  nozzle  and  because  of  cost  considerations  in  making 
such  measurements. 


The  program  was  divided  into  three  phases: 


Phase  I 
Phase  II 
Phase  III 


Study  and  Design 

Fabrication,  Test,  and  Evaluation 
Rocket  Nozzle  Static  Firings 


The  results  from  these  phases  are  summarized  in  Section  2.  A  more 
detailed  description  of  the  Phase  I  effort  is  described  in 
Sections  3,  4  and  5.  Section  3  presents  both  an  evaluation  of  the 
measurement  systems  studied  and  a  general  description  of  the  total 
and  radiative  systems  used  in  the  present  program.  Section  4  presents 
the  results  of  a  detailed  systems  study  and,  also,  the  design  which 
was  formulated  from  this  study.  The  systems  study  included  am  error 
analysis  of  the  primary  measurements  for  both  the  total  and  radiative 
heat  flux  systems.  In  addition,  this  section  describes  the  instru¬ 
ment  fabrication  and  installation  procedures,  and  the  accessory 
recording  instrumentation  required  to  make  the  primary  measurement. 
Section-  5  presents  the  data  reduction  procedures  for  both  measurement 
systems.  Section  6  presents  the  Phase  II  effort  which  consisted  of 
the  test  evaluation  of  the  measurement  systems  in  a  simulated  en¬ 
vironment  which  was  generated  by  the  Aerotherm  arc  plasma  generator 
(APG) .  The  Phase  III  effort  is  presented  in  Section  7.  The  per¬ 
formance  for  several  measurement  .systems  of  both  types  is  presented 
for  two  ablative  nozzles  that  were  test  fired  on  the  AFRPL  40-inch 
char  motor.  Finally,  conclusions  and  recommendations  are  presented 
in  Section  8. 


SECTION  2 
SUMMARY 


Total  heat  flux  and  radiation  heat  flux  measurement  systems 
were  developed  for  application^ in  the  severe  environment  of  a  solid 
rocket  nozzle  firing.  These  systems  are  installed  in  the  ablative 
components  of  a  nozzle  and  are  designed  to  perform  continuously 
throughout  a  firing. 


The  total  heat -flux  system  is  a  thermocouple  instrumented  plug 
of  the  same  material  and  with  the  same  layup  angle  as  the  ablative 
part  in  which  it  is  installed.  Tungsten  5  percent  rh'enium/tungsten 
26  percent  rhenium  and  Chrome 1 /A lumel  thermocouples  are  employed  at 
five  in-depth  locations  in  the  plug.  The  total  heat  rlux  is  deter¬ 
mined  by  a  computer  code  analysis  which  considers  the  measured  in- 
depth  temperatures  and  post-test  surface  recession  to  (1)  define  the 
surface  recession  and  surface  temperature  histories  and  (2)  define 
all  the  surface  heat  flux  terms  such  that  the  total  heat  flux  is  the 
only  unknown  in  the  surface  energy  balance.  The  total  flux  system 
may  be  used  in  any  material  for  which  the  thermal  and  ablative  prop¬ 
erties  are  well  characterized.  MX4926  carbon  cloth  phenolic  and 
MX2600-96  silica  cloth  phenolic  were  considered  in  this  program. 


The  radiative  heat  flux  system  is  a  narrow  view  angles  Gardon- 
type  radiometer.  The  narrow  view  angle  is  required  so  that  the 
radiometer  can  be  recessed  in  a  cavity  below  the  ablative  surface 
to  protect  it  from  the  combustion  products  and  to  allow  for  surface 
recession.  The  constantan  foil  sensor  is  protected  by  a  sapphire 
lens  which  is  in  turn  protected  by  a  nitrogen  gas  purge.  The  total 
(hemispherical)  radiation  heat  flux  is  determined  directly  from  a 
calibration  curve  relating  the  sensor  output  to  the  total  radia¬ 
tion  flux  of  a  uniform  source.  Tht  radiative  flux  system  may  be 
used  in  any  material  or  nozzle  configuration  providing  the  instal¬ 
lation  does  not  result  in  a  high  thermal  load  on  the  radiometer 
due  to  conduction  from  the  material.  The  radiometers  were  installed 
in  MX4926  carbon  phenolic  in  this  program. 


The  program  encompassed  detailed  accuracy,  desigrt,  fabrication, 
and  installation  analyses;  definition  of  the  data  reduction  tech¬ 
niques  and  procedures;  preliminary  evaluation  program  in  which  the 
measurement  systems  were  tested  under  simulated  rocket  environment 
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conditions  in  an  arc  plasma  generator;  and  final  evaluation  and 
demonstration  program  in  which  the  measurement  systems  were  tested 
in  two  40-inch  char  motor  test  firings.  This  program  sequence  ap¬ 
plied  to  both  measurement  system  types;  also  each  char  motor  nozzle 
in  the  final  evaluation  and  demonstration  program  incorporated  both 
the  total  and  radiative  flux  systems. 


For  the  total  flux  system,  the  analysis  of  potential  errors 
considered: 

•  Accuracy  of  the  indicated  temperature 

•  Accuracy  of  the  thermocouple  output 

•  Accuracy  of  the  thermocouple  location 

The  errors  associated  with  temperature  field  distortions  due  to  the 
presence  of  the  thermocouple  and  the  hole  that  accepts  the  thermo¬ 
couple  and  to  the  contact  resistance  between  the  "flat  plate"  plug  ^ - 
geometry  and  the  axisymmetric  geometry  of  the  part  in  which  it  is 
installed  were  defined.  The  errors  in  thermocouple  output  due  to 
electrical  shorting  across  the  conducting  char  and  due  to  thermo¬ 
couple  calibration  uncertainties  and  possible  high  temperature  cali¬ 
bration  shifts  were  defined.  Finally,  the  uncertainties  in  the  X-ray 
photography  measurement  of  the  thermocouple  locations  in  the  plug 
were  defined.  The  potential  uncertainty  in  the  total  heat  flux 
which  resulted  from  these  accuracy  analyses  ranged  from  -11  to  5  per¬ 
cent  for  the  thermocoupled  MX4926  plug  and  from  -20  to  13  percent 
for  thef  thermocoupled  MX2600-96  plug.  The  results  of  the  accuracy 
analyses  were  incorporated  into  the  plug  design. 


The  thermocouple  locations  and  thermocouple  types  are  selected 
based  on  the  anticipated  response  of  the  ablative  material  at  the 
installation  location.  Bare  wire  5  mil  thermocouples  are  used  and 
are  fed  through  a  strapped  hole  along  an  isothermal  plane  in  the  plug 
such  that  the  bead  (thermocouple  junction)  rests  against  the  step 
at  the  plug  centerline.  The  thermocouple  wires  are  fed  through 
beryllia  sleeving  along  each  side  of  the  plug  to  prevent  electrical 
shorting.  The  thermocouple  wire/lead  wire  connection  is  made  on 
pins  installed  in  the  bottom  of  the  total  flux  plug.  The  plug  is 
bonded  into  the  ablative  pa  t  with  epoxy  resin. 


For  the  radiative  flux  system,  the  accuracy  analysis  considered: 

•  Errors  due  to  thermal  effects 

•  Errors  due  to  optical  effects 

•  Calibration  uncertainties 


f 
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The  errors  associated  with  a  change  in  the  radiometer  body  tempera¬ 
ture  and  its  effect  on  the  sensor  calibration  and  radiation  inter¬ 
change  were  defined.  In  addition,  the  errors  due  to  the  wavelength 
response  of  the  sapphire'  protective  lens  and  the  high  emissivity 
sensor  coating  as  it  relates  to  wavelength  differences  between  the 
calibration  source  and  the  measured  source  were  defined.  Finally, 
the  uncertainties  in  the  sensor  calibration  and  in  the  assumption 
of  a  uniform  radiation  flux  in  the  rocket  nozzle  were  defined.  The 
uncertainty  in  the  measured  radiative  heat  flux  resulting  from  the 
accuracy  analysis  varied  from  -6  to  -2  percent  which  could  be  re¬ 
duced  to  -3  to  2%  percent  by  making  the  necessary  corrections  for 
the  body  temperature . 


The  radiometer  employs  a  heat  sink  body  to  maintain  an  accept¬ 
able  sensor  temperature  level  without  the  added  complication  of 
water  cooling.  The  unit'  is  nitrogen  gas  purged  to  protect  the  sap¬ 
phire  lens  from  convective  heating  and  oxide  particle  contamination. 

The  radiometer  is  installed  in  a  cavity  in  the  ablative  part  and  is 
recessed  below  the  exposed  surface  to  protect  it  from  the  solid 
propellant  combustion  products..  The  design  allows  for  minimum  thermal 
contact  between  the  radiometer  and  the  high  temperature  region  of  the 
ablative  part.  The  unit  is  potted  in  place  with  epoxy  resin  or  threaded 
into  the  nozzle  steel  shell. 


‘  Data  reduction  for  the  total  heat  flux  measurement  system  en¬ 
compasses  : 

•  Iterative  in-depth  conduction  solutions  using  the  Charring 
Material  Ablation  (CMA)  computer  code;  final  solution  pro¬ 
vides  a  close  match  of  the  measured  and  calculated  in-depth 
temperature  histories  and  defines  the  surface  recession, 
surface  temperature,  and  some  of  the  surface  heat  flux  terms. 

r 

•  Definition  of  the  thermochemical  surface  heat  flux  terms 
from  the  above  results  and  using  the  Aerotherm  Chemical 
Equilibrium.  (ACE)  computer  code. 

•  Summation  df  the  various  terms  of  the  surface  energy  balance 
to  define  total  heat  flux. 

The  final  result  is  the  total  heat  flux  history  at  the  measurement 
location  throughout  the  nozzle  firing.  Data  reduction  for  the  radia¬ 
tive  heat  flux  system  is  accomplished  by  simply  relating  the  mea¬ 
sured  millivolt  output  of  the«fadiometer  to  the  calibration  curve 
of  millivolts  vs  total  radiative  heat  flux.  The  final  result  is 
the  total  radiative  heat  flux  history  at  the  measurement  location 
throughout  the  nozzle  firing. 
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In  the  preliminary  evaluation  test  program,  six  total  flux  sys¬ 
tems  and  two  radiative  flux  systems  were  installed  in  carbon  phenolic 
and  silica  phenolic  two-dimensional  nozzles.  These  nozzles  were  test 
fired  in  an  arc  plasma  generator  under  conditions  which  simulated  a 
rocket  nozzle  firing  and  which  also'  Allowed  an  independent  measure¬ 
ment  of  total  and* radiative  heat  flux,  surface  temperature  history, 
and  suxface  recession  history.  The  measurement  systems  were  de¬ 
signed  to  closely  duplicate  those  used  in  an  actual  rocket  nozzle. 

The  units  performed  properly  from  both  functional  and  measurement 
standpoints. 


In  the  final  evaluation  and  demonstration  program,  six  total 
flux  systems  and  three  radiative  flux  systems  were  installed  in  two 
ablative  pozzies  which  were  test  fire<J  on  the  AFRPL  40-inch  char 
motor.  The  seguence  of  the  test  firings  was  such  that  the  second 
nozzle  firing  incorporated  modifications  suggested  by  the  first 
firing.  A  summary  of  system  locations  and  materials  in  both  nozzles 
is  presented  below: 


Nozzle 


System 
Type  " 


System 
Location , 
Area  Ratio 


Nozzle  Ablative 
Material 


Total 


Radiative 

Total 


.Radiative 


-2.0 

1.0 

2.0 

4.0 

1.0 


3.0 

1.0 

3.0 


MX4926  Carbon  Phenolic 


MX2600-96  Silica  Phenolic 
MX^926  Carbon  Phenolic 


A  highly  aluminized  propellant?,  APG  112D,  with  a  flame  temperature 
of  7020°R  at  1000  psia  was  us^d  in  both  test  firings.  The  maximum 
chamber  pressures  were  1030  and  930  psia  for  the  first  and  second 
firings,  ^respectively;  the  chamber  pressures  were  highly  regres¬ 
sive  due  to  the  surface  recession  of  the  ablative  throat  with  final 
pressures  in  both  firings  of  about  300  psia.  The  test  firing  times 
were  .just  under  one  minute  for  both  firings. 

•The  total  heat  flux  measurement  systems  performed  satisfactorily 
and  repeatably  in  both  test  firings  and  provided  an  accurate  defini¬ 
tion  of  the  total  heat  flux  histories.  In  the  first  firing,  a 
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failure  of  the  nitrogen  purge  gas  to  the  radiative  flux  system  caused 
a  malfunction  of  the  radiometer  early  in  the  firing.  Modifications 
were  made  to  the  purge  gas  system  to  eliminate  this  problem  in  the 
second  firing.  The  radiometer  at  an  area  ratio  of  3.0  in  the  second 
firing  performed  satisfactorily  throughout  the  firing  and  provied 
an  accurate  definition  of  the  total  radiation  heat  flux  history.  The 
view  of  the  radiometer  at  an  area  ratio  of  1.0  was  apparently  par¬ 
tially  obstructed  through  most  of  the  firing  and  therefore  the  re¬ 
sults  at  this  location  are  questionable. 


The  end  result  of  the  program  was  the  successful  development 
of  total  and  radiative  heat  flux  measurement  systems  for  use  in 
ablative  solid  rocket  nozzles. 


/ 

/ 

/ 

/ 

/  ■ 
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SECTION  3 


MEASUREMENT  SYSTEM  CONCEPTS 


The  proper  selection  of  total  and  radiation  heat  flux  measure¬ 
ment  systems  requires  consideration  of  the  response  of  ablative 
materials  such  as  the  cloth  phenolics  to  solid  propellant  combustion 
environments.  This  material  response  and  the  environmental  bound¬ 
ary  "conditions  to  which  the  material  is  exposed  are  illustrated 
schematically  in  Figure  3-1.  The  surface  recedes  due  to  surface 
thermochemical  reactions  and/or  thermomechanical  erosion.  The 
material  degrades  in  depth  forming  a  char  through  which  the  gase¬ 
ous  degradation  products  flow.  This  response  is  typically  tran¬ 
sient,  steady  state  ablation  conditions  being  achieved  only  late 
in  a  rocket  firing  if  at  all. 


RADIATION 


CHEMICAL  species 


FREE 

STREAM 


BOUNDARY 

LAYER 


CHAR 


PYROLYSIS 

ZONE 


VIRGIN 

MATERIAL 


BACKUP 

MATERIAL 


FIGURE  3-1  CHARACTERIZATION  OF  THE  SOLID  PROPELLANT  ROCKET 

ENVIRONMENT  AMD  THE  RESPONSE  OF  A  CHARRING  ABLATOR 
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As  shown  by  this  figure,  the  heat  flux  from  the  combustion  products 
to  the  nozzle  wall  is  a  combination  of  a  convective  flux,  a  chemical 
species  diffusion  flux,  and  a  radiative  flux.  The  convective  flux 
represents  sensible  energy  transfer  to  the  wall,  and  the  diffusion 
flux  represents  chemical  energy  transfer  to  the  wall.  These  two 
terms  are  generally  combined  and  referred  to  as  the  convective  flux. 

The  objective  of  the  instrumentation  systems  is  to  measure  both  the 
total  and  radiative  heat  flux  to  the  nozzle  wall.  The  difference 
between  these  two  heat  fluxes  represents  the  convective  flux  which 
is  a  combination  of  the  sensible  and  chemical  energy  transfer.  The 
ground  rules  for  selection,  the  approaches  considered,  and  the 
selected  approach  are  presented  in  Section  3.1  for  the  total  heat 
flux  system  and  in  Section  3.2  for  the  radiative  heat  flux  system. 


3.1  TOTAL  HEAT  FLUX  SYSTEM 

The  surface  response  of  an  ablative  material  is  indeed  compli¬ 
cated  and  the  events  which  make  up  this  response  are  interrelated 
(Figure  3-1) .  In  order  to  measure  the  total  heat  flux  to  an  ablative 
material,  the  measurement  system  must  therefore  represent  a  neg¬ 
ligible  disturbance  to  the  surface  and  in-depth  response  of  the 
material.  A  system  which  employs  a  plug  of  the  ablative  material 
that  is  instrumented  with  fine  wire  thermocouples  satisfies  this 
requirement.  The  in-depth  thermal  response  of  the  material  is 
measured  and  the  surface  response  including  total  heat  flux  is  deter¬ 
mined  from  this  measurement.  This  concept  best  satisfies  the  require¬ 
ment  of  an  accurate  measurement  of  total  heat  flux  to  an  ablative 
material  and  was  therefore  the  concept  selected.  The  system  advant¬ 
ages  are  summarized  as  follows: 

9 

•  negligible  disturbance  to  the  in-depth  and  surface 
response  of  the  material;  measured  temperature  is 
that  of  the  material 

g  data  are  obtained  continuously  through  a  firing;  total 
heat  flux  is  defined  throughout  a  firing 

•  design,  fabrication  technique,  data  reduction  procedure, 
and  calibration  procedure  are  within  the  current  state- 
of-the-art 

•  necessary  recording  instrumentation  and  data  reduction 
tools  are  readily  available 


The  measurement  system  design  is  illustrated  in  Figure  3-2, 
and  its  installation  in  a  rocket  nozzle  is  shown  schematically  in 
Figure  3-3.  The  measurement  system  consists  of  a  plug  of  the  same 
ablative  material  and  layup  angle  as  that  of  the  part  in  which  it 
is  installed.  The  plug  is  instrumented  with  five  thermocouples 
such  that  temperature  measurements  to  5400°F  are  possible.  A 
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SCHEMATIC  OF  TOTAL  HEAT  FLUX 
MEASUREMENT  SYSTEM 
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FIGURE  3-3  SCHEMATIC  OF  CHAR  MOTOR  NOZZLE  WITH 

measurement  system  installed 


continuous  measurement  of  the  in-depth  thermal  response  of  the 
ablative  material  is  therefore  made  throughout  a  firing.  An  in- 
depth  conduction  solution  which  is  based  on  these  measurements  and 
considers  the  response  mechanisms  depicted  in  Figure  3-1  then  pro¬ 
vides  the  required  total  heat  flux  to  the  material  surface.  The 
two  ablative  materials  selected  for  measurement  system  evaluation 
are  MX  4926  carbon  phenolic  and  MX  2600-96  silica  phenolic.  These 
materials  are  representative  of  those  of  interest  for  the  measure¬ 
ment  system  application  and  also  are  accurately  characterized  in 
terms  of  their  properties  and  response  mechanisms . 


Other  measurement  system  concepts  were  also  considered  but  none 
was  as  attractive  as  that  selected.  The  major  drawback  of  most  con¬ 
cepts  is  the  finite  size  of  the  sensor  and  the  resultant  disturbance 
of  the  surface  and  in-depth  response.  Measurement  techniques  in 
this  category  include  the  transient  slug  calorimeter,  the  steady- 
state  Gardon-type  gauge,  and  the  conduction  calorimeter  (based  on 
the  AT  across  a  low  conductivity  wafer) .  With  these  techniques ,  the 
total  heat  flux  measured  is  not  the  heat  flux  to  the  ablative  material 
and,  because  of  the  complicated  interaction  between  the  sensor  and 
the  material  in  which  it  is  installed,  the  relation  between  the 
measured  and  desired  fluxes  is  very  difficult  if  not  impossible  to 
define.  Some  of  the  problems  in  this  regard  are  the  discontinuity 
in  the  energy  and  mass  boundary  layers  associated  with  a  dissimilar 
material  in  the  ablative  part,  the  disturbance  to  the  thermal  con¬ 
duction,  pyrolysis  gas  flow,  and  thermal  decomposition  of  the 
ablative  material,  the  dissimilar  surface  recession  of  the  two 
materials,  and  the  availability  of  data  at  only  a  discrete  time 
or  times .  Note  that  these  comments  also  apply  to  the  use  of 
several  sensors  of  the  same  type  stacked  together  to  form  a 
measurement  system  plug.  Another  general  approach  considered  is 
the  direct  measurement  of  surface  temperature.  Optical  techniques 
are  not  possibly  because  of  the  optically  thick  particle  cloud  in 
solid  propellant  combustion  products.  The  various  techniques  for 
measuring  surface  temperature  using  a  thermocouple  approach  have 
proved  to  be  unreliable  (Reference  3-1)  and  therefore  are  not 
attractive.  The  thermocouple  instrumented  plug  of  the  ablative 
material  avoids  the  above  problem  areas  and  was  therefore  the 
logical  selection  for  the  measurement  system. 


3.2  RADIATIVE  HEAT  FLUX  SYSTEM 

The  total  radiative  heat  flux  to  the  surface  of  an  ablative 
material  in  a  solid  rocket  nozzle  results  from  the  aluminum  oxide 
particles,  from  various  gas  phase  species,  and  from  the  opposite 
walls.  This  radiative  flux  is  a  function  of  parameters  such  as 
particle  density,  particle  temperature,  particle  emissivity,  and 
particle  optical  depth,  and,  with  the  exception  of  particle  emissiv¬ 
ity,  these  parameters  vary  with  nozzle  axial  location.  This  radiative 
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intensity  is  therefore  a  complicated  function  of  the  direction  of 
the  line  of  sight.  Because  of  this  conq>lication  and  because  of  the 
variation  in  the  chamber  conditions  during  the  motor  firing,  the 
following  ground  rules  were  established  tentatively  for  the  design 
of  the  radiative  heat  flux  measurement  system. 

a.  Continuous  measurement  through  a  firing  since  the 
chamber  pressure  will,  in  general,  be  changing 
with  time. 

b.  Measurement  of  total  radiative  heat  flux  to  the  nozzle 
wall  (as  opposed  to  spectral  and  narrow  angle  flux) . 

c.  Measurement  location (s)  at  the  same  area  ratio (s)  as 
the  ablative  plugs. 

Unfortunately,  (a)  and  (b)  are  incompatible  in  that  a  total  radiation 
flux  sensor  must  be  mounted  at  the  surface  for  a  180°  view  but  the 
severe  heat  flux  and  the  surface  recession  preclude  a  continuous 
measurement  capability.  An  attractive  solution  is  (1)  the  use  of 
a  narrow  view  angle  radiometer  which  is  recessed  below  the  ablative 
surface  for  continuous  measurement,  (2)  the  use  of  a  wide  view  angle 
(180°)  radiometer  which  becomes  exposed  at  the  surface  for  an  instan¬ 
taneous  measurement  before  it  is  subsequently  destroyed, '  and  (3)  the 
use  of  the  two  measurements  at  the  identical  times  to  calibrate  the 
relationship  between  the  total  and  narrow  angle  radiation  fluxes. 

The  development  of  two  radiation  measurement  systems  was  beyond  the 
scope  of  the  subject  program,  however,  and  therefore  the  continuous 
measurement  system,  the  narrow  view  angle  radiometer,  was  selected 
for  development. 


The  design  of  the  narrow  view  angle  radiometer  is  shown  in 
Figure  3-4,  and  the  installation  within  the  nozzle  is  shown  sche¬ 
matically  in  Figure  3-5.  Frcxn  this  figure,  one  notes  that  the 
radiometer  is  placed  in  a  cavity  below  the  initial  surface  of  the 
nozzle,  and,  thus,  it  is  never  directly  exposed  to  the  solid  propel¬ 
lant  combustion  products  even  though  surface  recession  occurs 
during  the  motor  firing.  The  narrow  view  angle  restriction  is 
required  to  eliminate  the  possibility  of  the  radiometer  viewing  the 
side  walls  of  the  cavity. 


The  basic  components  of  the  narrow  view  angle  radiometer  shown 
in  Figure  3-4  are 

Gardon-type  heat  flux  gauge 

blackened  constantan  sensing  disk 
copper  mounting  heat  sink  body 

copper  lead  wire  connected  to  conatantan  disk  center 
nitrogen  purge  gas  flow  system 
sapphire  lens 
copper  body 

thermocouple  to  measure  body  temperature 
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FIGURE  3-4  SCHEMATIC  OF  RADIATION  HEAT  FLUX 
MEASUREMENT  SYSTEM 


A  complete  description  of  the  analytical  background  of  a 
Gardon- type  heat  flux  gauge  is  given  in  Reference  3-2.  Simply  stated, 
the  temperature  difference  between  the  center  of  the  constantan  foil 
sensor  and  the  copper  heat  sink  body  results  in  a  thermocouple  emf 
which  varies  linearly  with  incident  heat  flux.  The  nitrogen  purge 
gas  provides  convective  cooling  of  the  sensor  body  and  the  lens. 

The  lens  shields  the  foil. disk  from  convective  heating  (or  cooling). 

In  this  configuration,  the  millivolt  response  of  the  Gardon-type 
heat  flux  #auge  represents  only  the  incident  radiative  heat  flux. 


In  addition  to  the  radiometer  technique  for  radiative  flux 
measurement,  a  spectral  approach  was  also  considered.  This  technique 
employed  the  measurement  of  the  spectral  emittance  of  the  sources 
over  a  narrow  view  angle  using  a  reference  source  and  a  spectrometer. 
This  approach  was  rejected  because  it  requires  sophisticate and 
delicate  instrumentation  and  requires  a  basic  knowledge  of  the  various 
sources  of  radiation  and  the  directional  variation  which  is  currently 
unavailable. 


\ 
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SECTION  4 


SYSTEMS  STUDY  AND  DESIGN 

\  ■ 

The  general  approach  to  be  used  in  making  the  total  heat  flux 
and  radiation  heat  flux  measurements  were  defined^in  Section  3.  A 
systems  study  of  these  approaches  for  the  purpose  of  defining 
the  errors  in  the  primary  measurements  was  performed  and  a  detailed 
design  was  formulated  which  minimized  these  errors.  The  results  of 
this  study  and  design  are  in  Section  4.1  for  the  total  heat  flux 
measurement  s’;satem  and  in  Section  4.2  for  the  radiative  heat  flux 
system. 


4.1  TOTAL  HEAT  FLUX  SYSTEM 

The  total  heat  flux  measurement  system  is  dependent  on  the 
accurate  measurement  of  in-depth  temperatures  using  thermocouples. 
Errors  which  are  associated  tfith  the  use  of  thermocouples  to  it  isure  - 
in-depth  temperatures  include  distortion  of  the  temperature  field, 
thermocouple  response  lag,  thermocouple  electrical  shorting, 
thermocouple  calibration,  and  accuracy  of  thermocouple  location. 

These  sources  of  error  are  discussed  in  Section  4.1.1.  Once  the 
errors  were  defined,  a  design  and  fabrication  analysis  was  performed 
which  defined  an  economical  measurement  system  configuration  that 
minimized  the  measurement  error.  This  analysis  is  presented  in 
Section  4.1.2  and  includes  a  discussion  of  thermocouple  locations  and 
types,  thermocouple  installation,  system  configuration  and  installation, 
and  peripheral  recording  instrumentation.  Note  that  even  though  the 
eiyror,  design,  and  fabrication  analyses  are  presented  separately, 
these  analyses  were  performed  concurrently  and  several  iterations 
were  required  before  a  final  design  Was  adopted. 

4.1.1  Error  Analysis 

The  primary  measurement  required  in  the  total  heat  flux 
measurement  system  concept  (Section  3.1)  is  temperature.  Its 
measurement  in  ablative  materials  can  be  the  source  of  significant 
errors.  The  problem  areas  related  to  temperature  measurement  with 
thermocouples  are  1)  the  thermocouple  perturbs  the  temperature  dis¬ 
tribution  which  would  otherwise  exist  in  the  material,  2)  the  thermo¬ 
couple  junction  temperaturq  lags  the  temperature  of  the  material  due 
to  heat  conduction  and  thermal  capacity  effects,  3)  the  thermocouple 
wires,  except  along  an  isotherm,  must  be  electrically  insulated  from 
each  other  in  order  to  prevent  any  discrepancy  between  the  electro¬ 
motive  force  (emf)  generated  at  the  thermocouple  junction  and  the 
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measured  emf,  and  4)  the  location  of  the  thermocouple  must  be  accu¬ 
rately  known.  The  first  three  problems  were  considered  in  the 
selection  of  the  measurement  system,  and  two  of  the  more  important 
design  objectives  were  to  minimize  problems  one  and  two  and  to  elim¬ 
inate  problem  three.  The  design  of  the  thermocouple  plug  is  not 
directly  influenced  by  the  fourth  problem,  and  thus,  this  problem 
can  be  solved  after  the  design  of  the  thermocouple  plug  has  been 
finalized.  The  purpose  of  the  error  analysis  presented  in  this 
section  is  to  determine  whether  or  not  these  design  objectives  are 
satisfied  and  to  define  the  overall  accuracy  of  the  primary  mea¬ 
surement,  temperature,  required  by  the  measurement  system. 


The  accuracy  of  the  thermocouple  temperature  (problems  1  and 
2)  is  established  in  Section  4. 1.1.1  followed  in  Section  2 

by  an  analysis  of  the  accuracy  of  the  thermocouple  emf  measure¬ 
ment.  In  this  last  section,  an  analysis  is  presented  for  the 
possible  shunting  effect  across  the  thermocouple  wires  (problem  3) , 
and  this  effect  is  combined  with  the  guaranteed  accuracy  of  the 
thermocouple  and  lead  wires.  The  problem  of  accurately  defining  the 
in-depth  thermocouple  location  is  discussed  in  Section  4. 1.1. 3. 
Finally,  the  errors  associated  with  the  individual  problems  are  com¬ 
bined  and  summarized  in  Section  4. 1.1. 4,  and  an  overall  temperature 
measurement  accuracy  is  defined. 


4. 1.1.1  Accuracy  of  Indicated  Temperature 

An  ideal  thermocouple  which  is  located  within  a  material 
would  generate  an  emf  corresponding  to  the  temperature  which  would 
exist  at  the  same  location  in  the  absence  of  the  thermocouple.  This 
temperature  and  corresponding  emf  will  be  referred  to  hereafter  as 
the  unperturbed  temperature  or  emf.  From  a  detailed  investigation 
of  the  thermocouple  plug  design  shown  in  Figure  3-2,  one  notes  three 
regions  which  may  cause  a  difference,  assuming  no  emf  error,  between 
the  measured  and  unperturbed  temperatures;  namely  1)  the  region 
incorporating  the  cylindrical  cavities  in  which  the  thermocouple  is 
located,  2)  the  region  along  the  sidewalls  of  the  plug  assembly, 
including  the  bondline,  and  3)  the  region  within  the  cavity  which 
includes  both  the  thermocouple  junction* and  the  thermocouple  wires. 
The  first  two  regions  which  will  cause  distortions  in  the  temperature 
field,  and  the  third  region  which  will  cause  a  response  lag  in  the 
temperature  of  the  thermocouple  junction  are  analyzed  in  the  follow¬ 
ing  sections.  In  these  sections,  only  the  tungsten  5%  rhenium/ 
tungsten  26%  rhenium  (W5Re/W26Re)  thermocouple  nearest  the  surface 
is  analyzed  in  detail.  This  thermocouple  location  was  chosen  for 
analysis  because  it  experiences  the  most  rapid  temperature  response, 
and,  as  a  result,  the  errors  associated  with  measuring  the  in-depth 
temperature  are  maximized. 


* 


4. 1.1. 1.1  Distortion  of  Temperature  Field.  In  order  to  define  the 
distortion  of  the  i  *mperature  field,  the unperturbed  temperature 
distribution  was  first  determined.  The  nozzle  configuration  used 
for  this  analysis  was  shown  schematically  in  Figure  3-3,  and  the  rep¬ 
resentative  analysis  locations  of  interest  were  at  the  nozzle  throat 
and  at  a  supersonic  area  ratio  of  four.  The  material  at  the  nozzle 
throat  is  carbon  cloth  phenplic  (MX  4926)  at  a  60°  layup  angle 
with  respect  to  the  nozzle  centerline.  At  a  supersonic  area  ratio 
of  four,  the  material  is  silica  cloth  phenolic  (MX  “2600-96)  also 
at  a  60°  layup  angle  with  respect  to  the  nozzle  centerline.  The 
thermal  analysis  was  performed  using  the  Charring  Material  Abla- 
tion  (CMA)  and  Aerotherm  Chemical  Equilibrium  (ACE)  computer  pro¬ 
grams  (References  4-1,  4-2,  and  4-3).  (The  ACE  computer  program 
is  an  expanded  version  of  the  Equilibrium  Surface  Thermochemistry 
(EST)  computer  program  (Reference  4-4)  .)  The  thermal  analysis  of 
rocket  nozzles  using  these  programs  has  been  documented  in 
References  4-5  and  4-6  and  will  only  be  summarized  here. 


In  brief,  the  ACE  computer  program  was  used  to  generate  the 
required  in-depth  and  surface  thermochemical  information  for  input 
into  the  CMA  computer  program.  This  information  was  generated 
assuming  the  alumina  particles  do  not  penetrate  the  boundary  layer. 
For  carbon  phenolic,  the  surface  information  was  generated  based  on 
diffusion  controlled  thermochemistry.  For  silica  phenolic  both 
diffusion  controlled  thermochemistry  and  mechanical  melt  removal 
were  considered  in  generating  this  information.  The  CMA  program 
calculates  the  transient  one-dimensional  temperature  distributions 
using  the ' ACE  generated  information  together  with  1)  virgin  material 
and  char  thermal  properties,  2)  virgin  material  decomposition  coeffi¬ 
cients  assuming  three  Arrhenius  rate  equations,  3)  incident  radiative 
heat  flux,  and  4)  convective  and  mass  diffusion  coefficients.  Items 
1  and  2 'were  obtained  from  References  4-5  and  4-6.  Item  3  is  satis¬ 
fied  by  specifyipg  the  emissivity  and  radiating  temperature  of  the 
particle  laden  combustion  species,  the  absorptivity  of  the  wall 
«  material,  and  the  view  factor.  The  radiating  temperature  was  assumed 
to  be  the  local  static  temperature,  and  the  emissivity  of  the  cloud 
of  alumina  particles,  using  the  initial  chamber  pressure  of  740  psia, 
was  defined  to  be  0.71  at  the  nozzle  throat  and  0.24  at  a  supersonic 
area  ratio  of  four.  The  absorptivity  (and  emissivity)  used  for  both 
wall  materials  was  0.85,  and  the  view  factor  used  was  1.0.  The 
convective  heat  transfer  coefficient  (Item  4)  was  evaluated  as  a 
function  of  wall  enthalpy  using  the  energy  thickness  method  described 
in  Appendix  A  of  Reference  4-7  to  solve  the  governing  turbulent 
boundary  layer  equations.  Here  again,  the  initial  chamber  pressure 
of  740  psia  Was  used.  Since  previous  experience  (Reference  4-6) 
indicates  that  this  method  of  analysis  consistently  gives  results 
which  cause  the  surface  recession  rate  to  be  overpredicted  by 
approximately  25%,  they  were  multiplied  by  0.75  before  being  input 
to  the  CMA  program.  At  the  nozzle  throat  location,  the  CMA  program 
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modifies  the  initial  convective  heat  transfer  coefficient  according 
to  the  proportionality 


where 


PUCh 


,  1 . 8-0 . 2n 

(-|y)  1-n 


(4-1) 


puCjj  -  convective  heat  transfer  coefficient 
r*  -  nozzle  throat  radius 

n  -  propellant  burning  rate  exponent  (0.625) 


Equation  (4-1)  accounts  for  the  decrease  in  the  convective  heat 
transfer  coefficient  due  to  the  chamber  pressure  decay  resulting 
from  the  nozzle  throat  erosion.  Based  on  the  CMA  analysis  at  the 
nozzle  throat,  the  convective  heat  transfer  coefficient  at  the 
station  defined  by  an  initial  area  ratio  of  four  was  evaluated 
as  a  function  of  time,  accounting  for  the  variation  in  chamber 
pressure  and  local  area  ratio  caused  by  the  nozzle  throat  erosion. 
The  mass  transfer  coefficient  (Item  4)  for  both  locations  was 
evaluated  by  employing  the  relationship. 


CM  = 


2/3 


(4-2) 


where 

CM  -  non-dimensional  mass  transfer  coefficient 
,  CH  -  Stantort  number 

Le  -  Lewis  number 


Using  the  method  of  analysis  described  above,  transient  one-dimensional 
temperature  profiles  were  calculated  at  the  nozzle  throat  (MX  4926) 
and  in  the  exit  cone  at  an  area  ratio  of  four  (MX  2600-96) .  Repre¬ 
sentative  profiles  are  shown  in  Figures  4-1  and  4-2  for  times  corres¬ 
ponding  to  motor  burnout  and  to  60  seconds  following  motor  burnout. 

In  addition,  the  surface  recession,  excluding  the  effect  of  char 
swell,  are  tabulated  in  these  figures  as  a  function  of  motor  burn 
time.  Therefore,  these  figures  provide  information  required  for 
defining  the  thermocouple  locations  and  the  minimum  length  of  the 
thermocoupled  plug.  Since  the  preliminary  design  of  the  plug  had 
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five  thermocouples  and  since  it  is  desired  to  have  at  least  three 
thermocouples  recording  at  all  times  during  the  motor  firing,  the 
depth  of  the  thermocouple  nearest  the  surface  was  taken  to  be  ap¬ 
proximately  the  same  as  the  surface  recession  after  20  seconds  of 
motor  firing.  Therefore,  in  the  error  analysis  presented  here,  the 
disturbance  to  the  temperature  field  caused  by  the  presence  of  the 
thermocouple  is  evaluated  for  the  thermocouple  nearest  the  surface 
being  at  a  depth  of  0.15  inches  for  the  MX  4926  material  and  0.11  for 
the  MX  2600-96  material.  The  unperturbed  transient  temperature  at 
these  locations  plus  at  a  depth  of  0.30  inches  for  the  MX  4926  Eire 
shown  in  Figure  4-3.  These  temperatures  provide  the  baseline  for 
defining  differences  between  the  unperturbed  and  perturbed  tempera¬ 
tures;  thus,  the  errors  presented  are  based  on  the  percent  change  in 
the  variable  of  interest  between  the  value  calculated  from  the  unper¬ 
turbed  temperature  and  the  value  calculated  from  the  perturbed  temp¬ 
erature  . 


Since  the  presence  of  the  thermocouple  perturbs  the  tempera¬ 
ture  field  as  a  result  of  two-dimensional  heat  conduction  ef¬ 
fects,  a  transient  two-dimensional  heat  conduction  analysis  was 
performed.  The  computerized  analysis  procedure  employed  for  this 
analysis  was  the  Axi -Symmetric  Transient  Heating  and  Material  Abla¬ 
tion  (ASTHMA)  computer  program  (Reference  4-8) .  Since  the  ASTHMA 
program  does  not  allow  for  in-depth  charring,  a  set  of  constant 
density  thermal  properties  was  generated  based  on  the  relationship 
between  the  CMA  and  ASTHMA  nodal  energy  balances  Euid  on  the  results 
from  the  CMA  analysis.  These  effective  thermal  properties 
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are  shown  in  Figures  4-4  and  4-5.  Using  these  thermal  properties, 
the  impressed  temperature  and  surface  recession  defined  by  the  CMA 
analysis,  and  assuming  no  radiative  heat  transfer  across  the 
cavity,  ASTHMA  was  run  for  the  nodal  geometries  shown  in  Figures  4-6 
and  4-7.  The  void,  designated  as  the  null  node  (or  nodes),  simulates 
one  half  of  the  thermocouple  hole  cavity  as  the  heat  conduction  is 
symmetrical  about  a  line  which  is  both  perpendicular  to  the  surface 
and  passes  through  the  center  of  the  cylindrical  cavity.  The  y-dimen- 
sion  of  the  nodal  network  was  made  large  enough  so  that  any  disturb¬ 
ance  caused  by  the  null  node  (or  nodes)  would  not  influence  the 
temperature  distribution  at  the  left  hand  face  (y=L) .  The  disturb¬ 
ance  caused  by  the  void  is  then  the  difference  between  the  temperature 
distribution  at  y  =  0  and  at  y  =  L.  These  distributions  are  shown  in 
Figures  4-8  and  4-9  for  the  MX  4926  and  MX  2600-96  plugs,  respectively. 
The  errors  (designated  by  the  solid  lines)  are  plotted  as  a  function 
of  time  in  Figure  4-10.  As  noted  from  this  figure,  the  maximum  AT 
for  the  MX  4926  is  approximately  50°R  which  is  less  than  1%  of  the 
local  temperature.  The  maximum  error  for  MX  2600-96  is  approximately 
+200°R  or  8%  of  the  local  temperature.  The  abrupt  reversal  in  the 
sign  of  the  temperature  difference  occurs  in  the  temperature  range 
in  which  the  thermal  conductivity  goes  through  an  abrupt  increase 
(see  Figure  4-5) .  Due  to  the  large  error  which  occurred  in  the  MX 
2600-96  material,  another  look  was  taken  at  both  the  analysis  pro¬ 
cedure  and  the  results.  This  reinvestigation  showed  that,  particu¬ 
larly  for  t^he  MX  2600-96  material,  radiative  heat  transfer  across 
the  void  was  significant.  The  ASTHMA  analysis  for  the  MX  2600-96 
material  was  modified  to  approximately  account  for  this  radiative 
heat  transfer.  The  error,  AT,  resulting  from  this  analysis  is  also 
shown  in  Figure  4-10,  and  it  is  significantly  less  than  the  one 
obtained  when  radiative  heat  transfer  across  the  void  was  considered 
negligible.  The  maximum  error  is  +25°R  which  results  in  an  accuracy 
of  2^5  percent  for  local  material  temperatures  of  approximately 
1000 °R  and  less  than  -1  percent  for  local  material  temperatures  above 
2500°R. 


’The  only  remaining  temperature  distortion  caused  by  the  plug 
is  due  to  the  bondline  between  the  plug  and  the  remainder  of  the 
material.  To  accurately  assess  the  magnitude  of  this  distortion, 
one  would  require  an  exact  measurement  of  the  bondline  thickness, 
the  charring  characteristics  of  the  bonding  agent,  and  a  detailed 
stress  analysis  which  accurately  depicts  the  deformations  resulting 
from  pressure  and  thermal  stresses.  Since  this  degree  of  sophisti¬ 
cation  is  not  warranted,  an  analysis  was  performed  which  would  yield 
the  maximum  possible  distortion.  This  was  done  by  comparing,  at  the 
nozzle  throat,  the  temperature  profile  resulting  from  an  axisymmetric 
solution  to  the  temperature  profile  resulting  from  a  flat  plate 
solution.  For  the  flat  plate  solution  to'  be  valid,  an  infinite 
thermal  contact  resistance  would  have  to  exist  at  the  thermal  bond¬ 
line.  The  results  of  the  analysis  are  shown  in  Figure  4-11  which 
shows  that  the  maximum  distortion  occurs  in  the  temperature  range 
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SPECIFIC  HEAT,  BTU/LB  °F 


THERMAL  CONDUCTIVITY  xlO 


TEMPERATURE,  °R 


FIGURE  4-5  CONSTANT  DENSITY  THERMAL 

CONDUCTIVITY  WHICH  ACCOUNTS 
FOR  IN-DEPTH  CHARRING 
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TEMPERATURE  DIFFERENCE  (  T),°R 


FIGURE  4-10  TEMPERATURE  DISTORTION  AROUND' 
THERMOCOUPLE  CAVITY 
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from  1000°R  to  2000°R  at  the  later  firing  times.  Since  an  infinite 
resistance  is  certainly  not  realistic  and  since  the  thermal  con¬ 
ductivity  of  the  phenolic  resin  and  the  bonding  agent  are  similar, 
the  temperature  distortion  would  be  expected  to  be  less  than  1%. 


4. 1.1. 1.2  Response  Lag.  The  configuration  for  the  thermocouple 
wire  and  junction  in  the  plug  assembly  is  shown  in  Figure  4-12.  The 
^thermal  response  of  the  thermocouple  junction  ii  governed  by  the  net 
radiative  heat  transfer  tp  the  wire  and  by  the  conductive  heat 
transfer  from  the  material  to  the  junction  at  the  point  of  contact. 
The  maximum  thermal  lag  of  the  thermocouple  would  occur  when  there  is 
no  contact  between  the  junction  and  the  ablative  material.  This  is 
not  the  desired  configuration  but  is  investigated  here  to  determine 
whether  or  not  it  is  tolerable.  • 


The  net  radiative  transfer  to  the  thermocouple  wire  can  be 
obtained  assuming  two  concentric  cylinders  and  is  given  by 


oAi<t24-  tx4) 

q  =  l/ex  +  A1/A2U/e2  -  1) 


(4-3) 


where 

a  -  Stefan-Boltzmann  constant 
A  -  cylindrical  surface  area 
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'"e  -  emissivity 
T  -  temperature 
q  -  net  radiative  heat  flux 

and  subscripts  1  and  2  correspond  to  the  inner  and  outer  cylinder, 
respectively.  Equation  (4-3)  can  be  put  in  the  form 

’  ’  °Al£eff(Tl4  -  T24> 


(4-4a) 


where 


'eff 


1/ej^  +  a1/a2(1A2  -  1) 


(4-4b) 


Employing  the  assumption  that  the  thermal  content  of  the  ablative 
material  Adjacent  to  the  thermocouple  junction  is  large  compared  to 
tne  heat  flux  to  the  junction,  the  thermal  model  for  determining 
the  lag  of  the  thermocouple  can  be  put  in  the  form  shown  in  Figure 
4-12,  The  above  assumption  permits  the  driving  potential  for  the 
heatf  transfer  by  conduction  from  the  ablation  material  to  the 
thermocouple  junction  to  be  an  infinite  source  at  the  local  material 
temperature.  This  temperature  is  taken  to  be  the  unperturbed  material 
temperature  as  given  in  Figure  4-1.  Also  shown  in  the  thermal  model 
is  an  infinite  sink  which  is  held  at  the  initial  temperature  (ap¬ 
proximately  room  temperature) .  The  effect  of  this  sink  is  to  maxi¬ 
mize  the  conduction  losses.  The  radiative  flux  is  applied  along 
the  horizontal  portion  of  the  thermocouple  wire.  Using  the  thermal 


/ 
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model  of  Figure  4-12,  two  solutions  were  obtained  for  the  thermocouple 
nearest  the  surface;  namely  1) ^infinite  thermal  contact  resistance 
with  radiation  to  the  thermocouple  wire  and  2)  both  radiative  and 
conductive  heating  of  the  thermocouple  wire  assuming  zero  contact 
resistance.  The  thermal  properties  used  in  the  analysis  sire  listed 
in  Table  4-1  and  the  source  temperatures  at  the  two  different 
locations  are  from  Figure  4-3.  The  results  of  this  analysis  sire 
given  in  Figure  4-13.  For  the  case  of  both  conductive  and  radiative 
heating,  the  thermal  contact  surface  area  was  taken  as  17  percent  of 
the  surface  area  o£  the  bead.  From  the  results  of  Figure  4-13,  one 
concludes  that  1)  at  temperatures  below  3000°R,  contact  must  be  main¬ 
tained  between  the  material  and  the  thermocouple  junction,  and  2)  the 


thermal  conduction  error,  which  is  less  than  2.5  percent  at  tempera¬ 
tures  above  3000°R,  is  acceptable. 


TABLE  4-1 


THERMAL  PROPERTIES  OF  TUNGSTEN 


Density  lb/ft3 

1185.0 

Specific  Heat,  Btu/lb°F 

@  527or 

0.33 

@  2290°R 

0.41 

@  4100°R  r 

0.47 

|  Thermal  Conductivity,  Btu/hr  ft°F 

€  527°R 

96.5 

$  1390°R 

70.0 

@  2290°R 

65.4 

@  4100OR 

60.4 

4. 1.1. 2  Accuracy  of  Thermocouple  Output 

The  first  problem  of  the  measurement  system  concept  was 
the  distortion  of  the  temperature  field  which  was  discussed  in  the 
previous  section  (Section  4. 1.1.1).  The  second  problem  is  accurately 
defining  the  temperature  of  the  thermocouple  junction  from  the  emf 
output.  Factors  which  could  cause  errors  in  the  emf  measurement  are 
electrical  shorting  and  thermocouple  (including  lead  wire)  calibration 
The  first  item  is  associated  with  shunting  across  two  points  at 
diflerent  electrical  potential  resulting  in  secondary  electrical 
circuits  within  the  primary  thermocouple  electrical  circuit.  The 
second  item  is  associated  with  the  guaranteed  accuracy  of  the 
thermocouple  and  lead  wires  and  with  the  tungsten-carbon  reaction 
at  high  temperatures  which  can  cause  a  shift  in  the  emf-temperaturc 
characteristics  of  the  thermocouple  junction  for  certain  environ¬ 
mental  conditions ,  These  potential  sources  of  error  are  considered 
in  the  following  subsections. 
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4 . 1 . 1 . 2 . 1  Electrical  Shorting.  The  voltage  measured  in  a  thermo¬ 
couple  circuit  is  due  to  the  Seebeck  emf  which  is  a  combination  of 
the  Peltier  and  Thomson  emf's.  The  Peltier  emf  is  generated  at  the 
thermocouple  junction  and  the  Thomson  emf  can  be  visualized  as  being 
generated  in  the  thermal  gradient  region  (References  4-9  and  4-10) 
of  each  of  the  homogeneous  wires  making  up  the  thermocouple  circuit. 
These  thermoelectric  emf’s  are  evaluated  b£  using  the  thermocouple 
equations  in  the  form  (Reference  4-11)  j* 


E 


A,B 


d£ 


A,B 


dT 

d2E 


-T 


dT 


A,B 

T~ 


"a.B^T  ~  ^A ,B^Tq  +  T^  (oft 


0B)dT 


(4-5a) 

(4-5b) 


(4-5c) 


where 


E  -  Seebeck  voltage  for  metals  A  and  B 
'  for  a  temperature  difference  of  T  -  Tq 

T  -  temperature 
*  * 

S. ,S_  -  entropy  tr«msport  coefficient  for  metals 

A  and  B  respectively 


it,  B  -  Peltier  coefficient  for  the  combination 
'  of  metals  A  and  B 

0ft,  0g  -  Thomson  coefficients  for  metals  A  "and  B 
respectively 


If  the  relationship  between  g  311(1  T  is  known,  then  the  magnitude’ 
of  the  Peltier  and  Thomson  emf  *s  can  be  determined.  For  W5Re/W26Re, 
this  relationship  is  (Reference  4-12)  i. 


E  =  0.293341  +  9.01521  x  10_3  T  +  2.70124  x  10_5  T2 
-3.15783  X  10~8  T3  +  1.78838  x  10-11  T4 
-5.04944  x  10_15  T5  +  5.52473  x  10~19  T6  (4-6) 


The  emfs  resulting  from  this  equation  are  shown  in  Figure  4-14. 


An  electrical  circuit  for  one  W5Re/W26Re  thermocouple  in  the 
plug  assembly  is  shown  in  Figure  4-15.  The  W5Re/W26Re *was  chosen 
for  analysis  because  it  would  be  exposed  to  the  most  severe  tempera¬ 
ture  and  temperature  gradient.  The  Peltier  and  Thomson  emf's  are 
distributed  at  various  distinct  locations  with  only  the  Peltier  emf 
in  series  with  each  shunt  resistance  as  the  shunt  leg  is  assumed  to 
be  along  an  isothermal  plane.  The  number  of  locations  is  determined 
by  the  accuracy  desired  in  the  analysis.  It  is  obvious  from  the 
circuit  that  to  obtain  an  accurate  measurement  of  the  thermocouple 
junction  emf,  shunting  across  the  two  legs  of  the  thermocouple  wire 
must  be  negligible.  This  means  that  the  electrical  resistance  in 
each  of  the  shunt  legs  must  be  large  compared  to  the  electrical 
resistance  in  the  remainder  of  the  electrical  circuit.  The  electrical 
resistance  in  the  shunt  leg  is  made  large  by  encapsulating  the  lead 
wire  in  beryllia  (BeO)  or  alumina  (AI2O3)  insulation.  Beryllia 
insulation  was  chosen  because  of 

•  its  higher  melt  temperature  than  alumina  (4600°F 
compared  to  3700°F) 

•  its  compatibility  with  tungsten  rhenium  wire  at 

high  temperature  as  compared  to  thoria,  for  ex¬ 
ample,  which  affects  the  calibration  curve  of 
tungsten-rhenium  thermocouples.  \x 

'  \ 

*  \ 

With  the  electrical  insulator  chosen,  all  the  materials  used  in 
the  thermocouple  electrical  circuit  have  been  defined.  Before  pro¬ 
ceeding  with  the  analysis,  however,  the  circuit  was  simplified  by 
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combining  the  Peltier  and  Thomson  emf's  into  a  distributed  Seebeck 
emf  as  shown  in  Figure  4-16.  This  simplification  is  conservative 
in  that  it  will  slightly  magnify  the  shunt  effect  due  to  the 
Seebeck  emf  being  larger  than  the  Peltier  emf.  A  representative 
temperature  distribution  which  was  used  in  the  analysis  is  also  shown 
in  this  figure. 

t 

The  equations  describing  the  electrical  circyit  shown  in 
Figure  4-16  were  solved,  and  .this  solution  is  described  in  consider¬ 
able  detail  in  Reference  4-11.  In  summary,  it  showed  that  due  to 
the  Rjj  resistances  being  a kf actor  of  10 11  to  10*  larger  than  the 
other  resistances,  the  shunting  effect  is  negligible  and  the  measured 
emf  will  indeed  be  the  emf  generated  at  the  thermocouple  junction. 

At  temperatures  above  the  melt  temperature  of  beryllia  (greater  than 
5000°F)  the; shunt  resistance  would  still  be  expected  to  be  adequate 
(Reference  4-12)  provided  the  liquid  beryllia  remains  around  the  thermo¬ 
couple  wire. 
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ELECTRICAL  CIRCUIT 


4. 1.1. 2. 2  Thermocouple  Calibration.  The  remaining  inaccuracies  in 
the  thermocouple  circuit  analysis  are  the  possible  variation  in  the 
emf- temperature  relationship  duo  to  the  tungsten-carbon  reaction  at 
high  temperatures  and  the  guaranteed  calibration  accuracy  of  the 
thermocouple  and  lead  wire  as  supplied  by  the  manufacturer. 


A  literature  survey  was  performed  to  define  whether  or  not  the 
tungsten-carbon  reaction  for  short  durations  at  high  temperature 
would  affect  the  calibration  of  the  thermocouple.  Spooner  and  Hall 
(Reference  4-13)  performed  stability  tests  for  W/W26Re  thermocouples 
in  a  carbon  vapor  atmosphere  for  durations  up  to  30  hours  and  found 
a  maximum  calibration  change  of  approximately  +30°C.  These  reported 
data  are  for  a  different  environment  than  the  one  expected  to  be 
encountered  in  the  thermocouple  plug  assembly,  carbon  vapor  at  2100°C 
to  2500°C  for  30  hours  compared  to  direct  carbon  contact  at  tempera¬ 
tures  from  530°F  to  5500°F  within  60  seconds.  Further  investigation. 
Reference  4-14 r  resulted  in  the  data  shown  in  Figure  4-17.  These 
data  are  .for  a  W/W26Re  thermocouple,  0.020"  wire  diameter,  in  direct 
contact  with  a  carbon  heater  at  3182PF.  The  deviation  is  the  difference 
between  the  readings  from  a  thermocouple  and  an  optical  pyrometer. 


\ 
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-  If  the  differences  in  the  deviations  are  attributed  to  the  tungsten- 
carbon  Reaction,  then  the  maximum  error  is  0°C  for  an  exposure  of 
0.2  hours,  3°C  for  0.3  hours,  and  7°C  for  0.5  hours.  Since  the  de¬ 
viation  isfla  function  of  the  amount  of  material  affected?,  the  vojume 
of  the  thermocouple  junction  would  be  expected  to  influence  the 
deviation,  the  smaller  the  volume,  the  larger  the  deviation. 

Assuming  the  difference  in  deviation  is  inversely  proportional  to 
the  diameter  squared,  of  the  thermocouple  wire,  then  this  difference 
-is  only  48°C  after  0.3  hours  (18  minutes).  Based  on  these  data, 
the  effect  of  the  tungsten-carbon  reactiqn  is  negligible  for  temp¬ 
eratures  up  to  3182°F  (36420R) .  From  Reference  4-15,  these  data 
are  felt  tp  be  representative  for  temperatures  up  to  2500°C  (5000?R). 
at  which  temperature  carbon  and  tungsten  form  a  eutectic.  However, 
due  to  the  short  exposure  times  at  temperatures  above »5000°R,  less 
than  10  seconds,*  the  effect  of  the  carbon- tungsten  reaction  is  ex¬ 
pected  to  be  negligible  even  at  these  very  high  temperatures.  The  ( 
guaranteed  accuracy  (Reference  4-16)  of  the  tungstqn,  rhenium  wire 
and  the  lfead  wires  are  +1%  and  .+0.5%  respectively.  Since  these  are 
,  random  errors,  they  can  be  combined  (  by .root  mean  square)  to  give 
an  accuracy  of  +1.23%.  *' 


4.1. 1.3  -Accuracy  of  Thermocouple  Location 

%  •  i  j* 

The  location  of  the  thermocouple  bead  within  the  heat  flux 
sensor  is  an  important  parameter  which  must  be  accurately  .defined 
to  obtain  useful  thermal -response  data.  For  the  plugs;  fabricated 
in y  this  program1,  the  position  of  the  thermocouples  relative  to  the 
pre-installation  surface  were  found  by  measurement  of  x-'ray  photo-  -- 
graphs.  As  described  in  Section  4. 1.2. 3,  the  plugs  were  ^Installed 
in  the  nozzle  with  some  excess  ping  material  protruding  abov^  the 
flame  surface.  This  protuberance  was  machined  flush  with  the!  sur¬ 
rounding  area.  The  amount  machined  was  measured  and  subtracted;: 
from  the  x-ray  measurements  to  obtain  the  pre-fire^  thermocouple 
depths  below  the  flame  surface.  The  accuracy  analysis  presented 
here  considers  both  the  accuracy  of  the  x-ray  measurement  techniques 
and  the  accuracy  of  measuring  the  amount,  removed  during  the  final 
machining  procedure;  X-ray  measurement  uncertainties  are  discussed 
first,- followed  by  ^  brief  discussion  of  the  final  machining  uncer¬ 
tainty.  Finally,  these  uncertainties  are  combined  and  then  con-i  ~  - — 
N  Verted  to  an  uncertainty  in  the  measured  temperature .  / 


A  schematic  representation  of^the  heat  flux  senspr  positioned 
for  X-ray  exposure  is  given  in  Figure  4-lb.  The  .top  surface  of  ;the 
sensor  is  defined  on  the  X-ray  negative  by  butting  the  sensor 
against  a  steel  slug.  The  position  of  each  bead  in  the  material 
relative  to  the  defined' surface  can  therefore  be  determined-  A 
A  slight  correction  for  parallax  from  the  plane  of  the  thermocouples' 
to  the  film  plane  must  be  made  because  the  incident  X-ray  beams  are 
not  parallel.  For  the  case  where  the  sensor  is  placed  exactly  as 
shown  in  Figure  4-18,  the  thermocouple  depths  are  evaluated  as  follows 

;  \ 
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FIGURE  4-18  TOTAL  HEAT  FLUX  SENSOR  CORRECTLY 
POSITIONED  IN  X-RAY  BEAM 
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(4-7) 


In  order  to  evaluate  the  uncertainty  in  the  calculated  thermo¬ 
couple  depths,  the  method  of  Reference  4-17  was  utilized  in  this 
analysis.  Briefly,  the  method  considers  the  mathematical  propaga¬ 
tion  of  uncertaintes  in  primary  quantities  into  the  uncertainty  in 
the  calculated  result  (the  primary  quantities  being  those  variables 
which  must  be  directly  or  indirectly  measured) .  Two  sources  of 
error  exist  in  the  X-ray  technique: 

•  -uncertainties  in  the  measurements  of  the  bead  locations 

on  the  X-ray  negative 

•  uncertainties  in  the  positioning  of  the  sensor  slug 

in  the  X-ray  beam 

Figure  4-19  gives  an  exaggerated  schematic  of  a  sensor  slug  incor¬ 
rectly  placed  for  the  X-ray  photograph.  The  primary  geometrical 
quantities  which  need  to  be  measured  to  determine  the  thermocouple 
depths  in  such  a  set  up  are: 

•  distance  from  the  X-ray  source  to  the  film  plane,  Y 

•  distance  from  the  plane  of  the  thermocouples  to  the  film 

plane,  Y 

•  angular  relation  of  the  thermocouple  plane  to  the  film 

plane,  $ 

•  angular  displacement  of  the  sensor  surface  from  the  verticle 

ray,  6 


The  equation  to  evaluate  the  thermocouple  depths  from  the 
above  quantities  was  easily  derived  from  geometric  considerations 


w  cos ($—6) _ 

2(1  -  1  -v"/2  ) 


(4-8) 


Notice  that  Equation  (4-7)  is  the  same  as  Equation  (4-8)  except 
for  the  additional  term  containing  the  <}>  and  6  angles.  Since  4> 
and  8  are  very  close  to  90°  and  0°,  respectively,  the  last  term  is 
generally  very  close  to  zero  and,  for  the  purpose  of  calculating 
6  ,  may  be  ignored.  However,  to  determine  the  propagation  of 

angular  uncertainties  to  the  thermocouple  depth  result,  the  term 
must  be  considered. 


Based  on  the  method  of  Reference  4-17  and  Equation  (4-8) 
the  following  equation  for  the  uncertainty  in  the  thermocouple 
position  was  derived 


4-25 


s 


r 


+  to  2)  + 

X 

o 


( - 

2(1- 


w 

y  -  w/2 


Y 


•) 


2.2  2, 

<“*  +  V 


+ 


V 


"]■ 


(4-9) 


where  to  is  the  uncertainty  in  the  n*-*1  variable  and  where  <j>  and 
0  have  Been  evaluated  at  90°  and  0°,  respectively. 


Table  4-2  presents  the  results  of  applying  Equation  (4-9)  to 
thermocouples  at  approximately  0.1  and  0.6  inches  below  the  surface. 
The  parameter  values  and  uncertainties  given  in  the  table  are  based 
on  past  experience  with  the  X-r&y  measurement  technique.  Also  shown 
in  Table  4-“2  for  each  primary  quantity  is  the  thermocouple  location 
uncertainty  which  would  be  calculated  if  the  uncertainties  of  the 
other  primary  quantities  were  neglected.  The  resultant  uncertainty 
is  the  square  root  of  the  sum  of  the  squares  of  the  individual  un¬ 
certainties.  Table  4-2  shows  that  the  location  uncertainty  is  quite 
insensitive  to  the  thermocouple  depth  below  the  pre-installation 
surface  (±0.0016  to  ±0.0017  inch).  The  most  significant  sources  Of 
error  are  due  to  the  X-ray  negative  reading  uncertainty.  The 
angular  orientation  uncertainties  also  contribute  slightly  to  the  £ 
overall  location  uncertainty. 

After  installation  and  final  potting  of  the  thermocouple  plug 
within  the  nozzle,  the  amount  of  plug  protruding  from  the  flame 
surface  is  measured  either  by  a  vernier  calipers  or  from  dial  set¬ 
tings  during  machining.  Least  indicator  accuracy  in  either  case 
is  ±0.001  inches.  This  random  uncertainty  combines  as  the  square 
root  of  the  sum  of  the  squares  with  the  X-ray  dimensional  measure¬ 
ment  to  give  a  cumulative  location  uncertainty  of  ±0.002  inches. 

The  relation  of  thermocouple  location  uncertainty  to  an  equiva¬ 
lent  temperature  uncertainty  is  defined  by  the  slope  of  the  temper¬ 
ature  profile.  For  a  typical  MX  4926  carbon  phenolic  response  such 
as  shown  in  Figure  4-1,  a  ±0.002  inch  location  uncertainty  would 
present  a  ±30°R  temperature  uncertainty  or  about  ±1.2  percent  for 
local  temperatures  above  2500°R.  For  MX  2600-96,  however,  the 
slope  of  the  temperature  profile  is  larger  so  that  a  ±0.002  inch 
uncertainty  would  represent  an  uncertainty  of  ±40°R  or  ±1.5  percent 
for  local  temperatures  above  2500°R. 


4. 1.1. 4  Overall  Temperature  Measurement  Accuracy 

In  Sections  4. 1.1.1  through  4. 1.1. 3,  accuracies  were 
defined  for  those  items  which  would  cause  the  temperature  obtained 
from  the  measured  thermocouple  emf  and  for  the  measured  thermocouple 
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location  to  be  different  than  the  unperturbed  temperature  which 
would  be  obtained  from  an  ideal  precisely  located  thermocouple. 

In  this  section,  these  inaccuracies  will  be  combined  to  give  an 
overall  system  accuracy  for  the  in-depth  temperature  measurements . 
Since  the  temperatures  above  2500°R  are  the  ones  of  primary  signi¬ 
ficance  in  defining  the  thermal  heat  conduction  into  the  ablative 
material,  the  accuracies  presented  are  limited  to  this  temperature 
range . 


Table  4-3  presents  the  accuracies  defined  in  Sections  4. 1.1.1 
through  4. 1.1. 3  for  a  local  material  temperature  of  approximately 
2500°R.  If  the  accuracies  preceded  by  a  plus  or  minus  sign  (but 
not  both)  are  considered  to  result  from  fixed  errors  and  all  others 
to  result  from  random  errors ,  then  the  desired  overall  system  accuracy 
can  be  obtained  from 


6 

system 


l  l  fixed  errors 


random  errors 


(4-19) 


Substituting  the  numbers  listed  in  Table  4-3  into  Equation  4-19, 
the  overall  accuracy  for  a  local  material  temperature  of  about 
2500°R  is  +0.2  percent  and  -3.2  percent  for  MX  4926  and  -1.6  percent 
and  -5.4  percent  for  MX  2600-96.  For  higher  material  temperatures, 
the  magnitude  of  all  the  fixed  errors  in  terms  of  temperature  dif¬ 
ference  decrease  with  increasing  temperature;  thus  the  accuracies  on 
a  percentage  basis  will  be  better  at  the  higher  material  temperatures. 
Assuming  the  magnitude  of  the  fixed  errors  (and  the  percentage  of 
the  random  errors)  remains  the  same,  the  accuracies  considering  only 
the  effect  of  the  higher  material  temperature  are  +0.8  percent  and 
-2.7  percent  for  MX  4926  and  +0.3  percent  and  4.1  percent  for  MX 
2600-96,  based  on  a  material  temperature  of  4000°R. 


TABLE  4-3 


ACCURACIES  IN  MEASURED  TEMPERATURES 
FOR  A  LOCAL  MATERIAL  TEMPERATURE 
OF  ABOUT  2500°R 


Item 

Section 

Error  (Percent) 

MX  4926 

MX  2600-96 

Distortion  of 

Temperature  Field 

4. 1.1. 1.1 

+  1 

-1 

Response  Lag  of 

Thermocouple 

4. 1.1. 1.2 

-2.5 

-2.5 

Guaranteed  Accuracy 

» 

of  Thermocouple 

and  Lead  Wires 

4. 1.1. 2. 2 

+1.23 

+  1.23 

Thermocouple  Depth 

4. 1.1.3 

+1.2 

+1.5 

Overall 

— 

+0.2 

-1.6 

-3.2 

-5.4 
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The  accuracies  stated  in  this  section  are,  with  the  exception 
of  the  thermocouple  and  lead  wires,  based  on  analysis  using  con¬ 
servative  assumptions.  If  these  assumptions  are  correct,  which  there 
is  no  reason  to  doubt  at  the  present  time,  and  if  the  problems  con¬ 
sidered  cure  the  only,  significant  ones,  then  the  above  stated  accuracies 
are  realistic  and  attainable . 


4.1.2  Design  and  Fabrication  Analysis 

The  design  of  the  total  heat  flux  measurement  system  combines 
the  practical  requirements  of  fabrication  and  assembly  and  the 
accuracy  requirements  discussed  in  the  previous  section.  The 
design  and  fabrication  techniques  which  satisfied  these  requirements 
are  discussed  in  the  following  sections.  Section  4. 1.2.1  presents 
the  thermocouple  locations,  types,  and  installation  scheme.  Section 

4. 1.2. 2  presents  the  overall  measurement  system  configuration  andf 
Section  4. 1.2. 3  describes  its  installation.  Finally,  Section  4.l'.2.4 
discusses  the  recording  instrumentation  required  for  the  system 
measurements. 


4 . 1 . 2 . 1  Thermocouples 

The  considerations  related  to  the  thermocouple  installation 
were  the  location  and  number  of  thermocouples  in  the  measurement 
system  plug,  the  types  of  thermocouples  required,  the  configuration 
of  the  hole  which  accepts  the  thermocouples,  and  the  method  of  lead¬ 
ing  the  thermocouple  wires  down  the  side  wall  of  the  plug.  The 
measurement  system  design  in  these  areas  is  presented  in  the  follow¬ 
ing  paragraphs. 


4. 1.2. 1.1  Thermocouple  Locations  and  Types.  The  definition  of  the 
location  and  number  of  thermocouples  in  the  measurement  system  plugs 
was  based  on  Jthe  anticipated  response  of  the  material  in  which  it  is 
installed.  This  anticipated  response  was  calculated  using  the 
Aerotherm  Chemical  Equilibrium  (ACE)  and  Charring  Material  Ablation 
(CMA)  computer  programs  for  the  Phase  III  char  motor  nozzle,  and  was 
estimated  based  on  past  experience  and  calculations  for  the  Phase  II 
evaluation  tests.  The  calculated  response  of  the  char  motor  nozzle 
at  the  throat  (MX  4926  carbon  phenolic)  and  at  A/A*  =  4.0  (MX  2600-96 
silica  phenolic)*  is  presented  in  Figures  4-20  and  4-21  in  terms  of 
the  in-depth  temperature  distributions  at  three  times  through  the 
firing.  The  selected  thermocouple  locations,  based  on  this  anticipated 
response,  are  indicated  in  the  figures.  The  five  thermocouples  for 


♦The  two  other  measurement  system  locations  are  A/A*  =  -2.0  (subsonic) 
and  2.0,  both  being  MX  4926  carbon  phenolic. 
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each  material  and  location  provide  an  accurate  definition  of  the  in- 
depth  response  throughout  the  60  second  firing.  The  thermocouple 
pattern  of  Figure  4-20  is  also  proper  for  the  other  two  measurement 
system  locations,  MX  4926  carbon  phenolic  at  A/A*  =  -2.0  (subsonic) 
and  2.0. 


The  thermocouple  pattern  for  the  Phase  II  evaluation  tests  is 
presented  in  the  table  below. 


Depth  Below  Surface 
•(in) _ 


0.050 

0.125 

0.200 

0.300 

0.400 

This  pattern  was  used  for  both  the  MX  4926  carbon  phenolic  and 
the  MX  2600-96  silica  phenolic  materials,  both  of  which  were 
subjected  to  similar  test  conditions.  (See  Section  6  for  the  com¬ 
plete  outline  of  this  program.)  Since  the  surface  recession  was 
somewhat  less  than  that  anticipated  in  the  char  motor,  the  thermo¬ 
couples  are  closer  to  the  surface. 


The  types  of  thermocouples  to  be  used  in  the  plugs  were 
selected  based  on  the  maximum  expected  temperature  at  the  particular 
thermocouple  locations,  the  thermocouple  response  characteristics, 
and  the  fabrication  and  handling  characteristics  of  the  thermocouples. 
The  tungsten  rhenium  thermocouple  systems  provide  the  highest  tempera¬ 
ture  capability.  Standard  calibrations  are  available  to  4200°F 
for  both  tungsten/tungsten  26%  rhenium  (W/W26Re)  and  tungsten  5% 
rhenium/tungsten  26%  rhenium  (W5Re/W26Re)  as  shown  in  Figure  4-22.  ' 
Calibration  results  to  higher  temperatures  are  also  available;  results 
from  Reference  4-12  to  5400°F  for  W5Re/W26Re  are  also  shown  in 
Figure  4-22.  As  seen  from  the  figure,  W/W26Re  has  a  greater  output 
sensitivity  (millivolts  per  degree)  at  high  temperature  and  is  there¬ 
fore  preferable  from  this  standpoint.  Based  on  Aerotherm  experience 
and  the  experience  of  others  (e.g.  References  4-12  and  4-16),  pure 
tungsten  presents  problems  in  handling,  however.  Pure  tungsten  is 
very  brittle  and  the  thermocouple  junction  formed  is  very  delicate. 

The  W5Re/W26Re  system  reduces  these  problems  considerably  although 
care  in  handling  is  required  with  this  system  as  well.  On  the  basis 
of  the  above,  W5Re/W26Re  was  chosen  as  the  thermocouple  type  for  all 
locations  where  temperatures  above  2500°F  are  expected. 
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For  temperatures  below  2500°F,  the  Chromel/Alumel  (C/A)  system 
was  chosen.  This  thermocouple  type  combines  high  output  sensitivity 
(Figure  4-22)  and  moderately  high  temperature  capability  and  is  also 
rugged.  The  other  possible  candidates  for  this  temperature  range 
are  the  platinum  rhodium  systems;  however,  they  exhibit  a  quite  low 
output  sensitivity. 


The  thermocouple  wire  to  be  used  in  all  measurement  system  plugs 
was  calibrated  by  the  manufacturer  (Hoskins  Manufacturing  Co.),  this 
calibration  being  directly  traceable  to  the  National  Bureati  of 
Standards . 


Based  on  the  selected  thermocouple  locations  and  on  the  tempera¬ 
ture  capabilities  of  the  thermocouple  types,  the  locations  and  cor¬ 
responding  types  of  thermocouples  selected  for  the  Phase  III  chaf 
motor  nozzle  and  the  Phase  II  test  sections  are  presented  in  Table 
4-4.  Because  of  the  relatively  high  surface  recession  anticipated  in 
the  carbon  phenolic  part  of  the  char  motor  nozzle,  four  of  the  five 
thermocouples  at  the  three  measurement  system  locations  are  W5Re/W26Re 


/ 
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TABLE  4-4. 


THERMOCOUPLE  DEPTHS  AMD  TYPES  TO  BE  USED  IM 
MEASUREMENT  SYSTEM  PLUGS 


Char  Kotor 


MX  4926  Carbon  Phegol ic 
-  -  •  0,  '  - 


NX  2600-96  Silica  Phenolic 
A/A#  -  4.0 


Depth* 

( inch) 

Thermocouple 

Type 

Depth 

(inch) 

Thermocouple 

Type 

Depth 

( inch) 

* 

Thermocouple 

YYPe 

0.125 

W5ReA26R* 

ns 

M5Re/W26Re 

o.oso 

WSRe/V26Re 

0.250 

SUM 

0.125 

0.400 

mBSM 

0.200 

0.600 

0.450 

C/A 

0.300 

c/A 

0.800 

C/A 

0.600 

1 

0.400 

, 

Phaae  II  T  ..t  Section* 


HX  4926  Carbon  phanolic  and 
MX  2600-96  Silica  phenolic 
1.0 


a)  Depth  below  original  surface 

b)  Tung* ten  5%  rhenltaa/tungsten  26%  rheniue  or  Chroaael/Alunel 


In  the  other  cases,  it  is  necessary  that  only  three  be  W5Re/W26Re. 
These  requirements  are  .apparent  for  the  char  motor  nozzle  from 
Figures  4-20  and  4-21. 


4. 1.2. 1.2  Thermocouple  Installation.  The  configuration  of  the  holes 
which  accept  the  thermocouples  is  shown  in  the  sketch  below.  The  step- 
ped-hole  scheme  provides  intimate  contact  between  the  thermocouple  and 
material  and  provides  the  means  of  precisely  locating  the  thermocouple 
bead  on  the  plug  centerline.  The  0.010-inch  diameter  hole  size  was 
dictated  by  the  type  of  materials  of  interest.  The  cloth  phenolics 
are  very  hard  and  also  have  macroscopic  (relative  to  0.010  inch)  vari¬ 
ations  in  density  and  structure.  This  combination  represents  a  sig¬ 
nificant  hazard  to  small  drills  both  in  terms  of  dulling  and  breaking. 
Also,  once  broken,  the  removal  of  a  drill  is  usually  impossible.  The 
0.010-inch  hole  was  therefore  found  to  be  the  minimum  practical  size 
for  successful  drilling  of  these  materials. 
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The  0.010-inch  hole  diameter  dictated  the  selection  of 
0.005-inch  diameter  thermocouple  wire.  This  is  a  standard  and 
readily  available  size  and  is  also  a  convenient  size  to  handle  and 
install.  The  thermocouple  bead  size  which  results  is  about  0.015  inch 
in  diameter;  this  dictated  the  0.020-inch  diameter  hole  size  which 
accepts  the  bead. 


Because  the  drills  tend  to  follow  the  ply  direction  in  cloth 
reinforced  materials,  drilling  at  any  angle  not  parallel  to  the  plies 
greatly  increases  the  chances  of  drill  breakage.  All  thermocouple 
holes  were  therefore  parallel  to  each  other  in  a  direction  also 
parallel  to  the  plies  as  shown  in  Figure  3-2 .  in  order  to  keep  the 
thermocouple  locations  as  close  to  the  plug  centerline  is  practical, 
the  three  closest-to-the-surface  thermocouples  are  led  down  the  side 
walls  through  the  insulative  sleeving  and  the  two  other  thermocouples 
are  potted  into  the  groove  formed  by  the  two  adjacent  insulators.  The 
insulative  sleeving  used  is  beryllium  oxide  as  discussed  in  Section  4.1.1 
The  potting  material  selected  is  a  boron  nitride  slurry  which  cures 
to  a  hard,  high  temperature,  non-conducting  material.  The  sleeving 
diameter  selected  is  0.029  +.002  inch  because  it  is  a  convenient, 
commercially  available  size. 


4. 1.1. 2  System  Configuration 

The  overall  measurement  system  plug  configuration  was 
presented  previously  in  Figure  3-2.  a  more  detailed  description  is 
presented  in  the  following  paragraphs . 
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The  basic  plug  diameter  is  1/2  inch.  This  is  a  practical  upper 
limit  for  convenience  in  thermocouple  hole  drilling  and  a  practical 
lower  limit  in  terms  of  the  thermocouple  wire  conduction  loss  dis¬ 
cussed  in  Section  4. 1.1. 1.2.  The  overall  plug  length  for  the  char  motor 
nozzle  was  approximately  1/4  inch  less  them  the  thickness  of  the 
part  in  which  it  was  installed.  This  provides  sufficient  clearance 
below  the  back  wall  of  the  part  for  the  lead  wire  connections  and 
allows  adequate  potting  for  plug  retention. 


The  measurement  system  plugs  were  fabricated  from  billets  of  the 
two  materials,  MX  4926  carbon  phenolic  and  MX  2600-96  silica  phenolic, 
such  that  the  layup  angles  matched  those  for  the  locations  at  which 
they  were  installed.  The  surface  and  in-depth  response  of  the  plug 
therefore  matched  that  of  the  surrounding  material,  as  required.  The 
surface  contour  of  the  plugs  was  machined  to  match  the  contours  at 
the  respective  installation  locations. 


The  thermocouple  lead  wire  junctions  at  the  bottom  of  the  plug 
were  formed  on  pins  mounted  in  the  plug  as  shown  in  Figure  4-23.  . 

This  technique  eliminated  any  strain  on  the  thermocouple  wires  due  to 
handling  the  lead  wires.  The  pin  was  a  short  section  of  large  diameter 
thermocouple  lead  wire  of  the  same  material  as  the  thermocouple  lead 
itself.  This  eliminated  any  errors  due  to  a  possible  temperature 
gradient  between  the  thermocouple  connection  point  and  the  lead  wire 
connection  point.  The  end  of  the  connection  pin  that  was  mounted  in 
the  plug  was  dipped  in  epoxy  resin  before  installation  so  that  it 


FIGURE  4-23  SCHEMATIC  OF  THERMOCOUPLE  LEAD 
WIRE  CONNECTION  SCHEME 

*•  I4lb 
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was  insulated  from  the  plug  once  installed.  The  pins  were  also  held 
in  place  with  epoxy  resin.  The  thermocouple  lead  and  the  lead  wire 
were  spot  welded  to  the  pins. 


4. 1.2. 3  Installation 

The  measurement  system  plug  was  press-fit  into  its  mating 
hole  in  the  ablative  part  or  test  section.  The  back  of  the  plug  was 
slotted  for  a  screwdriver  as  shown  in  Figure  4-23  so  that  rotational 
and  depth  adjustments  could  be  made.  The  connecting  pin  pattern  pro¬ 
vided  ample  room  for  a  screwdriver  without  danger  to  the  thermocouple 
and  leads . 


The  hole  in  the  ablative  part  and  the  thermocoupled  plug  were 
match  machined  in  order  to  insure  a  snug  fit.  For  the  Phase  III  char 
motor  firings,  this  hole  was  counter-bored  as  shown  in  the  sketch  below 


to  a  depth  which  allowed  the  plug  to  protrude  above  the  internal  surface 
of  the  nozzle  by  approximately  0.050  inch.  The  amount  of  the  profusion 
was  accurately  measured  and  then  machined  flush  with  the  nozzle  contour. 


The  plug  was  bonded  to  the  ablative  part  using  a  low  viscos¬ 
ity  epoxy  resin  (Epon  815  and  DTA) .  The  bond  line  included  the 
external  surface  of  the  plug  with  the  exception  of  the  surface 
exposed  to  the  internal  surface  of  the  nozzle.  In  addition,  the 
cavity  at  the  backside  of  _the_plug  is  partially  filled  with  epoxy 
resin.  The  thickness  of  this  resin  is  approximately  equal  to  the 
length  of  the  strain  relieving  pins,  and,  thus,  it  protects  the 
wire  connections  made  at  the  back  of  the  thermocoup led  plug.  Dur¬ 
ing  the  installation  of  the  plug,  the  thermocouple  lead  wires  are 
bent  and  led  out  in  the  proper  direction  for  installation  of  the 
ablative  nozzle  into  the  steel  housing. 


4. 1.2. 4  Recording  Instrumentation 

Standard  recording  instrumentation  was  used  for  the  measure¬ 
ment  system.  In  the  Phase  II  evaluation  test  program,  the  thermo¬ 
couple  outputs  were  recorded  on  an  oscillograph  continuously  through 
each  test.  The  recording  circuit  used  for  all  thermocouples' is  pre¬ 
sented  schematically  in  Figure  4-24.  A  room  temperature  reference 
junction  was  used;  the  reference  temperature  was  irieasured  with  a 
mercury-in-glass  thermometer.  The  circuit  sensitivity,  millivolts  per 
inch  of  galvanometer  deflection,  was  adjusted  to  yield  a  maximum  of 
6  inches  of  deflection  (the  usable  range  of  the  galvanometer)  for  the 
maximum  temperature  anticipated  for  each  thermocouple.  The  galva¬ 
nometers  used  and  the  circuit  resistance  were  matched  to  yield  the 
proper  galvanometer  damping  characteristics. 


Each  thermocouple  circuit  was  calibrated  just  prior  to  test  using 
an  external  emf  source  (Figure  4-24).  The  calibration  consists  of 
imposing  four  or  five  emf  values  in  steps  up  to  the  maximum  anticipated 
in  the  test  and  simultaneously  recording  the  galvanometer  deflections. 
The  precise  calibration  emf  values  are  determined  with  a  potentiometer. 


In  the  Phase  III  char  motor  firings,  the  thermocouple  outputs 
were  measured  with  the  RPL  SEL-600  digital  voltmeter  system  and 
recorded  on  magnetic  tape.  A  150°F  oven  provided  the  reference  temp-  ' 
erature  for  all  thermocouples.  No  circuit  calibration  was  required 
since  direct  voltage  measurements  were  made. 
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FIGURE  4-24  PHASK  II  THERMOCOUPLE  tECOUIK 
«K>  CALIBHATIHG  CIRCUIT 


4.2  RADIATIVE  HEAT  FLUX  SYSTEM 

The  radiative  heat  flux  measurement  system  is  dependent  on 
measuring  the  temperature  difference  from  the  center  to  the  edge  of 
the  constantan  foil  in  the  narrow  view  angle  radiometer  (Section  3.2, 
Figure  3-4) .  This  temperature  difference  is  recorded  as  an  emf  out¬ 
put  which  is  converted  to  incident  heat’  flux  based  on  the  instrument 
calibration.  In  making  this  conversion,  potential  sources  of  error 
include  the  wavelength  characteristics  of  the  lens  and  sensor  surface 
coating,  the  response  time  of  the  constantan  foil  sensor,  the  procedure 
used  for  instrument  calibration, and  the  copper  heat  sink  temperature. 
These  potential  errors  are  discussed  in  Section  4.2.1.  Once  the 
potential  sources  of  error  are  defined,  the  instrument  design,  fabri¬ 
cation,  and  installation  are  defined  in  Section  4.2.2. 


As  mentioned  in  Section  1  and  described  in  more  detail  in 
Section  7,  two  40-inch  char  motor  tests  were  conducted  under  this 
contract.  Each  firing  had  a  narrow  view  angle  radiometer  at  the  nozzle 
throat  plane.  In  addition,  the  second  char  motor  firing  had  a  narrow 
view  angle  radiometer  of  a  different  configuration  at  a  supersonic  area 
ratio  of  3.0.  The  design  of  this  radiometer  was  dictated  by  the 
requirements  of  a  high  chamber  pressure  motor  (Reference  4-18) ,  and 
these  requirements  will  not  be  discussed  here.  Therefore,  the  error 
analysis  presented  in  Section  4.2.1  is  directed  to  the  radiometer 
located  at  the  nozzle  throat  plane.  However  both  radiometers  are 
discussed  in  Section  4.2.2. 
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The  errors  associated  with  measuring  the  narrow  view  angle 
radiation  flux  using  the  radiometer  shown  in  Figure  3-4  can  be 
categorized  as  those  resulting  from  thermal  effects,  optical  effects, 
response  time,  and  method  of  instrument  calibration.  The  thermal 
effects  include  the  variation  of  emf  output  with  the  copper  sink 
temperature  and  the  variation  of  the  elemental  view  factor  across  the 
constantan  disk,  and  these  effects  are  defined  in  Section  4. 2. 1.1. 

The  optical  effects  result  from  the  variation  of  the  lens  transmissivity 
and  of  the  sensor  surface  coatings  absorbtivity  with  wavelength  and  from 
the  coupling  of  this  variation  with  the  wavelength  characteristics  of 
the  radiative  source.  The  optical  effects  are  defined  in  Section 
4. 2. 1.2.  The  error  associated  with  the  time  constant  of  the  instru¬ 
ment  is  dependent  on  the  fluctuations  which  are  expected  in  the  emissive 
power  of  the  radiative  source,  and  this  is  discussed  in  Section  4. 2. 1.3. 
The  final  error  associated  with  measuring  the  radiative  heat  flux  is 
the  method  of  instrument  calibration,  and  this  is  discussed  in  Section 
4. 2. 1.4.  Finally,  the  errors  from  the  various  sources  are  combined 
to  give  the  overall  system  accuracy  in  Section  4. 2. 1.5. 


4. 2. 1.1  Thermal  Effects 

The  recorded  output  of  the  narrow  view  angle  radiometer  is 
an  emf  which  is  dependent  on  the  temperature  difference  across  the 
constantan  disk  which  is  in  turn  dependent  on  the  incident  radiative 
heat  flux.  The  temperature  distribution  for  the  constantan  foil 
sensor  shown  in  the  sketch  below  is  given  by 
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From  Equation  (4-10) »  one  notes  that  for  a  constant  heat  flux  the 
temperature  distribution  is  independent  of  the  copper  sink  temperature 
provided  the  thermal  conductivity  of  constantan  is  independent  of 
temperature.  Since  the  thermal  conductivity  of  constantan  increases 
with  increasing  temperature,  the  temperature  difference  across  the 
sensor  decreases  with  increasing  copper  sink  temperature.  If  the  •• 
thermal  conductivity  of  constantan  is  represented  by 

K  =  K  (1+aT) 
o 

T  -  is  the  temperature  above  the  reference  temperature . 
(In  this  analysis  T=Tr-;R  is  taken  as  the  reference 
temperature) . 

a  -  1.28  x  10_3./°F 

Kq  -  12.6  Btu/hr  ft°F  @  32 °F 

q  -  10  Btu/ft2sec 
s  -  0.0005  in 
R  -  0.0775  in 


where 


and  if 


then,  the  temperature  distributions  for  copper  sink  temperatures  of 
100,  200,  and  300°F  are  shown  in  Figure  4-25.  The  temperature  difference 
for  sink  temperatures  of  100  and  300°F  is  shown  in  Figure  4-26a  as 
a  function  of  heat  flux.  Due  to  the  variation  of  the  temperature  differ 
ence  (Tr=o  ~  Tr=R)  with  sink  temperature,  the  recorded  emf  output  also 
varies.  However,  one  notes  from  the  standard  emf  charts  for  copper/ 
constantan  that  the  measured  emf  for  a  constant  temperature  difference 
is  dependent  on  the  temperature  level.  This  dependency  is  shown  in  Fig¬ 
ure  4-26b.  When  the  results  of  Figures  4- 26a  and  4-26b  are  combined,  the 
effect  of  sink  temperature  on  the  relationship  of  incident  heat  flux 
vs.  recorded  emf  is  shown  in  Figure  4-27.  For  a  sink  temperature  vari¬ 
ation  of  200Pf,  the  error  in  heat  flux  is  -3  percent. 


Another  source  of  error  which  is  caused  by  an  increase  in  the  cop¬ 
per  sink  temperature  is  due  to  reradiation  from  the  sensor.  As  shown 
previously  in  Figure  4-24, the  maximum  sensor  temperature  increases 
with  increasing  sink  temperature  for  a  constant  incident  heat  flux. 

This  results  in  an  increase  in  the  effective  sensor  temperature  which 
is  defined  by 


T  =  T  +  0.75AT 
e  o 
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FIGURE  4~26a  VARIATION  OF  TEMPERATURE 

DIFFERENCE  WITH  TEMPERATURE  LEVEL 
AND  INCIDENT  HEAT  FLUX 


FIGURE  4- 2 6b  VARIATION  OF  RECORDED  EMF  OUTPUT  WITH 
TEMPERATURE  DIFFERENCE  AND  TEMPERATURE 
LEVEL 


For  the  maximum  sensor  heat  flux,  the  reradiative  flux,  defined  using 
the  effective  sensor  temperature,  is  shown  in  Figure  4-28  as  a  function 
of  the  heat  sink  temperature.  The  results  shown  in  this  figure" are 
just  from  the  blackened  face  of  the  sensor.  Due  to  the  small  emissivity 
(e=0.06)  of  constantan,  reradiation  from  the  opposite  face  of  the 
sensor  is  negligible.  For  the  situation  of  the  maximum  heat  flux, 
an  increase  in  sink  temperature  from  100  to  300°F  causes  a  -3  percent 
error  in  the  measured  heat  flux. 


4. 2. 1.2  Optical  Effects 

The  narrow  view  angle  radiometer  has  a  lens  which  transmits 
the  radiative  heat  flux  to  the  sensor  while  protecting  it  from  any 
convective  heating.  In  addition,  a  special  coating  is  put  on  the  N 
sensor  so  that  it  acts  as  much  like  a  black  body  as  possible. 

However,  both  the  transmissivity  of  the  lens  and  the  absorbtivity 
of  the  sensor  coating  exhibit  some  wavelength  dependency,  and  this 
dependency  is  shown  in  Figure  4-29  for  the  sapphire  lens  and  coating 
materials  used  by  Heat  Technology  Laboratory  in  their  fabrication  of 
narrow  view  angle  radiometers.  One  notes  from  this  figure  that  the 
absorbtivity  of  the  sensor  coating  is  approximately  0.90  and  increases 


nans  i-Ji  afproxjmats  radiation  uat  loss  from  sensor 
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slightly  with  wavelength  while  the  transmissivity  of  the  sapphire 
lens  is  approximately  constant  between  0,5  and  4.5  microns  and 
decreases  for  wavelengths  outside  of  this  range.  The  wavelength 
dependency  of  the  coating  and  of  the  sapphire  lens  can  introduce 
an  error  into  the  measurement  of  the  radiative  heati  flux  if  the 
emissive  power  wavelength  .characteristics  of  the  source  being 
measured  are  different  than  those  of  the  calibration  source.  The 
significance  of  this  error  must  be  evaluated  for  the  particular 
application  of  interest  and  is  evaluated  below  for  temperatures 
representative  of  those  encountered  in/ solid  rocket  motors. 

Since  the  calibration  source  simulates  a  blackbody  source,  its 
emissive  power  wavelength  characteristics  as  a  function  of  source 
temperature  are  adequately  described  by  the  blackbody  relation¬ 
ships,  and  these  are  shown  in  Figure  4-30  -for  source  temperatures 
of  3000,  5000,  and  7000°R.  One  notes  that  the  emissive  power  is 
concentrated  in  the  shorter  wavelengths  as  the  source  temperature 
increases.  The  effect  on  the  measured  heat  flux  of  the  wavelength 
characteristics  of  the  source  combined  with  the  wavelength  charac¬ 
teristics  of  the  sapphire  lens  and  sensor  coating  is  best  illus¬ 
trated  by  an  example.  In  ^he  first  case  of  this  example,  the  heat 
flux  incident  on  the  sapphire  lens  is  from  a/ blackbody  source  at 
400u°R,  and,  in  the  second  case,  jfche  same  heat  flux  is  generated 
by  a  6000°R  gray  body  source  with  an  emissivfity  of  0.135.  For 
these  cases,  the  heat  flux  absorbed  by  the  censor  is  3%  less  for 
the  blackbody  source  than  fqr  the  gray  body  source.  If  a  similar 
example  was  performed  which  considered  a  black  body  source  at 
3000°R  and  a  gray  body  source  at  6600°R  that  emitted  the  identical 
heat  flux,  the  difference  in  the  heat  flux  absorbed  by  the  .sensor 
would  be  greater  than  3%.  This  is  due  to  the  higher  attentuation 
of  the  black  body  source  by  the  sapphire  lens. 


4. 2. 1.3 


Response  Time 


The  response  time  o£_a^Gardon  type  radiometer 
erence  3-2  as 

_2 
cpR 

T  =  TIT 

* 

For  the  application  being  considered  here,  the  radius  of  the  sensor 
is  approximately  0.0775  inches.  Thus,  for  a  constantan  foil,  this 
gives  a  time  constant  of  approximately  0.10  seconds.  Whether  a  time 
constant  of  this  magnitude  is  acceptable  depends  on  the  characteristics 
of  the  heat  spurce  which  are  to  be  accurately  recorded.  As  described 
earlier,  the  source  viewed  by  the  narrow  view  angle  radiometer  is  a 
combination  of  gas  radiation,  particle  radiation,  and  opposite  wall 
radiation.  Since  start  up  transient  information  concerning  the 
radiancy  of  the  opposite  wall  and  the  radiancy  of  the  particle  laden 
combustion  products  would  be  valuable,  it  is  desirable  to  have  a 
time  constant  which  is  ^n  order  of  magnitude  less  than  the  time  period 
in  which  these  initial  start  up  transients  occur.  In  the-following 
paragraphs,  an  analysis  is  performed  which- shows  the  response  of  a 
Gardon  type  narrow  view  angle  radiometer  to  a  sinusoidally  time 
varying  heat  flux  input.  This  shape  was  chosen  for  the  heat  flux 
input  because  a  quarter  cycle  is  a  reasonable  approximation  for  the 
radiancy  of  the  opposite  wall  and  the  alumiha  particles  during  the 
motor  start  up.' 


As  described  in  Reference  3-2  the  rate  of  temperature  rise  of  the 
circular  foil  sensor  can  be  represented  by  an  exponential  law  such  as 

T  =  T*  (l-e~t/,T)  (4-12) 

>  Where  T  -  steady  state  temperature 

Since  the  temperature  of  the  sensor  is  approximately  linear  with  heat 
flux,  the  above  relation  also  implies  that 


is  given  in  Ref 

(4-11) 


[’  =  q  (l-e"t/T 


(4-13) 


where  q  -  actual  calorimeter  absorbed  flux 
q’  -  heat  flux  inferred  (recorded) 

The  sinusoidal  profile  selected  for  this  analysis  was 

q  =  q  sin  ~ 


(4-14) 


where 


qo  -  peak  heat  flux 

4  -  twice  the  motor  ignition 

•  transient  time 
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Now  the  differential  equation  whose  solution  is  Equation  (4-13)  is 


dt  4 


(4-15) 


By  substituting  Equation  (4-14)  into  Equation  (4-15)  and  solving  the 
differential  equation  with  the  appropriate  boundary  condition,  one 
obtains  the  results  shown  in  Figure  4-31  Due  to  the  large  free  volume 
of  che  char  motor  the  ignition  transient  is  on  the  order  of  five 
seconds  which,  for  a  time  constant  of  0.1  seconds,  corresponds  to  a 
t A  ratio  of  0.01.  For  this  ratio,  one  notes  from  Figure  4-31  that 
the  recorded  heat  flux  is  in  close  agreement  with  the  incident  heat 
flux. 

4. 2. 1.4  Calibration 

The  narrow  view  angle  radiometer  is  calibrated  by  exposing  it 
to  a  radiation  source  which  is  an  electrically  heated  graphite  block. 
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FIGURE  4-31  NON-DIMENSIONAL  CALORIMETER  RESPONSE  AS  A 
FUNCTION  OF  NON-DIMENSIONAL  TIME 


The  uniformity  of  the  heat  flux  from  the  graphite  block  has  been 
carefully  mapped  using  a  standard  view  angle  radiometer.  During 
the  calibration  procedure,  th6  transducer  to  be  calibrated  is  held 
in  a  water  cooled  holder,  and  the  aperture  of  the  transducer  is 
facing  the  graphite  block.  The  distance  between  the  aperture  and 
the  graphite  block  is  such  that  the  field  of  view  of  the  transducer 
is  located  within  the  uniform  heat  flux  region  of  the  graphite 
block.  The  standard  transducer  is  located  on  the  opposite  side  of 
the  uniform  temperature  block,  but  it  views  the  same  region.  The 
transducer  to  be  calibrated  has  the  lens  in  place  and  the  sensor 
coating  applied  prior  to  calibration.  The  upper  limit  on  the 
calibration  heat  flux  is  approximately  250  Btu/ft^sec.  This  means 
that  the  calibration  curve  for  transducers  whose  range  is  outside 
this  limit  must  be  obtained  by  extrapolating  the  available  calibra¬ 
tion  data.  For  Gardon  type  radiometers,  a  straight  line  extrapola¬ 
tion  of  the  available  calibration  data  is  used  for  the  higher 
(725  Btu/ft2sec)  heat  flux  range.  This  extrapolation  is  in  agree¬ 
ment  with  that  which  is  predicted  by  analytical  methods.  The  out¬ 
put  from  the  calibration  procedure  is  a  curve  which  relates  incident  a 
hemispherical  heat  flux  to  recorded  emf  output.  The  incident  heat 
flux  is  defined  using  the  output  from  the  standard  transducer. 


The  important  items  to  note  from  the  above  description  of  the 
calibration  procedure  are  1)  the  transducer  is  in  a  water  cooled 
holder  and,  thus,  its  sink  temperature  is  approximately  ambient 
temperature,  2)  the  lens  is  in  place  during  calibration,  3)  the  limit 
on  the  calibration  heat  flux  is  250  Btu/ft2sec  and  4)  the  calibration 
curve  relates  total  hemispherical  heat  flux  to  recorded  emf  output. 
The  uncertainty  in  the  calibration  procedure  is  stated  by  the  manu¬ 
facturer  to  be  +2%  which  includes  nonlinearity  effects.  This 
uncertainty  is  In  addition  to  the  errors  described  in  the  previous 
sections . 


4 . 2 . 1 . 5  Overall  Outpujt  Accuracy 

r  V 

In  Sections  4. 2. 1.1  through  4. 2. 1.4,  the  errors  associated 
with  the  use  of  a  narrow  view  angle  radiometer  to  measure  the  total 
radiancy  at  a  nozzle  throat  location  were  discussed.  These  errors 
which  are  summarized  in  Columns  1  and  2  of  Table  4-5  are  uncorrected 
errors  as  no  effort  was  made  to  reduce  their  magnitude  through  the 
use  of  either  analytical  techniques  or  experimental  data.  The  errors 
shown  in  Column  1  are  the  maximum  errors  and  would  occur  only  if  the 
transducer  sink  temperature  became  300°F  and  if  the  wavelength 
characteristics  of  the  source  simulate  a  black  body.  Since  the 
expected  maximum  sink  temperature  for  the  radiometer  being  used  in 
rocket  nozzle  is  200°Fr  the  probable  error  dtfe  to  this  effect  is 
1t1/2%.  The  probable  errors  resulting  from  radiation  heat  loss  from 
the  sensor  and  from  differences  in  the  temperature  and  wavelength 
characteristics  between  the  calibration  and  measurement  sources 
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TABLE  4-5 


RADIOMETER  ACCURACY 


["sources  of  Error 

Oncorrected 
Maximum  Error 
(Percent) 

Uncorrected 
Probable  Error 
(Percent) 

Corrected 

Error 

(Percent) 

Calibration  Shift  with 
Body  Temperature 

-3 

-1.5 

-0.5 

Radiation  Loss  from 
Sensor 

-7 

-3.5 

+0.5/-1 

Temperature  Differ¬ 
ence  between 
Calibration/ 
Measurement  Sources 

+3 

+1 

♦0.5 

^Calibration  Uncertainty 
(Nonlinearity  and 

+2 

+2 

♦2 

Error) 

Overall  Error 

-S/-9 

-2/-6 

+2. 5/- 3 

are  also  less  than  the  uncorrected  maximum  errors  as  shown  by  a 
comparison  of  Columns  1  and  2. 


Since  the  first  two  errors  shown  in  Table  4-5  result  from  a 
temperature  rise  of  the  transducer  sink  temperature,  a  knowledge  of 
this  temperature  would  permit  the  user  to  correct  the  recorded  emf 
output  for  these  errors,  and,  thus,  they  could  be  essentially  elimin¬ 
ated.  THe  requited  information  is  obtained  from  a  thermocouple  which 
is  attached  to  the  copper  sink  and  the  reduced  error  is  shown  by 
Column  3  of  the  table.  Even  with  this  sink  temperature,  the  errors 
cannot  be  entirely  eliminated  because  of  such  items  as  the  possibility 
of  temperature  gradients  within  the  copper  heat  sink,  the  possibility 
of  the  aperature  and  the  copper  barrel  increasing  in  temperature  and, 
thus,  influencing  the  heat  flux  to  the  sensor,  and  the  possibility 
of  higher  heat  conduction  loss  in  the  thermocouple  wire  with  increasing 
sensor  temperature. 


Although  the  radiancy  from  the  alumina  particles  in  a  solid 
rocket  plume  as  a  function  of  wavelength  and  temperature  is  not  well 
understood,  the  error  resulting  from  the  temperature  difference  between 
the  calibration  and  measurement  sources  can  be  reduced  assuming  the 
calibration  source  is  a  black  body  and  the  measurement  source  is  a 
gray  body.  Due  to  the  uncertainties  in  the  characteristics  of  the 
measurement  source,  this  reduction  in  error  can  be  considered  only 


a  first  order  correction.  More  knowledge  of  the  measurement  source 
must  be  available  before  any  further  corrections  could  be  made. 


The  overall  radiometer  accuracy  for  the  three  cases  studied 
(uncorrected  maximum  error,  uncorrected  probable  error,  and  corrected 
error)  is  determined  from  the  following  relationship 


error  =  £ (fixed  errors)  +  \E (random  errors) 


All  the  errors  are  fixed  errors  except  the  calibration  uncertainty 
which  is  a  random  error.  The  overall  accuracy  for  the  three  cases 
are  tabulated  on  the  last  line  of  Table  4-5 


4.2.2  Design  and  Fabrication  Analysis 

Narrow  view  angle  radiometers  have  been  used  extensively  to 
measure  the  total  radiative  heat  flux,  and  they  are  commercially 
available  as  off-the-shelf  or  pustcm-built  instruments.  Design  or 
configuration  options  include  size,  shape,  view  angle,  gas  purge, 
optical  lens,  and  water  cooling.  The  design  considerations  which  were 
evaluated  in  defining  the  transducer  requirements  are  discussed  in 
Section  4. 2. 2.1.  The  incorporation  of  the  transducer  requirements, 
design  considerations,  and  fabrication  techniques  into  an  overall 
system  configuration  is  discussed  in  Section  4. 2. 2. 2.  The  installation 
of  the  transducer  into  the  rocket  nozzle  is  discussed  in  Section 
4J2.2.3,  arid,  finally,  the  recording  instrumentation  and  associated 
hardware  required  by  the  transducer  are  defined  in  Section  4. 2. 2. 4. 


4.2.2. 1  Design  Considerations 

Once  the  requirements  of  a  narrow  view  angle  radiometer 
have  been  defined  the  radiometer  design  is  formulated  based  on  satis¬ 
fying  certain  geometrical  constraints  and  cooling  requirements.-  -Most- 
of  the  transducer  requirements  have  been  defined  in  Section  3.2  but 
they  are  briefly  summarized  below: 

•  Continuous  reading  for  the  duration  of  the 
rocket  motor  firing 

•  Must  not  jeopardize  the  integrity  of  the 
nozzle  or  cause  a  significant  degradation 
in  its  performance 

•  Must  be  within  the  current  state-of-the- 
art 
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In  addition  to  these  general  transducer  requirements,  specific 
requirements  are  that  the  recorded  millivolt  output  should  be  approxi¬ 
mately  10  mv  for  the  maximum  expected  radiative  flux,  and  that  the 
maximum  temperature  of  the  constantan  sensor  should  not  exceed 
approximately  700°F.  This  latter  requirement  is  due  primarily  to 
the  emf  falling  off  rapidly  above  this  temperature.  Based  on  these 
requirements,  the  solution  of  the  geometrical  and  cooling  problems 
are  discussed  separately  in  Section  4. 2. 2. 1.1  and  4. 2. 2. 1.2, 
respectively. 


4. 2. 2. 1.1  Geometrical  Constraints .  Before  a  preliminary  design 
of  the  radiometer  can  be  developed,  the  geometrical  constraints,* 
which  include  the  cavity  diameter,  depth  of  the  transducer  aperture 
below  the  original  nozzle  surface,  and  the  total  thickness  of  the 
ablative  component  at  the  measurement  location  must  be  defined. 

With  the  exception  of  the  latter  item  which  is  fixed  by  the  nozzle 
design,  there  is  some  latitude  in  defining  the  remaining  two  items. 
In  the  paragraphs  below,  the  general  guidelines  in  accomplishing 
the  radiometer  design  are  described,  and  the  geometry  of  the 
final  design  is  presented.  This  discussion  is  based  on  the  measure¬ 
ment  system  at  the  throat  of  the  char  motor  nozzle  (Section  7) . 


From  Figure  4-1,  one  notes  that  in  order  to  prevent  the  aperture 
of  the  radiometer  from  reaching  the  nozzle  surface  it  must  be 
recessed  greater  than  0.386  inches.  This  figure  also  shows  that 
to  prevent  excessive  temperatures  of  the  copper  shell  due  to  heat 
conduction  from  the  carbon  phenolic  ablative  component,  the  radi¬ 
ometer  should  be  located  approximately  0.6  inches  below  the  nozzle 
surface  at  motor  burnout.  In  addition,  one  notes  from  the  descrip¬ 
tion  of  the  char  motor  nozzle  design  presented  in  Section  7.1  that 
the  ablative  component  has  a  thickness  of  2.7  inches.  With  these 
constraints,  the  radiometer  would  be  approximately  1.7  inches  long 
and  its  aperture  would  be  approximately  one  inch  below  the  original 
nozzle  surface. 


Design  constraints  related  by  the  response  of  the  narrow  view 
angle  radiometer  are  dictated  by  the  following  relations: 


where 


.  .  2 
q  =  q  sin  (a/2) 
o 

a  -  view  angle 
q  -  hemispherical  radiation 
flux  emitted  by  source 


(4-16) 
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T  -  T2  -  A 


(4-17) 


E/-  =  16.7  - 
q  s 


(4-18) 


Tr=0  =  £6* 


(4-19) 


Equation  (4-16)  defines  the  heat  flux  incident. to  the  radiometer  as 
a  function  of  the  view  angle.  The  ratio  of  q/q  is  shown  in  Figure 
4-32  as  a  function  of  the  view  angle  for  the  range  of  interest  here. 
Also  shown  in  this  figure  is  the  variation  across  the  sensor  of  the 
incident  heat  flux.  Equation  (4-17)  defines  the  temperature  rise  » 
across  a  thin  circular  disk  sensor  for  a  specific  incident  heat 
flux;  Equation  (4-18)  defines  the  transducer  sensitivity. 
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FIGURE  4-32  EFFECT  OF  VIEW  ANCLE  ON  MEASURED  MEAT  FLOE 


and  Equation  (4-19)  defines  the  transducer  time  constant.  In 
addition  to  the  above  relationships,  two  other  factors  which  must  be 
considered  are  1)  the  relationship  between  the  view  angle  and  the 
effective  view  angle  as  shown  schematically  by  the  sketch  below  and 
2)  the  requirement  that,  if  only  10%  or  less  of  the  total  (hemispheri¬ 
cal)  emitted  flux  is  incident  on  the  sensor,  the  aperture  and  sensor 
diameters  must  be  approximately  equal.  The  first  factor  is  important 


in  defining  the  cavity  depth  below  the  ablative  surface  before  the 
sensor  sees  the  cavity  walls,  and  the  second  is  necessary  if  the 
elemental  view  factor  at  any  location  on  the  sensor  is  to  be 
essentially  constant,  thus  minimizing  any  nonlinear  effect  which  would 
have  been  introduced  into  the  transducer  design. 

Before  considering  the  transducer  design  using  the  above  infor¬ 
mation  one  should  note  from  Equation  (4-18)  that  the  transducer  sensi¬ 
tivity  is  inversely  proportional  to  the  sensor  thickness.  Since  it 
is  desirable  to  maximize  the  instrument  sensitivity,  the  sensor  thick¬ 
ness  should  be  minimized.  The  practical  minimum  thickness  is 
0.0005  inches.  This  bit  of  information  together  with  the  desire  to 
have  an  emf  output  of  10  mv  allows  one  to  solve  Equations  (4-16) 

(  (4-17)  ,  and  (4-18)  for  q,  T,  and  R  as  a  function  of  view  angle.  By 

iterating  using  the  view  angle  as  the  independent  variable,  a  narrow 
view  angle  radiometer  was  defined  which  satisfactorily  met  all  the 
requirements.  The  pertinent  characteristics  of  this  radiometer  are 
tabulated  in  Table  4-6.  With  the  effective  view  angle  shown  in  this 
table,  the  aperture  in  a  0.5  inch  diameter  cavity  can  be  recessed  no 
more  than  0.725  inches.  This  is  based  on  a  sensor  diameter  of  0.156 
inches.  One  also  notes  from  the  transducer  sensitivity  that  a  10  mv 
output  is  obtained  for  a  source  radiancy  of  630  Btu/ft^sec  which 
approximates  the  maximum  expected  radiancy  of  570  Btu/ft^sec.  In 
addition,  as  described  in  Section  4. 2. 1.3  the  time  constant  is  suffi¬ 
ciently  small  that  the  transducer  can  follow  the  ignition  trans¬ 
ients. 


As  described  previously,  the  requirements  for  the  radiometer 
located  at  an  expansion  ratio  of  3.0  were  dictated  by  the  requirements 
of  a  high  chamber  pressure  motor.  The  characteristics  of  this  radiom¬ 
eter  are  tabulated  in  Table  4-7. 


TABLE  4-6 


CHARACTERISTICS  OF  NARROW  VIEW  ANGLE 
RADIOMETER  LOCATED  AT  THE  NOZZLE  THROAT 


View  Angle,  Degrees 

12.1 

Effective  View  Angle,  Degrees 

26.4 

Sensitivity,  mv/ (Btu/f t2sec) 

9.0 (at  570  Btu/ft2sec) 

Time  Constant,  sec- 

0.194 

TABLE  4-7 

CHARACTERISTICS  OF  THE  NARROW  VIEW  ANGLE  RADIOMETER 
LOCATED  AT  AN  EXPANSION  RATIO  OF  3.0 


View  Angle,  Degrees 

12.2 

Effective  View  Angle,  Degrees 

24.0 

Sensitivity,  mv/ (Btu/ft2sec) 

9.0 (at  570  Btu/ft2sec) 

Time  Constant,  sec. 

0.190 

\  \ 

4. 2. 2. 1.2  Cooling  Requirements.  Once  the  basic  ^dimensions  of  the 
transducer  Rave  been  defined ,  the  next  problem  is  to  define  its  cool¬ 
ing  requirements.  Heat  is  transferred  to  the  transducer  by  radiation 
heating  of  the  sensor,  aperture,  and  inner  cylindrical  walls,  by  con¬ 
vective  heating  of  the  aperture  from  the  hot  combustion  chamber  gases 
in  the  cavity,  and  by  conductive  heating  of  the  copper  housing  from 
the  ablative,  component.  The  introduction  of  a  gas  purge  Vhich  flows 
around  the  sapphire  lens  and  exits  through  the  aperture  provides 
an  effective  method  of  eliminating  the  convective  heat  transfer  from 
the  combustion  products.  Another  advantage  of  this  purge  gas  is 
that  it  prevents  any  alumina  particles  or  ablative'1  surface  material 
residue  from  contaminating  the  lens.  This  purge  gas  flow  rate  was 
approximately  defined  by  making  an  energy  balance^  on  the  control 
volume  shown  in  Figure  4-33.  The  equation  which  describes  this  energy 
balance  is 


+ 


peUeCH^aw  hc^Ac 


(h.  -  h..  ) 

le  ix 


«c  *  •  *£> 


FIGURE  4-33  CONTROL  VOLUME  FOR  DEFINING 
GAS  FORGE  FLOW  RATE 
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.Solving  the  above  equation  for  purge  flow  rate  as  a  function  of  cavity 
temperature  gives  the  results  shown  in  Figure  4-34  for  both  the  char 
motor  firing  and  the  arc  plasma  generator  tests.  These  results  show 
that  an  average  cavity  temperature  of  1000°R  can  be  maintained  for 
purge  gas  flow  rates  oi^the  order  of  0.010  Ibs/sec.  In  order  to  be 
assured  that  this  «flow  rate  is  obtained,  the  stagnation  pressure  for 
the  purg*e  system ‘must  be  approximately  twice  the  nozzle  throat  static 
pressure.  For  the  motor  being  considered  here,  this  means  that  the 
purge  system  stagnation  pressure  must  be  on  the  order  of  1500  psia. 


The  conductive  heat  transfer  from  the  ablative  component  to  the 
transducer  is  a  result  of  heat  penetration  from  the  nozzle  surface. 
Previously,  it  was  mentioned  that  the  aperture  should  be  approximately 
1 . 0  .inches  below  the  original  nozzle  surface  in  order  to  keep  the  tip 
of  the  transducer  below  1000?R.  However,  geometrical  limitations  make 
i,t  impossible  to  have  the  transducer  recessed  this  deep.  To  solve  the 
problems  of  both  the  conductive  heating  and  the  geometrical  limitation, 
the  barrel  of  the  transduce!:  was  tapered  as  shown  in  Figure  4-35. 


With  the  convective  and  conductive  heat  transfer  to  the  transducer 
minimized,  acceptable  copper  temperatures  can  be  maintained  by  using  a 
heat  sink  cooling  approach  without  having  to  use  water  cooling.  The 
copper  flange  shown  in  Figure  4-17  was  sized  to  prevent  the  average 
copper  temperature  at  the  end  of  a  60  second  motor  firing  from  exceeding 
200°F.  Advantages  of  the  heat  sink,  rather  than  the  water  cooling 
approach,  are  that  it  simplifies  the  radiometer  design,  simplifies 
the  nozzle  assembly  and  on  site  supporting  hardware,  and  simplifies  the 
reduction  of  data.' 


The  external  configuration  for  the  high  pressure  radiometer  is 
shown  in  Figure  4-36.  The  components  for  this  configuration  are 
basically  the  same  as  those  described  for  the  radiometer  shown  in 
Figure  3-4.  The  configuration  was  dictated  primarily  by  the  installa? 
tion  and  high  purge  pressure  requirements  for  the  high  chamber  pressure 
motors  to  be  fired  under  Contract  F04611-69-C-0065 .*  , 


4.2. 2.2  System  Configuration 

Based  on  the  reasoning  described  in  Section  3.2,  a  narrow  view 
angle  radiometer  was  selected  to  measure  the  radiative  heat-  flux  from 
the  opposite  wall  and  the  combustion  products.0  These  transducers  are 
available  as  off  the  shelf  items  in  several  basic  configurations.  How¬ 
ever,  the  configurations  described  in  Section  4. 2. 2.1  are  not  standard 
configurations,  and,  thus,  these  transducers  had  to  be  fabricated  as  a 
specialty  item.  After  Aerotherm  had  developed  the  requirements  for  the 

*  11 

High  Chamber  Pressure  Nozzle  Material  Development,  United  Technology 
Center,  Sunnyvale,  California 

4-58 


/ 


l 


NITROGEN  PURGE  GAS  FLOW  RATE  LI/SE 


CHAR  WOT 01  Fill MG 


KJUUtOK  ANGLE  RADIOMETER 


■  7S0.ltUIIF.3A 


riGUl'E  4-3*  RICH  PRESSURE  RADIOMETER 


radiometers.  Heat  Technology  Laboratory  performed  the  final  design  and 
fabricated  the  transducer.  The  paragraphs  below  describe  the  system 
configuration. 


A  schematic  and  the  definition  of  the  external  configuration  of 
the  narrow  view  angle  radiometer  located  at  the  nozzle  throat  plane 
have  been  shown  in  Figures  3-4  and  4-35,  respectively;  the  external 
configuration  for  the  high  pressure  radiometer  was  shown  in  Figure 
4-36.  Since  these  radiometers  are  identical  in  basic  details,  the 
general  description  presented  below  is  applicable  to  both. 


The  size  of  a  narrow  view  angle  radiometer  is  dependent  on 
numerous  variables  which  include  environmental,  configurational,  and 
design  constraints.  For  measuring  the  radiant  heat  flux  in  solid 
rocket  motors,  a  practical  aperture  diameter  is  on  the  order  of  0.20 
inches  and  a  practical  diameter  for  the  cylindrical  barrel  is  on  the 
order  of  0.50  to  0.75  inches.  The  copper  heat  sink  is  sized  to 
satisfy  cooling  requirements  and  a  disk  of  copper  which  is  1.50  inches 
in  diameter  and  0.50  inches  thick  is  sufficient  to  satisfy  those 
requirements  which  were  considered  here.  The  sensormust  be  protected 
with  a  lens  to  eliminate  convective  heating  to  the  sensor  surface  and 
a  gas  purge  is  required  primarily  to  prevent  contamination  of  the  lens 
due  to  alumina  particles  or  ablation  products.  A  secondary  benefit  of 
th;^  gas  purge  is  to  cool  the  radiometer.  The  recorded  output  from 
the  radiometer  is  the  emf  output  corresponding  to  the  temperature 
difference  between  the  center  and  the  edge  of  the  sensor,  and  this  emf 
output  is  related  to  radiant  heat  flux  by  calibrating  the  instrument 
using  a  known  heat  flux  source. 


/ 
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The  housing  of  the  radiometer  as  well  as  the  heat  sink  is  made 
from  copper  primarily  because  of  the  high  thermal  conductivity 
of  this  material  which  insures  that  the  sink  is  at  a  uniform  tempera¬ 
ture.  The  sensor  is  composed  of  a  thin  blackened  constantan  disk. 

This  material  is  particularly  attractive  because  its  thermal  emf 
characteristics  are  known  and  because  its  sensitivity  is  approximately 
linear  over  the  temperature  range  of  interest.  The  sensor  surface  is 
blackened  so  that  it  absorbs  all  but  a  small  percentage  of  the  inci¬ 
dent  radiative  flux.  The  gas  purge  line  is  a  thin  walled  stainless 
steel  tube  which  is  brazed  to  the  copper  housing.  This  tube  must  be 
flexible,  and  it  must  be  adequately  joined  with  the  copper  housing 
so  that  a  moderate  amount  of  bending  can  be  tolerated  during  the 
assembly  of  the  radiometer  in  the  steel  housing. 


The  most  critical  item  in  assembling  the  radiometer  in  the  nozzle 
is  making  sure  that  a  good  seal  exits  between  the  radiometer  and  the 
mating  nozzle  component.  For  the  radiometers  installed  at  the  nozzle 
throat  plane,  the  epoxy  resin  bond  acts  both  as  the  seal  and  as  the 
retention  system.  For  the  high  pressure  radiometer  in  the  exit  cone, 
the  seal  is  made  using  an  "0"  ring  which  is  compressed  when  the  proper 
torque  is  applied  to  the  radiometer.  The  installation  of  the  radi¬ 
ometers  in  the  nozzle  is  described  in  more  detail  in  the  following 
section. 


4. 2. 2. 3  Installation 

The  narrow  view  angle  radiometer  at  the  nozzle  throat  is 
located  in  the  MX  4926  ablative  component  as  shown  in  Figure  4-37. 

The  radiometer  was  bonded  in  the  ablative  component  using  Epon  828 
resin  with  a  Versamid  140  hardener.  The  primary  bond lines  are  along 
the  threaded  portion  of  the  barrel  and  along  the  face  and  sides  of  the 
1.50  inch  diameter  heat  sink.  These  bonds  were  cured  prior  to  the 
installation  of  the  ablative  component  in  the  steel  shell.  During  the 
installation  of  the  ablative  component,  the  purge  line,  sensor  leads, 
and  thermocouple  leads  were  routed  through  the  opening  in  the  steel 
shell. 


The  installation  of  the  high  pressure  narrow  view  angle  radiometer 
in  the  steel  shell  is  shown  in  Figure  4-38.  Since  this  radiometer 
was  sized  by  the  high  chamber  pressure  motor  application,  considerable 
modification  of  the  steel  shell  was  required  in  order  to  achieve  the 
proper  installation.  This  modification  was  dictated  by  the  cavity 
diameter,  by  the  required  recessed  depth  of  the  transducer,  and  by 
the  method  of  transducer  attachment.  Since  it  was  considered  impracti¬ 
cal  to  drill  the  cavity  in  the  ablative  component  prior  to  installing 
it  in  the  steel  shell,  then  the  cavity  diameter  must  not  exceed  the 
diameter  of  the  transducer.  This  dictated  that  the  steel  shell  ue 
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counterbored  a  maximum  of  approximately  0 . 20  inches ,  and  that  the 
upstream  edge  of  the  cavity  be  tapered.  After  allowing  approximately 
0.10  inches  for  the  "0"  ring  groove,  this  leave  0.20  inches  of  thread 
engagement  which  is  sufficient  for  the  maximum  expected  static 
pressure  at  an  expansion  ratio  of  3.  The  counterbore  diameter  is 
2.90  inches  which  allows  1/2  inch  on  the  radius  for  clearance  of  the 
gas  purge  line  and  thermocouple  and  sensor  leads.  The  counterbore 
of  the  steel  shell  was  machined  prior  to  installing  the  ablative  com¬ 
ponent.  In  addition,  a  small  centering  hole  was  machined.  However  the 
hole  was  not  made  in  the  steel  shell  and  ablative  component  until 
after  the  nozzle  assembly,  including  curing  of  the  bondlines.  After 
completing  the  necessary  machining  the  radiometer  and  the  "0"  ring 
were  installed  using  a  torque  of  approximately  10  ft-lbs. 


4. 2. 2. 4  Recording  Instrumentation 

Since  the  wires  leading  out  of  the  radiometer  are  both  copper 
(one  being  attached  to  the  center  of  the  constantan  disk  sensor  and 
the  other  to  the  copper  housing) ,  an  emf  is  measured  which  corresponds 
to  the  temperature  drop  across  the  sensor,  and,  thus,  no  reference 
temperature  junction  is  required.  This  deletion  of  the  reference 
temperature  junction  was  the  only  difference  between  the  recording 
circuit  for  the  radiometer  and  that  described  in  Section  4. 1.2. 4  for 
the  thermocouple.  The  output  from  the  thermocouple  which  was 
attached  to  the  copper  heat  sink  was  recorded  using  the  instrumentation 
circuit  described  in  Section  4.1.2'. 4.  Therefore,  standard  recording 
instrumentation  was  used  for  the  radiometer  as  well  as  the  thermo- 
coupled  plug. 
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SECTION  5 


DATA  REDUCTION 


The  basic  data  measured  by  the  two  heat  flux  measurement  sys¬ 
tems  described  in  Section  4  must  be  properly  reduced  to  define  the 
desired  total  and  radiative  fluxes.  These  procedures  are  described 
in  this  section.  Section  5.1  describes  the  procedure  for  relating 
the  in-depth  thermocouple  output  from  the  total  heat  flux  measure¬ 
ment  system  to  a  total  surface  heat  flux.  Section  5.2  describes 
the  procedure  for  relating  the  millivolt  output  of  the  narrow  angle 
directional  radiometer  to  a  total  surface  radiative  heat  flux. 


5.1  TOTAL  HEAT  FLUX  SYSTEM 

The  data  reduction  procedure  for  the  total  heat  flux  measure¬ 
ment  system  requires  the  conversion  of  data  which  define  the  tem¬ 
perature  histories  at  five  in-depth  locations  in  the  material  to 
total  heat  flux  at  the  material  surface.  The  details  of  this  pro¬ 
cedure  are  presented  in  this  section.  Section  5.1.1  discusses  the 
reduction  of  the  recorded  data  to  temperature  histories  and  Section 
5.1.2  discusses  the  conversion  of  these  temperature  histories  to 
total  heat  flux.  Finally,  Section  5.1.3  discusses  the  accuracy  of 
these  total  heat  fluxes  as  related  primarily  to  uncertainties  asso¬ 
ciated  with  the  data  reduction  process. 

5.1.1  Definition  of  Temperature  Histories 

The  definition  of  the  temperature  histories  is  straightfor¬ 
ward  and  employs  standard  data  reduction  techniques.  The  continu¬ 
ously  recorded  thermocouple  output,  in  the  form  of  millivolts  or  in 
the  form  of  galvanometer  deflection,  is  converted  to  temperature 
through  appropriate  calibration  curves.  (In  the  latter  caise,  an 
intermediate  step  to  convert  galvanometer  deflection  to  millivolts 
is  typically  required.)  The  calibration  curve  is  a  fit  of  the 
thermocouple  temperature  -  emf  characteristic  which  normally  falls 
within  ±1  percent  of  the  actual  thermocouple  calibration.  Chromel/ 
Alumel  thermocouple  wire  typically  agrees  very  closely  with  the  NBS 
standard  and  therefore  a  fit  of  this  standard  calibration  may  be 
used.  Tungsten  5  percent  rheniura/tungsten  26  percent  rhenium  has 
no  NBS  standard  and  therefore  a  fit  of  the  actual  calibration  for 
the  particular  wire  used  or  the  manufacturers  accepted  calibration 
is  employed.  The  manufacturers  actual  or  accepted  calibration 
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typically  goes  up  to  4200°F  only  but  the  fit  may  be  extended  to 
5400°F  based  On  the  results  of  Reference  4-12. 

The  recorded  thermocouple  output  and  the  reduced  data  are  sub¬ 
jected  to  close  scrutiny  to  relate  thermocouple  failure  to  surface 
recession  and  to  relate  the  temperature  at  thermocouple  failure  to 
surface  .temperature.  Any  anomolous  behavior  is  noted  for  considera¬ 
tion  in  the  subsequent  data  reduction  process. 


5.1.2  Definition  of  Total  Heat  Flux 


The  basic  experimental  results  in  the  form  of  in-depth  tem¬ 
perature  histories  discussed  in  the  previous  section  provide  the 
basis  for  the  conduction  solution  and  definition  of  total  heat  flux. 
The  sequence  of  events  in  this  procedure  is  1)  the  definition  of 
certain  energy  flux  terms,  surface  recession,  and  surface  tempera¬ 
ture,  based  primarily  on  the  measured  temperatures,  2)  the  defini¬ 
tion  of  other  necessary  flux  terms  related  primarily  to  the  sur¬ 
face  thermochemical  response,  and  then  3)  the  definition  of  the 
desired  total  heat  flux.  The  background  behind  this  procedure  and 
the  mechanics  of  accomplishing  it  are  presented  in  the  following 
paragraphs . 

The  energy  balance  at  the  surface  of  a  charring  ablator  is 
shown  in  the  sketch  below  (see  also  Figure  3-1) .  These  flux  terms 
vary  with  time  due  to  the  transient  response  of  the  material  and 
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due  to  surface  recession.  The  flux  terms  in  the  sketch  are  defined 
as  follows:  ' 


convective  heat  flux 

T 

%onv  =  puCH0<hr  ‘  Vs  +  puCMo£  (zie  "  ziw)hiW  which  reduces 
to  the  more  familiar  puCE  {hf  -  hw)fc  for  unity  Lewis 
number  and  equal  diffusion  coefficients 

radiation  heat  flux  from  the  free  stream  and/or  opposite  wall 
^rad  in  =  eeffffTeff  w^ere  thi-S  equation  is  more  descriptive 
than  computationally  solvable 


•  heat  flux  conducted  into  the  material 


‘cond 


-(*81 


•  radiation  heat  flux  from  the  wall 

^rad  out  ~  ewaTw  u. 

y 

•  energy  flux  due  to  pyrolysis  gas  convection 

qpyrol  “  j*1  ghg^w 

>  . 

•  energy  flux  due  to  char  convection  ^ 

i' 

^char  ~  "'c^c  ~~  *^pc^c 

w  w 

I  ■ 

•  energy  flux  due  to  convection  of  the  ablation  products 


* 


! 


•  convective  heat  flux  decrement  due  to  blowing  at  the  wall 

qblow  =  (X  *  ^77)qconv  where  qconv  is  defined  above,  net 

convective  flux  to  the  wall  'is*  q  -  q.  ,  '  = 

^conv  ^blow 


5  _  ^  ^conv 

•  energy  flux  due  to  mechanical  removal  of  condensed  phase 
qmech  ~  m*^* 

w 


The  sum  of  the  first  two  terms  above,  qconv  +  qra<j  in*  *-s  the  de¬ 
sired  total  heat  flux  which  is  y^iven  by* 


/ 


j 

qconv  +  qrad  in,  ~  qcond  +  qrad  out  +  q 


abl 


qLlow  "  qpyrol  qchar  +  qmech 


(5-1) 


The  data  reduction  procedure  encompasses  the  evaluation  of  the 
terms£  on  the  right  of  Equation  (5-1)  from  the  primary  measurements 
of  internal  temperature . 


The  first  step  is  an  internal  conjunction  calculation  using  the 
Charring  Material  Ablation  (CMA)  computer  program  (Reference  4-1) 
wherein. the  surface  recession  and  surface  temperature  histories 
are  specified  (Option  2) .  This  calculation  is  performed  iteratively 
by  successive  inputs  of  surface  recession  and  surface  temperature 
histories  until  good  agreement  between  the  measured  and  calculated  \ 
in-depth  temperature  histories  is  achieved.**  The  initial  inputs 
are  based  on  the  measured  in-depth  temperature  histories  and  the 
post- test  measured  recession.  The  final  inputs  for  whiOh  good 
agreement  between  in-depth  temperatures  is  achieved  define  the 
actual  surface  and  in-depth  response  of  the  material  in  terms  of: 


■k  v 

The  theoretical  details  of  the  surface  energy  balance  and  the 
terms  therein  are  presented  in  References  4-1  and  4-2. 

* 

This  same  procedure  has  been  used  successfully  in  a  reverse  fas- 
ion  wherein  the  surface  boundary  conditions  were  known  and  the 
material  thermal  conductivity  was  determined  by  requiring  a 
match  of  the  measured  and  calculated  in-depth  temperatures  (Ref¬ 
erences  4-5,  4-6,  and  5-1).  * 
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conduction  heat  flax,  <Jcond 

radiation  heat  flux  from  the  wall,  <3rad  out 

pyrolysis  off  gas  energy  flux,  c[pyroi 

char  energy  flux,  qcha^. 

surface  recession,  s 

surface  temperature ,  ( 

pyrolysis  off  gas  rate,  m_ 

char  removal  rate ,,  m 


Note  that  all  these  terms  are  a  function  of  time  through  the  test 
or  firing.  The  solution  of  the  energy  balance.  Equation  (5-1) ,  is 
now  reduced  to 


‘Jconv  +  grad  in  =  E  (known  q's)  +  gabl  +  %low  +  Wh  (5"2) 


The  above  results,  together  with  a  parametric  calculation  of 
the  surface  thermochemical  response  using  the  Aerotherm  Chemical* 
Equilibrium  (ACE)  computer  program  (Reference  4-3) ,  allow  the  cal¬ 
culation  of  the  other  required  terms .  For  the  materials  of  in¬ 
terest  herein,  diffusion  rate  controlled  (equilibrium)  surface 
thermochemistry  is  valid  and  is  required  in  the  calculations.  The 
program  output  provides  the  following  information  as  a  function  of 
surface  temperature,  dimensionless  off  gas  rate  B' ,  and  dimension¬ 
less  char  removal  rate  B' :  g 

c 


i 

dimensionless  gas  phase  ablation  rate,  B  =  (  VA/puC 

/  a  /  w.  M 

dimensionless  condensed  phase  removal  rate , \ B*  =  m*/puCM 

wall  enthalpy,  hw 

recovery  enthalpy,  h  i  T 

chemical  enthalpy  potential,  E  (Z^e  -  Z^w)h^  . 

transport  properties,  specifically  Lewis  number,  Le,  where 


<VCH  '  Lei/’ 


Plots  of  this  output  in  the  form  shown  below  allow  the  calcu¬ 
lation  of  the  remaining  flux  terms. 
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The  energy  flux  due  to  convection  of  the  ablation  products 
from  the  surface  is  given  by  , 


q  ,  ,  =  (pv)  h  =  (m 
^abl  K  w  w  —  g 


mc) 


B* 

a 


B* 

c 


w 


(5-3) 


•  • 

where  rrig  and  mc  are  available  from  the  conduction  solution  and  the 
other  two  terms  are  available  from  the  surface  thermochemistry  so¬ 
lution  (see  the  above  plot) .  The  energy  flux  due  to  mechanical 
removal,  if  appropriate,  is  given  by 


b; 

Wh  =  m*h*  =  (mg  +  mc>  B1  +  B-  hw  (5‘ 

<  ^  g  c 

./ 

where  the  terms  are  available  as  above,  except  that  the  condensed 
phase  enthalpy  is  a  function  of  the  surface  temperature  cnly. 


The  calculation  of  the  convective  heat  flux  decrement  due  to 
blowing  requires  the  definition  of  the  heat  and  mass  transfer  coef¬ 
ficients.  The  dimensionless  char  removal  rate  B,1  is  first  deter-  . 
mined  from  a  plot  of  the  form  above.  The  mass  transfer  coefficient 
is  ehen  calculated  from  the  defining  relation  between  nj,  and  B^ 
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and  the  heat  transfer  coefficient  calculated  from  Colburn's  analogy 


„  puCjj  =  puC^Le“2/3 


(5-6) 


The  blowing  reduction  flux  is  given  by 


qblow  (E  -  eC  \  1)qconv  (X  '  a;  \  1)[P"C»0<hr  ‘ 


+  »uCM  Ei,Zie  -  Zi„» 


■I“]  <5-' 


where  c  =  2Xfpv)w/puCM  .  Note,  however,  that  puCH  and  puCMq  are 

4  ®  r  r 

not  Known,  rather  puC„  =  — — — ;  puCH  and  puCM  =  — -  puCM 

H  eC  -  1  °  M  eC  -  1  ° 

only  are  known.  Some,  simple*  algebraic  manipulation,  however,  yields 
a  solvable  expression  for  £; 


e1*  -  1  = 


2A (pv) 


=  2XB' 


(5-8) 


from  which  puCy  and  puCmq  can  be  calculated  and  Equation  (5-7) 
solved.  Note  that  the  chemical  enthalpy  term  in  Equation  (5-7) 
is  also  determined  from  a  plot  of  the  form  above. 


With  all  flux  terms  on  the  right  of  Equation  (5-1)  now  calcu¬ 
lated,  the  total  heat  flux  is  available  by  the  summation  of  these 
terms. 


5.1.3  Accuracy  of  Total  Heat  Flux 

The  sources  of  error  in  the  calculation  of  total  heat  flux 
by  the  method  described  in  Section  5.1.2  are  uncertainties  in  the 
primary  temperature  measurement,  in  the  internal  conduction  solu¬ 
tion,  and  in  the  surface  equilibrium  calculation.  The  severity  of 
errors  introduced  due  to  the  mathematical  modeling  of  physical 
events  is  difficult  to  assess.  Based  on  past  experience  with  the 
analytical  technique  utilized  in  the  heat  flux  analysis,  however, 
the  accuracy  of  calculated  results  can  be  approximated.  Sources 
of  uncertainty  in  the  temperature  measurements,  the  conduction 
solution,  and  the  thermochemical  model  are  discussed  below. 


The  overall  uncertainties  in  the  measured  in-depth  temperatures 
were  presented  in  Section  4.4.  Their  estimated  values  at  high  tem¬ 
perature  are  -75°R  and  +15°R  for  carbon  phenolic  and  -125°R  and 
+150°R  for  silica  phenolic.  These  potential  errors,  of  course,  pro¬ 
pagate  through  the  subsequent  data  reduction  procedures. 


The  accuracy  of  the  internal  conduction  solution  is  dependent 
upon  the  material  response  model  employed,  the  thermophysical  and 
degradation  properties  input,  the  experimental  in-depth  temperature 
response,  and  the  sensitivity  of  the  iterative  procedure.  Experi¬ 
ence  has  shown  that,  for  the  materials  utilized  in  this  program, 
the  response  model  and  material  properties  give  accurate  results. 
Also,  the  iterative  procedure  is  quite  sensitive  to  the  input  values 
of  surface  recession  and  surface  temperature.  Therefore,  the  final 
values  of  the  qcond»  ^char'  ^pyrol'  ^rad  opt'  s'  an<^  Tw  histories 
are  expected  to  be  within  the  accuracy  limitations  of  the  in-depth 
temperature  determinations  on  which  the  iteration  is  based. 


The  sensitivity  of  the  conduction  solution  to  the  in-depth  tem¬ 
perature  uncertainties  was  studied  by  parametrically  varying  the 
surface  temperature  and  surface  recession  rates  about  the  values 
calculated  in  the  NASA/Aerojet  260-SL-3  nozzle  prediction  presented 
in  Reference  4-6.  Figure  5-1  presents  a  graphical  display  of  some 
important  results  of  the  parametric  study  at  30  seconds  through  the 
firing.  It  defines  the  dependence  of  the  predicted  in-depth  tem¬ 
perature  profile  on  a  ±200°R  variation  of  surface  temperature.  Con¬ 
versely,  the  effect  of  in-depth  thermal  response  uncertainties  on 
surface  temperature  accuracy  can  be  evaluated  for  a  given  surface 
location  uncertainty.  Note  that  the  range  of  uncertainty  is  rep¬ 
resented  by  the  shaded  region.  Also,  the  uncertainty  values  chosen 
for  the  analysis  were  extreme  and  are  more  illustrative  than  real. 

It  is  apparent  that  the  uncertainty  in  the  surface  location  and 
in-depth  temperature  has  a  strong  effect  on  the  surface  tempera¬ 
ture.  For  the  combined  extremes  considered,  ±0.025  inch  and 
±200°R,  the  surface  temperature  is  within  ±400UR  of  the  true  value 
in  the  case  considered.  More  realistic  uncertainty  values  (Sec¬ 
tion  4.4  and  above)  yield  approximate  surface  temperature 
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uncertainties  of  -125°R  and  +25  R  for  carbon  phenolic  and  -200°R 
>v  and  +75°R  for  silica  phenolic  in  the  char  motor  nozzle. 


The  parametric  calculations  of  the  260-SL-3  nozzle  response 
were  also  utilized  to  assess  the  effects  of  conduction  solution 
uncertainties  on  the  calculated  total  heat  flux.  For  each  CMA  con¬ 
duction  solution  the  total  heat  flux  was  calculated  so  that  the 
effects  of  surface  temperature  or  recession  errors  could  be  evalu¬ 
ated.  The  variation  of  surface  recession  had  only  a  slight  effect 
on  the  various  heat  flux  terms;  however,  the  total  heat  flux  was 
quite  sensitive  to  the  variation  of  surface  temperature .  Figure 
5-2  quantitatively  relates  surface  temperature  uncertainty  to  a 
corresponding  uncertainty  in  total  heat  flux.  The  sensitivity  is 
primarily  due  to  the  dominance  of  the  radiation  and  conduction 
heat  flax  terms  (both  strongly  dependent  on  surface  temperature) 
in  the  surface  energy  balance.  A  ±4  percent  (±200°R)  uncertainty 
in  surface  temperature  represents  a  -13  percent  and  +18  percent 
uncertainty  in  total  incident  heat  flux.  If  the  results  of  Figure 
5-2  are  applied  to  the  more  accurate  uncertainties  of  the  previous 
paragraph,  the  resultant  error  band  in  total  heat  flux  is  -8  per¬ 
cent  and  +2  percent  for  carbon  phenolic  and  -15  percent  and  +8 
percent  for  silica  phenolic.  Note  that  the  results  for  silica 
phenolic  are  estimates  since  Figure  5-2  is  strictly  valid  only 
for  carbon  phenolic  response. 
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FICURE  5-2  SENSITIVITY  OF  TOTAL  HEAT  FLUX  TO 
SURFACE  TEMPERATURE  UNCERTAINTY 


The  above  uncertainties  assume  that  the  heat  flux  terms  based 
on  the  surface  thermochemical  response  (q^^f  qblow'  ^mech^  are 
completely  accurate.  Realistically  this  is  not  the  case,  and  ac¬ 
count  must  be  taken  of  their  contribution  to  the  error  in  total 
heat  flux.  These  terms  are  small,  however  (e.g.  see  the  table  of 
Section  5.1.2),  and  therefore  the  impact  of  any  errors  in  their 
values  on  the  total  heat  flux  is  also  small. 


Summarizing  the  various  sources  of  error,  the  overall  accuracy 
of  the  determined  total  heat  flux  is  estimated  in  the  table  below. 


Carbon  Phenolic  Silica  Phenolic 

-11  percent  -20  percent 

+  5  percent  +13  percent 


The  error  analysis  was  done  on  a  conservative  analytical  basis  and, 
therefore,  these  values  are  felt  to  be  outside  extremes;  the  actual 
accuracy  is  probably  better  them  indicated. 
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5.2  RADIATIVE  HEAT  FLUX  SYSTEM 

The  data  reduction  procedure  requires  the  conversion  of  the 
emf  output  from  the  sensor  to  a  total  hemispherical  radiative  heat 
flux.  For  this  conversion  to  be  meaningful,  one  must  understand 
the  simplifying  assumptions  which  are  made,  and  these  are  dis¬ 
cussed  in  Section  5.2.1.  The  accuracy  of  the  measured  radiative 
heat  flux  is  discussed  in  Section  5.2.2. 


5.2.1  Definition  of  Radiative  Heat  Flux 


The  factors  considered  in  selecting  a  narrow  view  angle  radi¬ 
ometer  to  measure  the  total  radiative  heat  flux  to  the  nozzle  sur¬ 
face  were  defined  in  Section  3.2.  The  topic  to  be  considered  here 
is  the  assumptions  which  were  made  to  relate  the  output  from  the 
narrow  view  angle  radiometer  to  a  total  hemispherical  radiative 
heat  flux. 


The  total  radiative  fliu?  to  an  elemental  area  on  the  nozzle 
surface  is  a  complicated  function  of  the  absorbing  and  scattering 
characteristics  of  the  particle  laden  combustion  products.  It 
is!  complicated  further  by  the  geometry  of  the  internal  flow  con¬ 
tour  (Figure  5-3) .  Even  though  the  cross  sectional  area  of  the 
nozzle  throat  and  exit  cone  is  axisymmetric,  a  surface  elemental 
area  has  a  field  of  view  which  varies  in  a  complicated  manner  with 
the  axial  location.  Since  the  basic  features  of  the  governing 
radiative  transport  equation  are  more  easily  understood  by  consid¬ 
ering  a  simplified  geometrical  model,  the  model  considered  here  is 
represented  by  two  planar  parallel  slabs  separated  by  an  absorbing, 
scattering  media.  The  slabs  simulate  the  nozzle  surface,  and  the 
media  is  representative  of  the  particle  laden  combustion  products 
for  a  solid  propellant.  A  schematic  of  this  simplified  model  is 
shown  in  Figure  5-4.  The  physical  model  for  developing  the  govern¬ 
ing  radiative  transport  equation  id  shown  in  Figure  5-5.  Note  that 
the  temperature  of  the  media  is  assumed  to  be  only  a  function  of  y 
and  that  the  total  intensity  of  radiation,  I^(y,6),  is  a  result  of 
the  contributions  in  the  positive  y  direction,  lt(y,0),  and  in  the 
negative  y  direction,  lj[(y,0).  If  the  elemental  surface  area  is 
on  surface  2,  the  total  intensity  of  radiation  is  given  by  I^(y,6). 
If  one  assumes  that  the  media  is  stationary  and  that  the  particle 
scattering  is  isotropic  and  coherent,  then  the  governing  radiative 
transport  equation  is  (Reference  5-2) 


+  +  ~tA/vj 

IA(toA'*j)  =  1  <°'^e  ’  + 


Jo  A 

*  /  *r[ 


KAebA(t) 


T  ga  <«] 


-(T0A-t)/Hdt 
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where 


e^  -  Planck's  monochromatic  blackbody  emissive  power 


G^  -  scattering  function 

-  absorption  coefficient 

-  scattering  coefficient 

-  extinction  coefficient 
A  -  wavelength 

p  -  cosO 

xQ^  -  optical  depth 


The  first  term  in  Equation  (5-8)  represents  the  energy  emitted  by 
surface  1  which  is  incident  on  surface  2.  Note  that  the  energy 
leaving  surface  1  has  been  attenuated  due  to  the  absorption  and 
scattering  of  the  medium.  The  integral  term  in  Equation  (5—8) 
represents  the  effect  of  the  emitting  and  scattering  medium.  !rhe 
variables  eb^,  G^,  8^ ,  and  xQ^  are  given  by 


Because  the  intensity  of  radiation  is  under  the  integral  sign  in 
Equation  (5-10),  Equation  (5-8)  is  difficult  to  solve. 


The  total  flux  absorbed  by  the  sensor  of  the  narrow  view  angle 
radiometer  is  obtained  by  integrating  Equation  (5-8)  with  respect 
to  the  appropriate  angle  and  with  respect  to  wavelength.  This  is 
represented  by 


/  'xU1! 


(ToA,w)c°s0dwdX 
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where  .  ''  j 

w  -  solid  angle 

«  / 

w^  -  the  integration  limit  on  the  solid  angle 
eA  -  emissivity  of  sensor  coating 
t’  -  transmissivity  of  sapphire  lens 


When  the  narrow  view  angle  radiometer  is  exposed  to t the  calibration 

source,  the  radiation  flux  absorbed  by  the  sensor,  qr  ,  is  given  by 
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Since  the  total  hemispherical  radiative  flux  from  a  blackbody  source, 
is  given  by 

00 

qr  -  »  / xbxdX  (5-15) 

cb  J 

A=0  / 


then,  if  one  assumes  that  ex  and  tJ  are  independent  of  wavelength 
Equations  (5-14)  and  (5-15)  can  be  combined  to  give 
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Note  that  the  ratio  of  the  integral  term  to  ir  is  the  view  factor  , 
that  the  sensor  has  of  the  radiation  source. 

*  i 

i v  1  5 

During  the  calibration  of  the  narrow,  view  angle  radiometer ,  an 
emf  is  measured,  and  this  emf  is  related  'to  a  total  hemispherical 
radiative  heat  flux  which  is  measured  using  a  calibrated  standard". 
As  a  first  approximation,  this  calibration  procedure  has  the  effect 
of  measuring  qr  which  is  then  related  to  the  total  blackbody  hemi¬ 
spherical  flux, qr  ,  using  Equation  (5-16). 


In  the  data  reduction  procedure,  the  totar  hemispherical  radia¬ 
tive  flux  is  rjeport^d .  However,  this  flux  can  be  related  to  a  par¬ 
ticular  radiation  sohrce  only  if  certain  restrictive  assumptions 
are  made.  For  example,  if  one  assumes  that  1)  the  radiation  emitted 
by  the  particle  laden  combustion  products  is  negligible,  2)  the 
scattering  coefficient  for  these  products  is'  zero,  and  3)  the  oppo¬ 
site  wall  emits  as  a  constant  teng>erature  gray  body,  then  the  re¬ 
ported  total  hemispherical  heat  flux  is,  a  true  indication  of  the 
radiative  heat  flux  from/the  opposite  wall.  Another. set  of^qssump- 
tions  which  isolate  the  Radiation  source  are  1)  the  Scattering  coef¬ 
ficient  for  the  combustion  products  is  zero,  ,2)  the  radiation  inten¬ 
sity  is  Independent  of  direction,  and  3)  the  particle  laden  combus¬ 
tion  products  are' optically  thick.  These  latter  sft  of  assumptions 
are  more  realistic  for  calculating  the  radiative  flux  from  the  com¬ 
bustion  products  of  a  solid  rocket  motor,  and,  when  usfd,  the  total 
heat  flux  emitted  by  -the  slab  of  combustion  products,  qr  ,  is 
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where  E, (t)  is  an  exponential  integral  function.  If  the  monochro¬ 
matic  form  of  this  heat  flux  is  se£.  equal  to 
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where  e„,  is  a  hemispherical  monochromatic  emittance.  Then,  assum¬ 
ing  a  constant  temperature.  Equations  (5-17)  and-  (5-18)  can  be  com¬ 
bined  to  give 
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For  'an  optically  thick  medium,  the  exponential  integral  in  Equation 
(5-19)  is  zero  which  results  in  the  emittance  being  equal  to  one. 

In  this  case,  the  combustion  products  radiate  as  a  blackbody.  For 
finite  values  of  the  monochromatic  emittance  which  are  independent 
/of  wavelength,  the  combustion  products  radiate  as  a  gray  body.  In 
either  case,  the  measured  radiative  heat  flux  is  a  true  indication 
of  the  amount  oi  radiation  being  emitted  by  the  combustion  products 
(within  the  assumption  of  a  constant  temperature) . 

Since  a  detailed  partitioning  of  the  measured  total  hemispheri¬ 
cal  radiative  flux  into  components  from  the  various  radiation  sources 
in,  a  ^Jolid  propellant  environment  is  beyond  the  scope  of  the  present 
study,  the  data  reduction  procedure  is  simply  a  conversion  of  the 
measured  emf  output  from  the  narrow  view  angle  radiometer  to  a  total 
hemispherical  radiative  flux.  This  conversion  is  performed  using 
the  calibration  supplied  by  the  manufacturer.  As  discussed  above, 
this  is  strictly  valid  only  for  the  two  sets  of  limiting  assumptions 
presented  above  -  essentially  either  a  transparent  stream  or  an 
opaque  stream. 


i 

5.2.2  Accuracy  of  Radiative  Heat  Flux 

Based  on  the  assumptions  that  the  c&mbustion  products 
viewed  by  the  radiometer  1)  are  an  optically  thick  radiation  source 
which  is  at  constant  temperature  and  2)  have  an  emissivity  which  is 
independent  of  wavelength,  then  the  accuracy  of  the  hemispherical 
radiation  heat  flux  is  identical^io  that  given  in  Table  4-5  (page 
4-51)  for  the  heat  flux  absorbed  by  the  sensor.  The  accuracy  of 
the  radiometer  system  stated  in  this  table  was  +2.5/L3  percent. ^ 


o 


5-16 


'J 


I 

('- 


REFERENCES 


‘Baker,  D.  L. ,  Wool,  M.  R. ,  and  Schaefer,  J.  W. :  A  Dynamic 
Technique  for  Determining  the  Thermal  Conductivity  of  Charring 
Materials.  Paper  presented  at  Eighth  Conference  on  Thermal 
Conductivity. 

Sparrow,  E.  M.  and  Cess,  R.  D. :  Radiation  Heat  Transfer. 
Breaks /Cole  Publishing  Co.,  August  1967  (Second  Printing). 

♦ 


/ 


5-17 


/ 


SECTION  6 

TEST  EVALUATION  PROGRAM 


The  total  heat  flux  and  narrow  view  angle  radiometer  measure¬ 
ment  systems  were  evaluated  under  laboratory  test  conditions  typi¬ 
cal  of  a  rocket  nozzle  application  in  the  Phase  II  test  program. 

The  purpose  of  this  program  was  to  demonstrate  the  measurement  sys¬ 
tems  prior  to  subjecting  them  to  the  char  motor  firings  at  AFRPL. 
The  Phase  II  test  program  is  described  in  Section  6.1.  The  results 
of  this  program  are  presented  in  Section  6.2  for  the  total  heat 
flux  system  and  in  Section  6.3  for  the  radiometer  system. 


6.1  PHASE  II  TEST  EVALUATION  PROGRAM 

i 

The  objectives  of  the  Phase  II  test  evaluation  program  were  to 
verify  the  measurement  system  design  and  fabrication  technique  and 
to  verify  the  data  reduction  procedure  and  firmly  define  the  total 
heat  flux  accuracy.  The  tests  were  performed  under  laboratory  con¬ 
ditions  typical  of  a  solid  rocket  nozzle  and  for  which  the  total 
heat  flux  and  other  important  variables  which  are  calculated  in 
the  measurement  system  data  reduction  procedure  are  known. 


The  Aerotherm  arc  plasma  generator  (APG) ,  Figure  6-1,  was  used 
to  provide  the  solid  propellant  exhaust  simulation,  and  a  two- 
dimensional  nozzle  test  section.  Figure  6-2,  was  used  as  the  rocket 
nozzle.  In  the  APG,  energy  is  added  electrically  to  the  test  gas 
by  a  continuous  arc  discharge.  Two  test  gases  were  used  in  this 
test  program;  both  closely  duplicate  the  specific  heat  of  typical 
aluminized  solid  propellants  but  one  is  chemically  inert  (Mixture  5) 
and  the  other  is  chemically  reactive  (Mixture  4).  In  Mixture  4, 
the  oxidation  potential  of  typical  solid  propellants  is  duplicated 
(Reference  4-7) .  The  chemical  compositions  of  these  test  ga&es  are 
presented  in  the  table  below  and  further  background  on  their  selec¬ 
tion  is  presented  in  References  4-5  through  4-7.  For  Mixture  5, 
no  surface  recession  occurs  for  MX  4926  carbon  phenolic;  this  allows 
the  evaluation  of  the  measurements  system  for  the  simpler  case  of 
no  surface  recession  and  the  identification  of  any  basic  problems 
that  might  have  otherwise  been  obscured.  For  Mixture  4,  the  mate-, 
rial  response  closely  duplicates  that  in  an  actual  nozzle.  This 
provides  the  complete  evaluation  of  the  measurement  system  at  condi¬ 
tions  typical  of  its  application. 
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ELEMENTAL  MASS  FRACTIONS 

Mixture  4  Mixture  5 

Helium  0.228  0.224 

Nitrogen  0.619  0.776 

Oxygen  0.153 


In  addition  to  providing  $  test  section  which  has  a  relatively 
constant  heat  flux,  the  two-dimensional  nozzle  eliminates  the  effect 
of  surface  curvature  in  the  evaluation  of  the  thermocouple  plug. 
Since  these  plugs -will  be  used  in  nozzle  which  have  a  local  radius 
which  is  lirge  compared  to  the  radius  of  the  plug,  the  tWo-dimen- 
sional  nozzle  provides  a  better  geometrical  simulation  than  an 
axisymmetric  nozzle  with  the  same  throat  area.  The  test  sections 
were  fabricated  by  rocket  nozzle  techniques  to  insure  that  the  fab¬ 
rication  process  was  not  a  variable.  The  measurement  system  thermo- 
coupled  plug  or  radiometer  was  located  at  the  throat  in  all  the 
evaluation  tests  since  the  test  conditions  can  be  more  accurately 
defined  at  this  location.  The  layup  angle  in  all  test  sections  was 
60°,  equal  to  or  close  to  that  for  the  measurement  system  locations 
in  the  char  motor  nozzle.  The  test  section  opposite  the  measure¬ 
ment  system  test  section  waa  the  same  material  as  the  measurement 
system  test  section  or  was  water-cooled  copper.  In  the  former  case 
the  radiation  flux  to  the  measurement  system  (qr£{j  in)  was  finite 
and  in  the  latter  case  it  was  zero.  This  provided  another  boundary 
condition  variable  which  'allowed  the  identification  of  any  basic 
problems  which  might  have  otherwise  been  obscured.  The  plug  design, 
fabrication  technique,  and  installation  procedure  used  in  these 
tests  were  identical  to  that  used  in  the  char  motor  except  for  neces¬ 
sary  changes  in  plug  length  and  thermocouple  locations.  These  items 
are  discussed  in  more  detail  in  Sections  6.2  and  6.3. 


The  nominal  test  matrix  for  the  Phase  II  test  evaluation  pro¬ 
gram  is  presented  in  Table  6-1.  All  of  the  test  variables  discus¬ 
sed  above  were  included  in  the  test  matrix  and  two  sets  of  dupli¬ 
cate  tests  were  included  to  check  repeatability  of  results.  Note 
that  since  surface  recession  would  be  expected  to  occur  for  silica 
phenolic  in  both  the  inert  and  reactive  environments ,  only  the 
more  complete  simulation  case.  Mixture  4,  was  considered  for  this 
material.  In  addition  to  the  tests  indicated  in  the  table,  cali¬ 
bration  tests  were  performed  at  each  test  condition  to  measure  the 
cold  wall  hr  at  flux  at  the  nozzle  throat.  For  those  tests,  the 
material  test  section  was  replaced  by  a  water-cooled  calorimeter 
test  section  for  measurement  of  convective  flux. 


The  chamber  conditions  achieved  during  the  APG  tests  together 
with  the  material  ablative  performance  are  summarized  in  Table  6-2. 
By  comparing  the  results  of  Tests  1317  and  1318  which  used  carbon 
phenolic  test  sections  and  Tests  1320  and  1321  which  used  silica 
phenolic  test  sections,  one  notes  that  approximately  identical 
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chamber  conditions  produced  identical  material  ablative  performance, 
thus  indicating  repeatability  of  performance.  Also,  one  notes  that 
negative  material  recession  was  obtained  for  Test  1319,  and  since 
the  environment  was  chemically  inert,  this  negative  recession  was 
due  to  char  swell.  The  radiometer  was  exposed  to  the  APG  environ¬ 
ment  in  Tests  1324  and  1325.  Originally,  these  tests  were  supposed 
to  duplicate  the  nominal  test  conditions  shown  by  Tests  8  and  9  in 
Table  6-1.  However,  the  radiometer  output  for  Test  1324  was  lost 
due  to  a  malfunction  of  the  instrumentation  recording  equipment. 
Since  radiometer  results  from  a  test  with  a  hot  carbon  cloth  phe¬ 
nolic  opposite  wall  would  be  of  more  value  than  results  for  a  cold 
copper  opposite  wall,  the  test  conditions  for  Testri325  were  altered 
to  repeat  those  of  Test  1324.  The  cold  opposite  wall  test  was  not 
conducted . 


The  thermocouple  and  emf  data  from  the 
heat  flux  systems  are  presented  in  Sections 
tively.  These  sections  present  the  reduced 
of  the  measurement  systems. 


total  and  radiative 
6.2  and  6.3,  respec- 
data  and  an  evaluation 


6.2  TOTAL  HEAT  FLUX  SYSTEM 

The  total  heat  flux  system  was  evaluated  in  seven  APG  tests. 
The  design  and  installation  of  the  sensor  for  these  tests  is  de¬ 
scribed  in  Section  6.2.1.  The  test  results,  including  both  the 
measured  data  and  the  reduced  data  using  the  data  reduction  tech¬ 
niques  described  in  Section  5,  are  presented  in  Section  6.2.2 

6.2.1  Design  and  Installation 

The  thermocoupled  plug  design  employed  in  the  test  evalua¬ 
tion  program  was  shown  in  Figure  4-2,  and  the  installation  in  the 
two-dimensional  test  section  was  shown  schematically  in  Figure 
6-2.  The  thermocouples  were  installed  at  the  nominal  depth  below 
the  surface  of  0.050,  0.125,  0.200,  0.300,  and  0.400  inches,  and 
the  thermocouple  plug  was  0.800  inches  long,  excluding  the  strain 
relief  pins.  Since  the  test  section  had  a  one-inch  depth  at  the 
nozzle  throat,  this  left  a  0.20  inch  cavity  at  the  backside  of  the 
plug.  This  cavity,  which  contained  the  strain  relief  pins  and  the 
thermocouple  lead  wire  connections,  was  potted  after  installing  the 
plug  in  the  test  section.  As  showp  in  Figure  6-2,  the  instrumented 
test  section  formed  one  half  of  the  APG  nozzle,  and  a  dummy  model, 
usually  fabricated  from  the  same  material,  formed  the  other  half. 


After  installing  the  instrumented  test  section  in  the  APG  noz¬ 
zle,  the  thermocouple  lead  wires  were  connected  to  the  measuring 
circuit  as  described  in  Section  5.  In  addition,  an  optical  pyrom¬ 
eter  was  focused  on  the  throat  of  the  instrumented  test  section. 
This  pryometer  measured  the  brightness  surface  temperature  at  a 
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wavelength  of  0.8  microns.  This  temperature  was  converted  to  a 
true  surface  temperature  using  the  material  eraissivity. 


6.2.2  Test  Results 


The  recorded  data  from  the  Phase  II  tests  of  the  total  heat 
flux  sensors  which  are  pertinent  to  the  measurement  of  this  flux 
include  thermocouple  and  pyrometer  data  and  pre-  and  post-test 
model  geometrical  measurements.  In  addition,  sufficient  parameters 
affecting  the  performance  of  the  APG  are  monitored  so  that  the 
properties  of  the  high  temperature  gases  can  be  defined.  The  gas 
properties  of  importance  include  enthalpy  and  pressure,  and  these 
have  been  reported  in  Table  6-2  of  Section  6.1.  Also,  reported 
'  in  this  table  is  the  total  surface  recession  and  the  average  sur¬ 
face  recession  rate  for  each  test.  In  the  following  paragraphs, 
an  example  of  the  in-depth  thermocouple  data  from  a  heat  flux  sen-  , 
sor  is  presented.  This  is  followed  by  a  presentation  of  the  com¬ 
parison  between  the  calculated  surface  temperature  (calculated  using 
the  in-depth  thermocouple  and  surface  recession  data)  and  the  mea¬ 
sured  temperatures  (defined  by  the  optical  pyrometer) .  Finally, 
the  total  heat  flux  obtained  from  the  data  reduction  procedure 
is  presented  for  each  of  the  tests,  and  this  heat  flux  is  converted 
to  a  convective  heat  transfer  coefficient.  The  intermediate  data 
reduction  steps  are  described  qualitatively.  These  steps  have  been 
described  quantitatively  in  Section  5  and  are  described  further  in 
Section  7  by  using  a  heat  flux  sensor  in  the  char  motor  firing  as 
an  example. 

Figure  6-3  presents  the  measured  in-depth  temperatures  as  a 
function  of  test  time  for  Test  1318  and  the  surface  temperature* 
calculated  using  the  in-depth  thermocouple  data  and  the  measured 
surface  recession.  As  described  in  Section  5.1.2,  the  surface  tem¬ 
perature  is  determined  by  an  iterative  calculation  using  the  CMA 
computer  program.  The  predicted  in-depth  temperatures  for  the 
final  iteration  on  the  surface  recession  rate  and  surface  tempera¬ 
ture  are  also  presented  in  Figure  6-3.  Note  in  this  figure  that 
the  nominal  thermocouple  depths  are  presented.  However,  the  actual 
depths  used  in  the  analysis  were  defined  using  X-ray  photographs. 


The  determined  surface  temperature  for  the  carbon  phenolic 
tests  are  compared  with  those  measured  using  an  optical  pyrometer 
in  Figures  6-4  through  6-6.  Similar  results  are  presented  in 


* 

The  surface  temperature  calculated  using  the  m-depth  thermcouple 
data  is  referred  to  as  the  determined  surface  temperature. 
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TEMPERATURE 


Figures  6-7  and  .6-8  for  the  silica  phenolic  tests-  These  figures 
illustrate  the  good  agreement  obtained  between  the  determined  and 
measured  surface  temperatures. 


The  final  iterated  values  of  the  determined  surface  tempera¬ 
tures  and  the^ surface  recession  rate^  together  with  the  in-depth 
pyrolypis  gas  generation  rate  (defined  by  the  CMA  program)  form 
the  basis  for  calculating  the  total  heat  flux.  The  ACE  computer 
program  generates  the  thermbchemical  information  used  in  the  data 
reduction  procedure.'  Because  the  optical  depth  of  the  gases  at 
the  nozzle  throat  is  very  small,,  the  gases  do  not  have  to  be  con¬ 
sidered  in  the  radiation  model,  and,  thus,  the  radiative  heat  flux 
to  the  sensor  surface  is  only  that  being  'emitted  by  the  opposing 
wall.  Therefore,  the  convective  heab  flux  is  easily  obtalined  from 
the  total  heat  flux.  For  the  heat  flux  sensor  tests  in  the  Phase 
II  program,  both  the  measured  total  and  convective  fluxes  are ‘pre¬ 
sented  in  Figures  6-9  through  6-11.  The  fcaih  item  to  note  in  these 
figures  is  that  the  measured  heat  fluxes  lor  carbon  cloth  phenolic  f 
in  Tests  1317  and  1318  are* essentially  identical.  Since  the  arc 
operating  characteristics  for  these  tests  were  approximately  iden¬ 
tical,  these"  measured  results  demonstrate  the  “repeatability  of  the 
measurement|  system.  In  general,  this  same  conclusion  is  demon¬ 
strated  in  Tests  1320  and  1321,  for  the  silica  cloth  phenolic. 
However,  the  measured  results  for  these  latter  tests  show  moife  V 
variance  than  for  the  carbon  cloth  phenolic  tests.  Since  the  mea¬ 
surement  system  accuracy  for  silica  cloth  phenolic  is  less  than  for 
carbon  cloth  phenolic  (Section  4)  ,  this  larger1  variance  ijs  to  be 
expected.  . 


The  convective  heat  fluxes  shown  in  Figures  6-9  through  6-11  % 
were  converted  to  a  convective  heat^jtransfer  coefficient  by,  dividing 
the  difference  between  the  free  stream  and  wall* enthalpies .  These 
convective  heat  transfer  coefficients  are  shown  in  Figures  6-12  and 
6-13  as  .a  function  of  test  time.  Also  shown  in  these  f igure^s~-is__^ 
the  predicted  variation  of  the  cohvective  heat  transfer  coefficient 
based  on  the  0.8  power  of  the  pressure  ratio.  Note  that  both  the 
predicted,  and  measured  convective  heat  transfer  coefficients  were 
assumed  identical  at  the  instant  the  measured  value  was  a  maximum. 
One  notes  from  these  comparisons  that  both  the  measured  and  pre¬ 
dicted  coefficients  decrease  with  test  duration.  However,  the 
measured  values  decrease  more  rapidly. 


RADIATIVE  HEAT  FLUX  SYSTEM 


The  radiative  heat  flux  system  was  evaluated  in  two  APG  tests. 
A  brief  description  of  the  design  and-  installation  of  this  system 
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is  presented  in  Section  6.3.1.  The  test  results  from  these  tests 
are  presented  in  Section  6.3.2. 


6.3.1  Design  and  Installation 

The  narrow  view  angle  radiometer  used  for  the  APG  tests  was 
identical  to  the  one  described  in  Section  4.  This  radiometer  was 
installed  in  a  carbon  cloth  phenolic  test  section  similar  to  that 
used  for  the  heat  flux  sensor.  However,  the  depth  of  the  test  sec¬ 
tion  had  to  be  increased  in  order  to  prevent  overheating  of  the 
aperture  of  the  radiometer.  JThe  modified  test  section  is  shown  in 
Figure  6-14.  Note  that  the  maximum  depth  is  1.65  inches  which  means 
the  aperture  of  the  radiometer  is  0.45  inches  below  the. initial  sur¬ 
face  of  the  test  section.  The  flange  and  threaded  section  of  the 
radiometer  were  bonded  to  the  test  section.  Note  that  the  threads 
on  the  radiometer  and  the  groove  in  the  test  section  insure  an  ef¬ 
fective  seal  of  the  possible  gas  path.  After  installing  the  instru¬ 
mented  test  section  in  the  APG,  the  thermocouple  and  sensor  lead 
wires  were  connected  to  a  recording  circuit  similar  to  the  one  de¬ 
scribed  in  Section  4.  In  addition,  the  purge  line  was  connected  to 
a- 1200  psia  nitrogen  source.  (The  flow  system  for  this  purge  sys¬ 
tem  was  similar  to  that  shown  in  Figure  7-4  for  the  first  char  motor 
firing.)  The  sonic  orifice  in  the  flow  system  was  sized  to  give  a 
purge  flow  rate  of  0.007  ILs/sec.  With  this  flow  rate,  the  pre¬ 
dicted  cavity  temperature  using  the  method  outlined  in  Section  4 
would  be  identical  to  that  predicted  for  a  0.01  lb/sec  purge  flow 
rate  in  the  char  motor  firings. 


The  anticipated  test  conditions  for  the  two  APG  tests  have 
been  presented  previously  in  Table  6-1  (page  6-4) .  As  explained 
earlier,  a  malfunction  of  the  recording  equipment  resulted  in  the 
requirement  that  both  tests  be  run  with  the  opposing  test  section 
being  carbon  cloth  phenolic.  The  arc  operating  characteristics 
for  these  test^  are  summarized  in  Table  6-2  (page  6-4) . 


6.3.2  Test  Results 

The  primary  test  data  from  the  narrow  view  angle  radiometer 
are  the  measured  emf  across  the  cons tantan  sensor  which  is  related 
to  total  hemispherical  radiative  heat  fluxes.  The  results  of  this 
measurement  are  shown  in  Figure  6-15  where  they  are  compared  with 
calculated  values  based  on  the  predicted  surface  temperature  of  the 
opposing  test  section.  The  maximum  variation  between  the  measured 
and  calculated  radiative  heat  flux  is  approximately  20  percent 
which  occurs  approximately  half  way  through  the  test.  However,  this 
variation  is  considerably  less  than  20  percent  during  the  majority 
of  the  test  and  is  negligible  during  the  initial  transient.  Also, 
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this  discrepancy  is  more  likely  due  to  an  inaccuracy  in  the  tempera¬ 
ture  indicated  by  the  pyrometer  than  an  inaccuracy  in  the  radiometer 
output . 

During  both  tests,  the  gas  purge  system  functioned  perfectly 
and  did  not  affect  the  operation  of  the  radiometer.  In  addition, 
visual  post-test  observations  showed  that  the  sapphire  lens  was 
clear  and  that  the  copper  housing  and  aperture  did  not  experience 
any  significant  heating.. 
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SECTION  7 
CHAR  MOTOR  TESTS 


The  total  and  radiative  heat  flux  measurement  systems  were  evalu¬ 
ated  through  two  nozzle  test  firings  which  were  performed  on  the  AFRPL 
40-inch  char  motor.  The  nozzle  designs  and  firing  conditions  were 
selected  to  test  these  systems  under  conditions  typical  of  their  appli¬ 
cation  in  rocket  nozzle  firings.  An  uncured  propellant,  APG  112D, 
which  has  a  flame  temperature  of  7010°R  at  1000  psia  was  used  in 
these  firings. 


The  first  nozzle  which  was  fired  on  August  13,  1969  contained  four 
total  heat  flux  ablative  plugs  and  one  radiation  sensor.  Three  of  the 
ablative  plugs  were  constructed  of  MX  4926  carbon  phenolic  which  formed 
the  throat  approach  and  the  nozzle  throat  region  while  the  other  was 
constructed  of  MX  2600-96  silica  phenolic  which  formed  the  nozzle  exit 
cone.  Apparent  failure  of  the  radiometer  purge  flow  supply  resulted 
in  only  minimal  radiometer  data.  The  total  and  effective  firing 
durations  were  54  and  48  seconds,  respectively,  with  the  chamber 
pressure  for  the  2.30  inch  diameter  nozzle  reaching  a  peak  of  1080 
psia  at  5  seconds  into  the  firing.  An  AFRPL  test  report  was  pub¬ 
lished  (Reference  7-1) . 


The  second  nozzle  which  was  fired  on  February  6,  1970  contained 
two  MX  4926  carbon  phenolic  total  heat  flux  measurement  systems  and 
two  radiometers.  One  radiometer  was  identical  to  that  used  in  the 
first  firing  while  the  other  was  identical  to  a  radiometer  scheduled 
for  a  firing  under  Contract  F04611-69-C-0065 .  A  molybdenum  transient 
slug  calorimeter,  designed  by  Atlantic  Research  Corporation  was  also 
included  in  the  nozzle  but  only  qualitative  information  was  obtained 
from  this  device.  All  other  measurement  systems  performed  as  ex¬ 
pected.  Total  and  effective  firing  durations  were  58  and  46  sec¬ 
onds,  respectively,  with  the  chamber  pressure  for  the  2.40  inch 
diameter  nozzle  reaching  a  maximum  of  930  psia  at  11  seconds  into 
the  firing. 


Detailed  descriptions  of  the  designs,  fabrication  technique  and 
firing  test  conditions  are  presented  in  Section  7.1.  The  overall 
ablative  and  structural  responses  of  the  two  nozzles  are  described  in 
Section  7.2.  The  response  of  the  total  heat  flux  measurement  systems 
is  discussed  in  Section  7.3.  Finally  the  measured  response  of  the 
radiative  heat  flux  measurement  systems  are  described  in  Section  7.4. 
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7-1  CHAR  MOTOR  NOZZLES 


Two  ablative  throat  char  motor  nozzles  were  designed  and  fabri¬ 
cated  under  this  contract.  These  nozzles  were  designed  to  incorpor¬ 
ate  the  heat  flux  measurement  systems  described  in  Section  4  at  a 
minimum  risk  to  the  nozzle  structural  integrity.  In  addition  the 
nozzle  components  were  designed  to  simplify  their  assembly  into  the 
steel  shell  in  order  to  minimize  the  risk  of  damaging  the  heat  flux 
measurement  systems.  The  nozzle  designs,  the  nozzle  component  fabri¬ 
cation,  and  the  assembly  of  the  nozzle  components,  including  the 
heat  flux  measurement  systems,  are  described  in  Section  7.1.1.  The 
environmental  conditions  achieved  in  each  of  the  char  motor  firings 
are  described  in  Section  7.1.2. 


7.1.1  Nozzle  Design,  Fabrication,  and  Assembly 

The  ablative  throat  nozzle  designs  for  the  two  40-inch  char 
motor  firings  are  shown  in  Figures  7-1  and  7-2.  The  internal  nozzle 
contour  for  both  of  these  nozzles  was  similar,  with  the  only  differ¬ 
ence  being  that  the  second  nozzle  had  a  throat  diameter  of  2.40 
inches  while  the  first  had  a  throat  diameter  of  2.30  inches.  The 
same  nozzle  structural  steel  shell  was  used  for  both  motor  firings, 
but  it  was  reworked  after  the  first  firing  to  accommodate  additional 
heat  flux  measurement  systems. 


For  both  nozzle  designs,  the  throat  approach  component  was  fabri¬ 
cated  from  MX  4926  carbon  cloth  phenolic  oriented  at  a  90°  layup  angle 
with  respect  to  the  nozzle  centerline.  This  material  was  chosen 
because  it  would  minimize  the  surface  recession  in  this  region  and, 
thus,  provide  a  relatively  smooth  entrance  to  the  nozzle  throat  region. 
The  90°  orientation  was  selected  in  order  to  minimize  fabrication  costs. 
This  orientation  did  result  in  the  laminates  being  parallel  to  the  sur¬ 
face  at  the  entrance  to  the  throat  approach.  However,  this  was  not 
considered  to  be  a  problem  because  of  the  relatively  low  convective 
heat  flux  in  this  region.  The  maximum  diameter  of  this  component  was 
dictated  by  the  inner  bore  diameter  of  the  flange  on  the  aft  closure 
of  the  char  motor.  The  throat  entrance  was  bonded  to  the  steel  shell 
and  the  nozzle  throat  assembly.  Epon  913  adhesive  was  used. 


The  nozzle  throat  ablative  component  was  also  fabricated  from 
MX  4326  carbon  cloth  phenolic  but  at  a  layup  angle  of  60°.  This 
material  was  chosen  because  it  was  a  well  characterized  material  and, 
thus,  met  the  requirement  for  the  total  heat  flux  measurement  system. 
The  layup  angle  was  chosen  to  coincide  with  that  typically  used  in 
nozzle  ablative  throats.  The  exit  cone  for  the  first  nozzle  was 
MX  2600-96  silica  cloth  phenolic  (A/A*  >.  3.0) (Figure  7-1).  This 
provided  an  evaluation  of  the  total  heat  flux  system  for  a  typical 
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FIGURE  7-1  NOZZLE  DESIGN  FOR  FIRST  CHAR  MOTOR  FIRING 


MOTOR  RXRXNO 


exit  cone  material.  For  simplicity,  carbon  phenolic  was  used 
throughout  for  the  second  nozzle  (Figure  7-2)  .  Also,  the  additional 
instrumentation  required  a  homogeneous  material  at  A/A*  =  3.0. 


The  nozzle  throat  diameter  was  chosen  so  that  the  maximum 
anticipated  chamber  pressure  history  was  740  psia.  The  predicted 
chamber  pressure  history  based  on  the  nozzle  throat  thermal  pre¬ 
dictions  described  in  Section  4  is  presented  in  Figure  7-3.  The 
highly  regressive  pressure  trace  is  due  to  the  high  throat  surface 
recession  associated  with  the  ablative  throat.  A  considerably 
higher  than  anticipated  maximum  chamber  pressure  was  achieved  in 
the  first  nozzle  firing  because  of  an-  increase  in  the  proptellant 
burning  surface  area  due  to  the  erosion  of  the  insulation  liner  in 
previous  motor  firings  and  because  of  a  higher  than  assumed  propel¬ 
lant  burning  rate.  As  a  result,  the  throat  diameter  for  the  second 
nozzle  was  increased  to  2.40  inches. 


The  backside  of  the  nozzle  throat  had  a  5-1/2°  taper  which 
insured  retention.  In  addition,  the  backside  of  the  throat  was 
bonded  to  the  steel  shell  with  Epon  913.  This  bondline  had  two 
purposes:  to  provide  a  seal  between  the  steel  shell  and  the  abla¬ 
tive  component  and  to  carry  the  shear  load  at  the  interface.  Note 
also  that  O-rings  were  used  to  seal  any  possible  gas  paths.  This 
redundancy  was  employed  because  of  the  difficulty  in  "assuring  that 
the  bond  would  be  a  100  percent  effective  gas  seal.  The  MX  2600-96 
silica  cloth  phenolic  exit  cone  of  the  first  nozzle  was  fabricated 
at  a  layup  angle  of  60°.  The  thickness  of  the  silica  phenolic  was 
defined  based  on  the  thermal  predictions  presented  in  Section  4. 
This  part  was  bonded  to  the  carbon  cloth  phenolic  throat  and  the 
steel  shell  using  Epon  913. 


For  the  first  nozzle  (Figure  7-1) ,  three  total  heat  flux 
sensors  were  located  in  the  nozzle  throat  region  corresponding  to 
area  ratios  of  -2.0,  1.0,  and  2.0.  In  addition,  a  narrow  view 
angle  radiometer  was  located  at  the  nozzle  throat  plane.  A  four'h 
total  flux  sensor  was  located  in  the  exit  cone  at  an  area  ratio  of 
4.0.  As  shown  in  Figure  7-2,  the  instrumentation  for  the  second 
nozzle  was  located  at  the  throat  and  at  an  area  ratio  of  3.0.  Each 
location  included  a  total  heat  flux  sensor  and  a  narrow  view  angle 
radiometer.  Wherever  possible,  the  measurements  were  separated 
radially  by  90°.  This  was  done  primarily  to  prevent  any  flow  per¬ 
turbations  generated  at  one  measurement  location  from  interfering 
with  the  measurement  system  at  another  location. 


With  the  exception  of  the  narrow  view  angle  radiometer  in 
the  exit  cone  of  the  second  nozzle,  all  the  measurement  systems  were 
bonded  to  the  ablative  component  in  which  they  were  installed.  This 
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was  also  true  of  the  molybdenum  slug  calorimeter  located  at  an  area 
ratio  of  3.0  in  the  second  nozzle.  The  narrow  view  angle  radiometer 
located  in  the  exit  cone  was  retained  by  threads,  and  an  O-ring  was 
used  to  provide  the  seal. 


As  shown  in  Figures  7-1  and  7-2,  slots  were  machined  in  the 
ablative  components.  In  most  locations,  these  slots  contained  either 
thermocouple  lead  wires  from  the  total  heat  flux  sensors  or  leads  and 
gas  purge  line  from  the  narrow  view  angle  radiometer.  (In  the  exit 
cone,  the  slot  was  required  only  in  the  assembly  of  the  nozzle.) 


The  structural  shell  was  fabricated  from  1022  steel.  A  conserva¬ 
tive  approach  was  taken  in  the  design  of  this  shell  due  to  the  numer¬ 
ous  cutouts  required  for  the  instrumentation.  The  size  and  location 
of  the  bolt  holes  for  mounting  the  nozzle  to  the  motor  case  were 
dictated  by  the  location  of  the  mating  holes  on  the  flange  of  the 
aft  closure.  As  shown  in  Figure  7-1,  the  steel  shell  had  a  flange 
at  the  exit  which  was  an  integral  part  of  the  shell.  However,  in 
order  to  remove  the  ablative  components  after  the  first  firing,  this 
flange  had  to  be  machined  off.  For  the  second  firing,  this  flange  was 
bolted  on  as  shown  in  Figure  7-2. 


The  fabrication  of  the  instrumentation  systems  was  described  in 
Section  4  and  was  performed  uy  Aerotherm  for  the  total  heat  flux 
sensors  and  by  Heat  Technology  Laboratory  for  the  narrow  view  angle 
radiometers.  Haveg-Reinhold  fabricated  the  steel  and  ablative  com¬ 
ponents.  After  the  ablative  components  were  fabricated,  they  were 
bonded  together  and  then  the  final  machining  was  performed.  This 
assembly  was  shipped  to  Aerotherm  for  the  installation  of  the  instru¬ 
mentation  systems.  Since  the  cylindrical  barrel  of  the  total  heat 
flux  sensors  was  intentionally  made  longer  than  the  mating  hole,  the 
sensors  were  machined  flush  with  the  surface  after  they  had  been  bonded 
in  place  and  after  the  protuberance  had  been  carefully  measured.  This 
measurement  was  required  in  order  to  accurately  define  the  distance 
of  the  in-depth  thermocouples  from  the  nozzle  surface.  After  assembly 
of  the  instrumentation,  the  ablative  component  was  returned  to  Haveg- 
Reinhold  where  it  was  bonded  into  the  steel  shell.  Aerotherm 
personnel  assisted  in  this  final  assembly  and  bonding  operation.  After 
curing  the  bond,  the  nozzle  assembly  was  shipped  to  AFRPL  where  it 
was  installed  on  the  40-inch  char  motor. 


The  necessary  thermocouple  lead  wire  and  gas  purge  line  hookups 
were  made  after  installing  the  nozzle  assembly  on  the  motor.  A  schematic 
of  the  gas  purge  system  is  shown  in  Figure  7-4  for  the  first  nozzle 


FIGURE  7-4  NITROGEN  G AS  PURGE ^  SYSTEM  FOR  THE  FIRST  CHAR  MOTOR  FIRING 


FIGURE  7-5  NITROGEN  CAS  I’URGF.  SYSTEM  FOR  THE  SECOND  CHAR  MOTOR  FIRING 


firing  and  in  Figure  7-5  for  the  second  nozzle  firing*  Both  sys¬ 
tems  contained  a  high  pressure  nitrogen  source,  and  the  flow  of 
the  nitrogen  gas  to  the  radiometers  was  controlled  by  one  or  more 
sonic  orifices  in  the  line.  The  flow  system  lines  for  the  first 
firing  consisted  of  high  pressure  flexible  tubing.  However,  be¬ 
cause  of  the  line  failure  caused  by  radiant  heating  from  the  high 
temperature  exhaust,  these  lines  were  changed  to  stainless  steel 
tubing  for  the  second  firing.  In  addition,  a  15  micron  filter  was 
installed  in  the  line  for  the  second  firing  to  insure  that  parti¬ 
cles  would  not  plug  the  sonic  orifices.  Also,  note  that  pressure 
transducers  were  located  both  upstream  and  downs tream  of  the  sonic 
orifices  for  the  second  firing. 


7.1.2  Char  Motor  Firing  Conditions 

The  propellant  description  and  firing  conditions  for  the 
two  char  motor  firings  are  presented  in  this  section  and  define 
the  environmental  conditions  to  which  the  heat  flux  measurement 
systems  were  exposed.  The  uncured  propellant,  APG  112D,  produced 
by  Atlantic  Research  Corporation  was  employed  in  both  test  firings 
The  elemental  composition  of  this  propellant  is  given  in  Table  7-1 
Under  the  assumptions  of  chemical  equilibrium,  the  thermodynamic 
state  of  this  elemental  composition  is  dependent  on  the  system 
pressure  and  temperature.  For  a  flame  temperature  of  7020°R  and 
a  chamber  pressure  of  1000  pisa,  the  molecular  composition  of  the 
combustion  products  is  given  in  Table  7-2. 


The  chamber  pressure  versus  time  traces  are  shown  in  Figure 
7-6.  Two  pressure  transducers  were  used  for  each  firing;  each 
agreed  within  3  percent  of  the  curves  shown  in  Figure  7-6  for  the 
entire  duration  of  the  respective  firings.  The  approximately  five 
seconds  prior  to  pressure  rise  in  the  second  firing  was  due  to  a 
delay  in  ignition.  However,  since  the  initial  throat  diameter  was 
slightly  larger  in  the  second  firing,  peak  pressure  was  lower  and, 
after  the  start-up  transients,  the  pressure  response  for  the  two 
firings  were  essentially  identical.  Pertinent  chamber  pressure 
information  from  the  two  char  motor  firings  are  compared  with  the 
anticipated  values  in  Table  7-3.  This  table  also  presents  propel¬ 
lant  ballistic  and  nozzle  throat  geometric  information. 


7.2  THERMAL  AND  ABLATION  NOZZLE  RESPONSE 

The  ablation  and  thermal  responses  of  the  MX  4926  carbon 
phenolic  and  MX  2600-96  silica  phenolic  nozzle  parts  to  the  con¬ 
ditions  presented  in  Section  7.1.2  are  described  in  this  section. 
Qualitative  descriptions  of  post-test  nozzle  appearances  and  over¬ 
all  response  trends  are  presented  and  then  measured  ablation  re¬ 
sults  are  reported.  In  general  the  response  of  the  nozzles  was 
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CHAR  MOTOR  FIRINGS 


consistent  with  the  pre-test  calculation  considering  that  the  peak 
chamber  pressures  were  higher  than  anticipated  due  to  the  problem 
discussed  in  Section  7.1.1. 


The  primary  purpose  of  the  firings  was  to  establish  convective 
and  radiative  heat  fluxes  sufficient  to  test  the  two  measurement  con¬ 
cepts.  This  was  certainly  satisfied.  The  surface  recession  and  in- 
depth  decomposition  of  the  MX  4926  carbon  phenolic  and  the  MX  2600-96 
silica  phenolic  components  represented  the  type  of  ablative  response 
which  the  two  measurement  systems  would  encounter  in  typical  solid 
rocket  motor  applications.  Moreover,  since  the  char  motor  chamber 
conditions  achieved  during  each  firing  were  quite  similar,  a  valuable 
indication  of  the  repeatability  of  the  measurements  was  provided. 


The  presence  of  the  numerous  sensors  embedded  at  various  loca¬ 
tions  in  the  ablative  nozzles  appeared  to  have  only  trivial  effects 
on  the  overall  nozzle  ablation  response.  The  most  serious  deviation 
occurred  downstream  of  the  holes  provided  for  the  narrow  view  radiom¬ 
eters.  At  these  locations  the  flow  contour  discontinuity  resulted  in 
a  local  gouge  in  the  ablative  material.  This  effect  was  most  signifi¬ 
cant  at  the  throat  of  the  first  nozzle  for  which  the  radiometer  nitro¬ 
gen  purge  gas  supply  system  had  failed.  Other  small  local  asymmetries 
occurred  in  the  two-nozzles,  but  these  did  not  appear  to  be  related  in 
any  way  to  the  locations  of  the  total  or  radiative  heat  flux  measure¬ 
ment  systems.  It  was  apparent  that  the  presence  of  the  imbedded  in¬ 
strumentations  did  not  impair  the  performance  of  the  rocket  nozzle. 


The  net  surface  recession  at  the  throats  of  the  two  nozzles 
was  slightly  higher  than  anticipated,  consistent  with  the  higher 
chamber  pressures.  The  throat  surface  recession  and  char  depth  and 
their  asymmetry  are  shown  in  Figures  7-7a  and  7-7b  for  the  respective 
firings.  The  figures  show  measured  post-test  surface  and  char  depth 
radii  at  seven  angular  locations.  Also  shown  are  the  pre-test  diameter 
and  post-test  average  surface  and  char  depth  diameters.  These  figures 
show  that  the  maximum  asymmetry  in  surface  recession  is  generally  less 
than  10  percent  of  the  respective  average  diameter. 


Measurements  of  recession  and  char  depth  at  various  nozzle 
locations  are  shown  in  Figures  7-8  and  7-9.  Of  particular  note  in 
these  figures  is  the  discontinuous  change  in  recession  which 
occurred  about  three  inches  upstream  of  the  throat  in  both  nozzle 
firings.  The  abrupt  onset  of  ablation  is  apparently  due  to  a  flow 
field  anomaly  and  not  a  material  problem.  Candidate  explanations 
include  transition  from  laminar  to  turbulent  flow,  a  separation/re¬ 
attachment  shock,  and  an  unusual  sonic  region  shock  pattern.  The 
actual  explanation  has  not  been  defined  and  was  beyond  the  scope  of 
this  program.  A  summary  of  the  ablation  response  shown  in  Figures 
7-8  and  7-9  is  given  in  Table  7-4. 
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7.3  TOTAL  HEAT  FLUX  MEASUREMENT  SYSTEM 

The  operation  of  the  six  total  heat  flux  measurement  systems 
tested  in  the  two  rocket  nozzle  firings  is  described  in  this  section. 
A  brief  summary  of  the  installation  and  fabrication  details  of  the 
sensors  which  are  significant  to  the  data  reduction  procedure  is 
discussed  in  Section  7.3.1.  The  response  of  each  sensor  is  then  pre¬ 
sented  in  Section  7.3.2,  and  the  results  of  the  data  reduction  are 
reported  in  Section  7.3.3. 


7.3.1  System  Design  and  Installation 

The  basic  design  of  the  total  heat  flux  measurement  system  has 
been  described  in  Section  4.  Specific  details  about  the  particular" 
sensors  utilized  in  the  two  char  motor  firings  of  this  program  are 
reported  in  Table  7-5. 


For  the  first  nozzle  firing  the  various  sensor  locations  were 
selected  so  that  the  total  heat  flux  measurement  systems  would  be 
exposed  to  a  wide  range  of  convective  and  radiative  boundary  con¬ 
ditions.  In  addition  the  effects  of  material  type  and  reinforce¬ 
ment  tape  layup  angle  on  measurement  system  operation  could  be  examined 
since  the  layup  angle  (referenced  to  the  nozzle  contour  tangent) 
varied  with  area  ratio.  The  orientation  of  each  sensor  was  shown  as 
a  part  of  the  overall  nozzle  response  diagrams  (Figures  7-8  and  7-9)  . 

As  was  mentioned  previously,  the  angular  location  of  each  sensor  was 
selected  to  minimize  the  possibility  of  disturbances  in  the  local  flow 
resulting  from  upstream  sensors.  Based  on  the  post-fire  appearance  of 
the  nozzles,  this  precautionary  measure  was  probably  not  necessary 
since  only  small  local  disturbances  could  be  attributed  to  the 
sensors . 


The  total  heat  flux  measurement  plugs  utilized  in  the  second 
char  motor  firing  were  positioned  to  confirm  results  of  the  first 
firing  and  obtain  data  for  an  additional  nozzle  location.  Since 
results  of  the  first  firing  were  satisfactory,  the  basic  sensor 
design  was  not  altered  for  the  second  firing.  The  slight  changes  in 
thermocouple  depths  shown  in  Table  7-5  between  the  nozzle  firings 
(and  between  each  sensor  plug)  resulted  from  machining  tolerances 
allowed  during  sensor  fabrication. 


The  fabrication  and  installation  techniques  employed  to  position 
each  sensor  plug  in  the  two  ablative  nozzles  were  as  described  in 
Sections  4.1,  4.2,  and  7.1.  In  the  assembly  of  the  steel  shell  and 
ablative  liner,  special  care  was  taken  to  insure  that  the  lead  wires 
remained  in  their  respective  slots  in  the  ablative  -liner  and  were 
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properly  aligned  and  fed  through  the  holes  in  the  steel  shell  as  it 
was  slipped  on  to  the  liner.  Even  though  the  installation  process 
was  successful  and  the  same  process  was  used  for  both  nozzles,  some 
thermocouple  circuits  were  lost  during  this  installation  for  the 
second  nozzle.  These  included  three  thermocouple  circuits  going  to 
the  total  heat  flux  sensor  at  an  A/A*  =  3.0  and  the  circuit  to  the 
Atlantic  Research  Corporation  (ARC)  calorimeter.  These  failures  were 
detected  by  continuity  checks  of  the  circuits.  The  continuity  checks 
were  made  continually  through  the  fabrication  and  installation  proc¬ 
ess.  If  these  checks  had  revealed  a  significant  number  of  circuit 
failures,  the  sensor  would  have  been  replaced.  However,  any  repair 
of  the  sensor  was  impossible  after  assembly  of  the  ablative  part  in 
the  steel  shell. 


7.3.2  Sensor  Response  - 

The  response  of  the  total  heat  flux  sensors  was  consistent 
with  the  requirements  for  effective  measurement  of  the  heat  flux 
boundary  conditions.  One  of  the  primary  uncertainties  in  the 
expected  response  of  the  sensors  was  the  degree  to  which  flow  dis¬ 
continuities  would  be  generated  during  the  ablation  of  the  heat  flux 
sensors.  The  design  was  directed  towards  minimizing  such  disturb¬ 
ances,  and  the  post-firing  appearances  of  the  sensor  locations  indi¬ 
cated  that  this  design  requirement  was  satisfied.  Photographs  of 
five  of  the  six  plug  locations  are  given  in  Figure  7-10.  The  other 
location  was  not  photographed  because  the  photographer  could  not  find 
the  sensor /nozzle  interface.  The  most  significant  deviation  from 
ideal  performance  occurred  at  the  throat  location  of  the  first  nozzle 
firing  where  slight  gouging  occurred  at  the  thermocouple  lead  wire 
slots  on  the  side  of  the  sensors  (Figure  7-10) .  This  gouging  caused 
premature  failure  of  the  thermocouples  as  ablation  occurred.  None 
of  the  other  sensors  including  the  one  at  the  throat  of  the  second 
nozzle  exhibited  such  significant  gouging. 


As  discussed  earlier  in  the  descriptions  of  the  measurement 
system  concept,  all  information  about  surface  heat  flux  boundary 
conditions  are  derived  from  the  subsurface  thermocouple  measurements. 
Moreover,  since  the  evaluation  of  the  surface  events  from  these 
measurements  are  extremely  sensitive  to  experimental  error,  the  data 
recording  system  and  the  thermocouple  calibrations  must  be  accurate. 

In  the  test  firings,  the  AFRPL  digital  data  recording  system  recorded 
thermocouple  millivolt  output  continuously  and  it  was  sampled 
approximately  every  10  milliseconds  in  the  data  reduction  process. 
Standard  NBS  calibration  curves  were  utilized  for  the  Chromel-Alumel 
thermocouples.  This  standard  was  checked  against  manufacturers  speci¬ 
fications  for  the  particular  wire  lot  used  for  the  total  heat  flux 
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FIGURE  7-10  PHOTOGRAPHS  OF  NOZZLE  SURFACES  SHOWING 
ABLATIVE  PLUG  SENSORS 


-ip* 
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sensor  and  was  found  to  be  accurate.  For  the  tungsten-5  percent 
rhenium/tungsten-26  percent  rhenium  thermocouples  (W5Be/W26Re) t  no 
NBS  standard  is  available.  Manufacturer's  specifications,  available 
to  4200°F,  were  curve-fit  and  utilized  to  reduce  the  data.  A  differ¬ 
ent  wire  lot  was  used  for  each  char  motoj:  nozzle;  hence  two  calibra¬ 
tion  curves  were  required.  For*  temperatures  from  4200°F  to  5400°F, 
data,  of  Reference  4-12  were  curve-fit.  The  calibration  equations 
utilized  to  reduce  the  measured  thermocouple  millivolt  outputs  are 
summarized  in  Table  7-6. 


The  reduced  thermocouple  data  are  presented  in  Figure  7-11.  Also 
shown  are  thermal  response  predictions  to  be  described  in  Section 
7.3.3.  Some  anomalies  are  apparent  from  Figure  7-11.  The  most  sig¬ 
nificant  .is  the  shift  in  the  rate  of  temperature  increase  which 
occurred  for  practically  all  thermocouples  at  3000°R.  An  error  in  a 
calibration  equation  was  suspected  so  the  millivolt  output  was 
plotted  for  the  A/A*  =  2 JO  sensor  of  the  first  firing.  This  is  pre¬ 
sented  as  Figure  7-12  and  shows. the  same  trend.  Therefore,  either 
the.  trend  is  real  or  some  unexpected  calibration  shift  or  electrical 
grounding  of  the  circuit  occurred  to  compromise  the  accuracy  of 
measurements  above  3000OR.  Despite  these  anomalies,  the  accuracy  of 
the  heat  fluxes  evaluated  from  the  overall  system  response  was  not 
significantly  impaired  because  sufficient  data  are  available  at 
temperatures  below  3000°R. 


The  anomalous  responses  of  the  three  deepest  thermocouples  of 
the  throat  sensor  in  the  first  nozzle  (Figure  7-11)  are  believed  to 
be  a  result  of  the  gouging  in  the  thermocouple  lead  slots  described 
above.  Sufficient  data  were,  however,  obtained  in  this  case  to 
provide  a  basis  for  extrapolating  heat,  flux  to  the  end  of  the  firing. 
In  conclusion,  although  a  few  anomalies  in  the  thermocouple  data  did 
occur,  enough  redundancy  was  provided  in  the  measurement  system  con¬ 
cept  and  design  to  allow  the  partial  definition  of  heat  flux  histor¬ 
ies  for  all  of  the  locations  where  a  sensor  was  placed. 


7.3.3  Total  Heat  Flux  System  Evaluation 

In  the  previous  section  the  response  of  the  total  heat  flux 
sensors  was  described  and  the  temperature  history  measurements  were 
presented.  The  evaluation  of  the  total  heat  fluxes  from  the  firing 
condition  information  and  sensor  output  is  presented  in  this  section. 
The  basic  data  reduction  technique  was  described  previously  in 
Section  5.  Specific  details  about  the  data  reduction  for  the  ablative 
materials  tested  in  the  two  char  motor  firings  are  presented  in  the 
following  paragraphs.  In  addition  the  analysis  is  concluded  with  a 
comparison  of  predicted  heat  flux  to.  che  measured  values  and  a 
discussion  of  measurement  system  performance. 
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where  U  =  (X(mv)  -  XO)/XI 
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FIGURE  7-12  TYPICAL  THERMOCOUPLE  MILLIVOLT  OUTPUTS  SHOWING 
TEMPERATURE  RISE  RATE  DECREASE  NEAR  3000°R 


The  steps  in  the  evaluation  of  total  incident  heat  flux  may  be 
summarized  as  follows: 

•  Analysis  of  overall  consistency  between  the  tempera¬ 
ture  histories,  thermocouple  depths,  ablative  material 
thermophysical  properties,  and  the  one-dimensional  thermal 
response  model. 

•  Estimation  of  first  guess  surface  temperature  and 
recession  histories. 

•  Iterative  evaluation  of  surface  temperature,  surface 
recession  and  in-depth  decomposition  events. 

•  Consideration  of  char  swell. 

\  i 

•  Evaluation  of  surface  thermochemical  events  from 
parametric  surface  mass  balance  calculations. 

•  Evaluation  of  surface  total  heat  flux. 

Discussions  of  the  techniques  employed  to  accomplish  each  of  these 
steps  are  given  in  Sections  7. 3. 3.1  through  7. 3. 3. 6.  Concluding  remarks 
are  then  given  in  Section  7. 3. 3. 7. 


7. 3. 3.1  Verification  of  Heat  Conduction  Model  and  Temperature 

Response  Data. 

One  advantage  of  calculating  heat  flux  from  multiple  thermo¬ 
couple  temperature  histories  is  that  redundant  information  is  avail¬ 
able.  For  an  ideal  one-dimensional  slab  of  ablative  material  being 
heated  on  one  surface  and  insulated  on  all  others,  the  surface  tem¬ 
perature  and  heat  flux  conducted  into  the  material  can  be  defined 
from  a  single  in-depth  temperature  history,  provided  the  pertinent 
material  thermal  and  chemical  properties  are  accurately  known.  A 
well  characterized  ablative  material  used  as  the  exposed  surface  ma¬ 
terial  in  a  rocket  nozzle  is  a  good  approximation  of  this  ideal  model. 
The  measurement  of  temperatures  at  more  than  one  in-depth  location 
provides  sufficient  information  to  estimate  how  closely  the  measured 
response  approaches  that  predicted  by  the  ideal  model.  This  compari¬ 
son,  if  satisfactory,  between  the  predicted  and  measured  in-depth 
material  thermal  response,  using  the  material  thermal  and  chemical 
properties  for  MX  2600-96  and  MX  4926  described  in  Section  4,  pro¬ 
vides  the  verification  of  the  theoretical  model. 


For  the  temperature  history  data  obtained  in  this  program,  these 
verifications  of  data  consistency  were  performed  by  inputting 
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thermocouple  locations  as  pseudo-surfaces  and  by  inputting  thermocouple 
data  as  pseudo-surface  temperature  into  the  Aerotherm  Charring  Material 
Ablation  (CMA)  computer  code.  The  code  then  calculated  the  transient 
thermal  response  of  the  material  below  the  pseudo-surface.  The 
degree  to  which  the  measured  in-depth  temperatures  could  be  predicted 
was  an  indication  of  the  quality  of  the  data.  A  typical  example  of 
the  results  of  one  of  these  calulations  is  given  in  Figure  7-13. 

Notice  that  since  the  thermocouple  closest  to  the  surface  was  lost 
during  the  firing,  a  shift  in  the  driver  was  required.  This  shift 
was  accomplished  by  impressing  a  large  recession  rate  for  a  short 
time  so  that  the  surface  moved  to  the  location  of  the  third  thermo¬ 
couple.  At  the  same  time  the  pseudo-surface  temperature  was  de¬ 
creased  to  the  measured  value  of  the  third  thermocouple.  This 
rapid  shift  caused  the  anomaly  seen  in  the  prediction  for  the  fourth 
thermocouple  at  20  seconds.  Despite  anomalies  of  this  type,  these 
calculations  showed  that  in  general  the  thermocouple  depths  and  cor¬ 
responding  temperature  data  obtained  from  the  sensors  were  consis¬ 
tent  with  the  thermal  properties  and  heat  conduction  model  utilized 
in  the  data  reduction  procedure. 


7. 3.3.2  Estimation  of  Surface  Response  Boundary  Conditions 

The  second  step  in  the  evaluation  of  the  total  heat  flux  was 
to  estimate  the  surface  temperature  and  recession  responses  of  each 
ablative  plug.  This  information  is  the  input  boundary  conditions  for 
the  computation  of  the  overall  material  thermal  response.  The  reces¬ 
sion  response  as  a  function  of  firing  time  is  defined  quite'  well  from 
pre-  and  post-test  diameter  measurements  and  from  the  times  at  which 
thermocouples  become  exposed  at  the  surface.  Plots  of  recession .versus 
time  for  each  total  heat  flux  sensor  plug'  are  shown  in  Figure  7-14. 
Since  the  thermocouple  typically  failed  before  reaching  the  surface, 
these  points  fall  to  the  left  (or  above)  the  actual  recession  his¬ 
tory.  The-  initial  estimate  of  the  transient  surface  recession  was 
based  on  this  interpretation  of  thermocouple  failure  and  on  the 
post-test  measurement.  Figure  7-14b  also  shows  the  results  of  pre¬ 
dictions  based  on  equilibrium  surface  thermochemistry  which  were 
performed  as  a  part  of  a  heat  flux  measurement  evaluation  analysis. 


Estimation  of  surface  temperature  histories  was  somewhat  more 
difficult,  especially  for  those  plugs  where  surface  recession  rate 
was  large.  The  technique  used  was  to  cross  plot  the  transient  re¬ 
sponses  against  thermocouple  location  at  several  specified  times 
and  to  extrapolate  the  resulting  temperature  profiles  to  the  surface 
locations  defined  as  described  above.  Due  to  the  relatively  high 
surface  recession  rate  and  the  relatively  high  surface  heat  flux  ex¬ 
perienced  by  several  of  the  plugs  in  these  nozzle  test  firings,  the 
profiles  were  so  steep  that  only  two  or  three  of  the  thermocouples 
could  be  used  in  establishing  this  profile.  Moreover,  the  anomalies 
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DATA:  DEPTH  BELOW  ORIGINAL 

THERMOCOUPLE  SURFACE  (INCH) 


FIGURE  7-13  ASSIGNMENT  OF  THERMOCOUPLE  MEASUREMENTS  TO  VERIFY  THE 

TEMPERATURE  DATA  AND  THE  IN-DEPTH  HEAT  CONDUCTION  MODEL 


RECESSION,  INCHES 


observed  above  3000°R  certainly  compromise  this  technique.  An  exam¬ 
ple  of  the  application  of  this  technique  is  shown  in  Figure  7-15. 

The  definition  of  the  shapes  of  the  temperature  profiles  are  assisted 
by  the  profiles  computed  during  the  data  verification  procedure.  These 
calculations,  along  with  the  data  points  available,  provided  sufficient 
information  to  estimate  the  surface  temperature  histories  for  the  ini¬ 
tial  input  to  the  iterative  thermal  response  evaluation  procedure. 


A  tradeoff  exists  between  the  accuracy  of  the  first  guess  sur  - 
face  temperature  and  the  first  guess  recession.  If  recession  is  high, 
many  thermocouples  wi,ll  be  broken  providing  numerous  data  points  for 
recession  history  estimation.  If  recession  is  low,  the  thermocouples 
will  provide  numerous  temperature  profile  data  points,  and  the  accu¬ 
racy  of  the  estimated  surface  temperature  will  be  improved.  Finally, 
since  these  are  only  first  guesses  for  the  subsequent  iterative  eval¬ 
uation  procedure,  uncertainties  in  these  estimates  do  not  signifi¬ 
cantly  affect  the  accuracy  of  the  final  heat  flux  measurement.  They 
only  increase  the  difficulty  of  converging  rapidly  to  the  answer. 

7. 3. 3. 3  Iterative  Evaluation  of  Surface  and  In-Depth  Responses 

Once  the  first  guesses  for  the  ablative  plug  surface  bound¬ 
ary  conditions  have  been  defined,  iterative  solutions  and  adjustments 
of  the  input  values  of  surface  recession  and  surface  temperature  are 
performed.  For  these  calculations  t^e  CMA  computer  code  is 
utilized  in  the  mode  where  the  heated  surface  response  (surface  re¬ 
cession  rate  and  temperature)  are  input  directly.  The  code  computes 
the  overall  heat  conduction  and  thermal  decomposition  response  of 
the  material,  including  predictions  of  thermocouple. measurements. 
Iteration  on  surface  boundary  condition  input  is  continued  uptil  the 
overall  transient  response  of  the  material  is  predicted  as,  closely 
as  desired. 


For  each  of  the  total  heat  flux  measurement  systems  tested  in 
the  two  char  motor  firings,  an  average  of  four  iterations  was  per¬ 
formed.  Figure  7-llf  (page  7-32)  presented  typical  iterative  values 
of  surface  temperature  and  compared  the  computed  in-depth  response 
to  thermocouple  data.  For  this  sensor  only  two  tuermocouple  cir¬ 
cuits  remained  after  nozzle  fabrication,  but  the  resulting  surface 
temperature  was  still  well  defined.  More  thermocouples,  as  were  y 
available  for  all  of  the  other  sensors,  decrease  the  importance  ofr 
errors  in  thermocouple  depths  or  in  thermocouple  data  readings  be¬ 
cause  such  anomalies  would  be  apparent  or  averaged  out.  The  final 
surface  temperature  histories  and  the  corresponding  in-depth  thermo¬ 
couple  response  predictions  £or  all  sensors  are  compared  to  the 
measured  thermocouple  data  in  Figures  7-lla  through  7-llf. 

In  addition  to  defining  the  surface  response  of  the  ablative 
plug,  the  iterative  evaluation  of  in-depth  response  defines  several 
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of  the  ablation  and  heat  transfer  results  required  in  the  final  sum¬ 
mation  of  surface  heat  fluxes.  These  quantities  have  been  described 
in  Section  5.  Some  of  the  representative  values  of  the  variables 
for  t,he  six  ablative  sensors  tested  in  this  program  are  given  in 
Table  7-7.  For  each  plug,  data  are  given  four  times  during  the  re-  \ 
spective  firings. 


7. 3. 3. 4  Char  Swell  Considerations 

The  energy  flux  terms  associated  with  chemical  reactions  at 
an  ablating  surface  are  strongly  dependent  on  the  rate  of  surface 
mass  removal  and  hence  on  the  surface  recession  rate.  The  surface 
recession  history  determined  by  the  technique  described  in  the  pre¬ 
vious  sections  accounts  for  two  surface  movement  effects: 

•  Surface  regression  due  to  chemical  reactions 

•  Surface  advancement  due  to  char  swelling 

These  response  mechanisms  are  illustrated  in  Figure  7-16.  The  actual 
surface  movement  defined  in  Section  7. 3. 3. 3  is  less  than  that  due  to 
thermochemical  attack  by  the  net ^amount  of  char  swell.  That  is, 

S  =  St  -S' 
m  tc  swell 

where 

S  =  measured  surface  recession  (Section  7. 3.3.3) 
m 

Sfcc  f  net^  thermochemical  regression 

S  , ,  =  net  surface  movement  due  to  swell 
swell 

t. 

This  swelling  has  been  found  to  occur  in  MX  49?6  carbon  phenolic 
but  not  in  MX  2600-96  silica  phenolic.  Th^  causes  of  the  phenomenon 
are  discussed  in  Reference  4-6.  In  that  study  it  was  found  that  the 
amount  of  char  swell  was  proportional  to  the  char  thickness  an^  that 
the  proportionality  constant  for  the  carbon  phenolic  is  0.135.  That 
is. 


(7-1) 


Sswell 


0.135  t 

c 


(7-2a) 


where  :  is  the  char  thickness.  Based  on  the  reasonable  assumption 
that  the  proportionality ; factor  is  constant  with  time.  Equation 
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(7-2a)  says  that  the  rate  of  surface  movement  due  to  char  swell  is 
proportional  to  the  rate  of  char  thickness  increase.  That  is 


S  0.135 

swell 


(7-2b) 


where  the  dot  superscripts  indicate  derivatives  with  respect  to  time 
(fates) . 


To  account  for  this  swelling  in  the  evaluation  of  chemically 
related  energy  flux  terms,  a  mass  rate  term  has  been  defined  as 
follows 


m  =  p  S  , , 

c  ,,  Kc  swell 
swell 


(7-3) 


where  p  is  the  char  density  before  swelling.  This  mass  rate  has 
been  adaed  to  the  mass  rate  defined  by  the  technique  of  Section 
7. 3. 3. 3  to  obtain  the  total  effective  thermochemical  mass  ablation 
rate.  Values  obtained  from  Equation  (7-3)  are  reported  in  Table 
7-7.  In  addition,  the  heat  flux  associated  with  this  ablation  rate 
is  given. 


7. 3. 3. 5  Surface  Thermochemical  Response  Evaluation 

The  procedure  described  in  Section  5  was  followed  for  the 
evaluation  of  the  quantities  associated  with  the  surface  chemical 
events.  The  assumption  that  the  surface  erosion  of  the  plug  is 
controlled  by  heterogeneous  equilibrium  chemistry  which  is  incorpo¬ 
rated  in  this  procedure  has  been  substantiated  by  the  results  pre¬ 
sented  in  Section  7. 2. 3. 2.  Moreover,  since  mechanical  removal  of 
surface  char  did  not  appear  to  be  a  significant  recession  mechanism 
in  either  the  MX  4926  carbon  phenolic  plugs  or  the  MX  2600-96  silica 
phenolic  plug,  the  required  evaluation  was  straightforward.  This 
fact  was  somewhat  surprising  because  usually  the  silica  in  the  MX 
2600-96  forms  a  liquid  layer.  In  this  case,  however,  the  surface 
temperature  was  apparently  high  enough  that  an  SiC>2  melt  did  not 
form  but  rather  the  silica  decomposed  into  gaseous  SiO  and  O2.  This 
was  confirmed  by  the  higher  than  expected  recession  rate. 


As  explained  in  Section  5,  the  independent  variables  for  defin¬ 
ing  chemically  related  heat  flux  are  obtained  from  the  in-depth  ther¬ 
mal  response  evaluation  and  from  the  rocket  motor  firing  conditions. 
These  variables  are: 


/ 
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ratio  of  B1  to  B* 
g  c 

B7B'  -=  m  /m 
gr  c  g  c 

surface  temperature 

T 

w 

Local  static  pressure 

P 

e 

propellant  and  material  compositions 

Typical  plots  of  the  several  dependent  variables  needed  for  the  total 
heat  flux  evaluation  are  presented  in  Figures  7-17  through  7-20  as  a 
function  of  the  independent  variables.  A  summary  of  the  figures  and 
the  variables  obtained  from  the  figures  is  given  below. 


Figure 

Number 

Independent 

Variable 

Interpolation 

Variable 

Dependent 
Variable (s) 

7-17 
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T 

B'  ,p  U  Cu 

e  g  e  g 
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7-18 
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w 

9 

w 

7-19 

T 

B’ 

T 

Z*  h.  w 

w 

9 

Xw  1 

7-20 

T 

T 

K.  h.  w 

7-20 

w 

T 

Xe  1 

T 

v*  u  w 

K.  n  . 

w 

i  i 
e 

Surface  thermochemistry  maps  of  the  type  shown  in  Figures  7-17 
through  7-20  were  generated  for  at  least  two  static  pressures  at 
each  ablative  plug  location.  The  results  for  intermediate  pressure 
were  obtained  by  interpolation.  The  results  for  four  times  during 
the  respective  firings  are  presented  in  Table  7-8  for  each  ablative 
plug.  Some  other  properties  of  the  propellant/surface  interaction 
which  were  required  to  obtain  the  values  of  q  .  .  and  q.  .  are  also 
shown  in  Table  7-8.  abl 
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Enthalpy  Potential  -  fh^EI^  hj  I  +  ^  x  ^rz^  w  -  IZi  hx  w^;  C^/C,,  -  0.70 


I 


7. 3. 3. 6  Evaluation  of  Total  Heat  Flux 

Values  of  the  various  component  heat  flux  terms  defined  in 
Section  5  are  given  in  Tables  7-7  and  7-8.  The  total  incident  heat 
fluxes  were  evaluated  by  the  surface  energy  balance  relation  given 
as  Equation  (5-1) .  Although  Tables  7-7  and  7-8  only  show  represen¬ 
tative  flux  values  for  four  times  during  the  nozzle  firing,  more 
points  were  calculated  to  provide  continuous  knowledge  of  the  varia¬ 
tions  in  total  flux.  Figure  7-21  shows  these  curves. 


7. 3.3.7  Discussion  of  Total  Heat  Flux  Measurements 

Several  trends  in  the  total  heat  flux  measurements  confirm 
the  validity  and  accuracy  of  the  measurement  concept.  The  total 
incident  heat  flux  histories  for  the  various  locations  in  the  two 
char  motor  nozzles  have  been  presented  as  Figure  7-21.  The  clearest 
and  probably  the  most  indicative  trend  is  apparent  in  the  responses 
of  the  MX  4926  carbon  phenolic  sensors  shown  in  Figure  7-2 la.  All 
of  these  sensors  recorded  a  decreasing  .total  incident  heat  flux  with 
firing  time.  This  trend  was  expected  since  the  heat  transfer  coef¬ 
ficient  for  turbulent  flow  is  roughly  proportional  to  the  0.8  power 
of  the  chamber  pressure.  Since  a  significant  radiative  heat  flux 
component  is  part  of  the  total  incident  flux,  the  overall  drop  in 
the  total  was  not  in  general  proportional  to  the  0.8  power  pressure. 


Another  significant  trend  indicated  in  Figure  7-21a  is  the  vari¬ 
ation  of  total  heat  flux  with  area  ratio.  For  the  subsonic  location 
(A/A*  =  -2.0)  and  the  throat  location,  the  total  heat  flux  histories 
are  definitely  higher  than  those  obtained  at  downstream  supersonic 
locations.  The  inlet  section  is  subjected  to  higher  radiation  heat¬ 
ing  and  recovery  enthalpies  and  the  throat  location  is  subjected  to 
larger  heat  transfer  coefficients  than  downstream  locations.  Hence 
this  trend  certainly  seems  reasonable. 


Comparisons  between  results  for  the  two  char  motor  firings  pro¬ 
vide  additional  verification.  At  the  throat  location  where  a  direct 
comparison  was  possible  because  firing  conditions  were  essentially 
identical,  the  difference  in  total  heat  flux  is  less  than  18  percent 
(i.e.  within  the  expected  accuracy  of  the  system).  The  total  heat 
flux  measurements  are  repeatable  and  the  quoted  accuracies  are 
reasonable . 


The  total  heat  flux  history  shown  for  the  A/A*  =  3.0  location 
of  the  second  char  motor  nozzle  is  anomalous.  For  this  sensor,  only 
two  of  the  five  original  thermocouple  circuits  remained  in  tact 
after  nozzle  fabrication,  and  therefore  the  results  are  compromised. 
Figure  7-llf  shows  that  the  most  useful  data  from  this  sensor  was 
taken  between  12  and  27  seconds  of  the  firing.  In  this  range  the 


\ 
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TOTAL  INCIDENT  HEAT  FLUX.  ATU/FTZ-SEC 


a)  MX4926  CARBON  PHENOLIC  TOTAL  HEAT  FLUX  SENSORS 


FIGURE  7-21  TOTAL  HEAT  FLUX  DETERMINED  FROM  MEASUREMENT  SYSTEM  DATA 


b)  HX2600-96  SILICA  1-IH.N'Jl.IC  TOTAL  III  AT  I  IJJX  :.l  MlHIH  A/A*  =  4.0 


hc.uiii:  7-21  coticumni 
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reported  heat  flux  values  agree  quite  well  in  magnitude  with  those 
taken  at  the  A/A*  =  2.0  location  of  the  first  nozzle.  It  would  be 
expected  that  these  two  locations  have  approximately  the  same  read¬ 
ings  since  no  drastic  changes  in  boundary  conditions  occur  between 

A/A*  =  2.0  and  3.0. 

}  .  ] 

\ 

Similarly,  the  total  heat  flux  determined  from  the  response  of 
the  MX  2600-96  silica  phenolic  plug  is  less  accurate  for  firing 
times  greater  than  20  to  25  seconds.  This  plug  location  experienced 
a  much  larger  recession  rate  than  expected  and,  therefore,  the  thermo¬ 
couples  became  exposed  apd.  were  lost  early  in  the  firing.  The  total 
heat  flux  history  presented  in  Figure  7-21b  shows  the  results  of 
extrapolated  data  (past  25  seconds)  as  dashed  lines.  The  total  in¬ 
cident  heat  flux  values  determined  from  the  silica  phenolic  plug  are 
also  potentially  in  error  because  the  equilibrium  thermochemical  re¬ 
sponse  of  the  silica  phenolic  is  not  well  understood  for  the  high 
surface  temperatures  obtained  in  this  case.  No  conclusions  about 
the  relative  magnitude  of  the  silica  phenolic  and  carbon  phenolic  * 
flux  results  are  possible  since  the  total  incident  heat  flux  is  a 
strong  function  of  the  significantly  different  surface  chemical 
events .  « 


In  an  attempt  to  further  verify  the  accuracy  of  the  total  heat 
flux  values  reported  for  the  MX  4926  carbon  pheinolic  sensors,  abla¬ 
tion  and  thermal  response  predictions  (using  the  ACE  and  CMA  com¬ 
puter  codes)  for  the  plug  locations  in  the  second  char  motor  firing 
were  performed.  (The  resultant  surface  recession  predictons  are  in¬ 
cluded  in  Figure  7-14b.)  For  the  throat  location,  the  predicted 
values  of  total  incident  heat  flux  and  incident  convective  heat  flux 
are  compared  to  measured  total  flux  in  Figure  7-22.  This  prediction 
required  an  assumption  about  the  incident  radiation  heit  flux.  Based 
on  the  radiation  measurements  described  in  the  following  section,  it 
was  assumed  that  the  incident  radiation  was  the  surface  absorptance 
time ,  the  average  of  the  black  body  radiation  from  the  free  stream 
temperature  and  the  blao^  body  radiation  from  the  opposite  wall. 

Thio  came  to  500  Btu/ft^sec  and  yas  assumed  constant  throughout  the 
firing.  Figure  7-22  shows  that  the  determined  and  predicted  heat 
flux  values  are  comparable.  In  conclusion,  it  is  apparent  that  the 
total  incident  heat  (lux  measurement  systems  provide  data  which  are 
self-consistent  and  accurate  to  within  ±20  percent. 


7.4  RADIATIVE  HEAT  FLUX  SYSTEM 

The  performance  of  each  of  the  three  narrow  view  angle  radiom¬ 
eters  in  the  two  char  motor  firings  is  described  tin  this  section. 

In  Section  7.4.1,  the  details  of  the  design  and  installation' which 
are  pertinent  to  interpreting  the  performance  data  are  presented. 

In  Section  7.4.2,  the  results  from  these  radiometers  are  presented, 
and,  in  Section  7.4.3,  the  interpretation  of  these  results  is  pre¬ 
sented.  '  -! 
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7.4.1  Design  and  Installation 


A  narrow  view  angle  radiometer  as  defined  by  Figure  4-35  was 
installed  at  the  throat  plane  for  each  of  the  char  motor  nozzles. 

For  the  second  char  motor  firing,  the  radiometer  as  defined  by  Fig¬ 
ure  4-36  was  installed  at  an  A/A*  =  3.0.  Each  of  these  radiometers 
was  connected  to  a  nitorgea  gas  purge  flow  system.  In  addition,  the 
thermocouple  measuring  the  copper  sink  temperature  and  the  leads 
measuring  the  emf  across  the  constantan  sensor  were  connected  to  the 
recording  system  described  in  Section  4. 


The  gas  purge  plumbing  system  was  designed  for  a  nitrogen  gas 
flow"  rate  of  approximately  0.01  lbs/sec  for  all  the  radiometers. 

The  basic  consideration  in  selecting  this  flow  rate  was  that  it  pro¬ 
vides  ..  n  acceptable  cavity  temperature  (540°F)  as  determined  in  the 
Phase  II  APG  tests.  Since  no  overheating  of  the  radiometer  occurred 
in  the  APG  tests,  this  selection  of  purge  gas  flow  rate  for  the  char 
motor  firings  implies  that  no  overheating  would  occur  in  the  char 
motor  firings.  Before  eacn  test,  a  flow  check  was  performed  on  the 
nitrogen  purge  gas  system  to  assure  that  this  system  was  operational. 

7.4.2  Test  Results 

The  measured  results  from  the  narrow  view  angle  radiometer 
located  at  the  nozzle  throat  plane  of  the  first  char  motor  firing 
are  presented  in  Figure  7-23.  These  results  include  the  radiometer 
sink  temperature,  the  radiometer  millivolt  output,  and  the  gas  purge 
supply  pressure.  A  decrease  in  the  required  nitrogen  purge  gas  flow 
apparently  occurred  at  3  seconds  into  the  firing  followed  by  a  com¬ 
plete  loss  of  this  flow  at  15  seconds.  The  failure  at  15  seconds 
was  due  to  the  thermal  decomposition  of  the  supply  hose  caused  by 
the  high  radiation  flux  from  the  exhaust  plume.-  The  explanation  for 
the  event  at  3  seconds  is  not  so  clear-cut,  however.  Several  explan¬ 
ations  are  possible  as  outlined  below: 

•  The  unexpected  high  chamber ’pressure  resulted  in  a 
decrease  in  the  design  purge  gas  flow  and  therefore 
a  thermal  failure  of  the  radiometer  view  restrictor 
and  sapphire  window  (Figure  3-4).  The  nitrogen  purge 
gas  flow  is  controlled  by  a  sonic  orifice  so  that  the 
variations  in  local  static  pressure  (or  chamber  pressure) 
do  not  influence  the  flow  rate.  This  requires  a  supply 
pressure  slightly  more  than  twice  that  of  the  maximum 
local  static  pressure  (to  also  account  for  the  pressure 
drop  downstream,  of  the  sonic  orifice)  .  The  maximum 
local  static  pressure  was  620  psia  which  therefore 
indicates  a  nitrogen  purge  gas  supply  pressure  slightly 


FTGURE  7-23  RADIOMETER  OUTPUT  AND  ASSOCIATED  INFORMATION 
FOR  THE  FIRST-’  CHAR  MOTOR  FIRING  " 


greater  than  1240  psia  is  required.  The  actual  supply 
pressure  was  1250  psia,  indicating  the  potential  for  a 
decrease  in  the  purge  gas  flow. 

•  The  pyrolysis  products  from  the  purge  gas  hose  which 
failed  15  seconds  into  the  firing  plugged  up  the  sonic 
orifice  early  in  the  test  firing  and  thus  caused  radiom¬ 
eter  failure.  This  hose  probably  started  pyrolyzing 
very  early  in  the  firing.  These  pyrolysis  products 
could  have  condensed  on  the  orifice,  plugging  it  up 
partially  or  completely.  Actually,  the  orifice  was 
almost  completely  plugged  after  the  firing  but  this 
could  also  have  been  due  to  the  reverse  flow  after 

hose  failure. 

•  Dirt  in  the  purge  gas  supply  line  plugged  up  the  sonic 
orifice  early  in  the  firing  and  thus  caused  radiometer 
failure-  The  purge  gas  supply  lines  could  have  been 
dirty  and  some  of  this  dirt  could  have  been  carried 

to  the  orifice  and  plugged  it. 

i 

I 

It  is  not  possible  to  determine  which  of  the  above  explanations  is 
proper.  Based  on  these  possibilities,  however,  the  following  guide¬ 
lines  were  established  for  the  second  char  motor  firing: 

i 

•  The  purge  gas  supply  pressure  should  preferably  be 
twice  the  anticipated  maximum  chamber  pressure. 

•  All  lines  in  the  vicinity  of  the  char  motor  should 
be  stainless  steel. 

•  A  filter  should  be  installed  as  close  to  the  sonic 
orifice  as  practical. 

The  total  hemispherical  radiation  heat  flux  reduced  from  the  results 
shown  in  Figure  7-23  are  shown  in  Figure  7-24.  One  notes  that  the 
magnitude  of  the  measured  radiative  flux  prior  to  the  loss  of  the 
nitrogen  purge  was  on  the  order  expected  and  that  the  transient 
response  of  the  radiometer  was  similar  to  the  chamber  pressure  trans¬ 
ient  which  was  also  expected. 


The  above  mentioned  recommendations  were  incorporated  in  the  gas 
purge  system  for  the  second  char  motor  firing,  and  the  measured 
results  for  the  two  radiometers  are  shown  in  Figures  7-25  through 
7-27.  The  hemispherical  radiation  heat  flux  obtained  from  the  data 
reduction  procedure  is  presented  in  Figure  7-28.  [The  only  suspicious 
c-ata  from  these  results  are  the  measured  emf  outpiit  from  the  radi¬ 
ometer  located  at  the  nozzle  throat  plane.  One  notes  from  Figure  7-27 
that  the  initial  emf  output  from  both  radiometers  have  the  same  trend 
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as  the  chamber  pressure.  However,  at  approximately  eight  seconds 
through  the  firing,  the  output  from  the  throat  radiometer  decreased 
but  then  started  to  increase  again  at  approximately  twelve  seconds 
into  the  firing.  The  output  took  another  dip  at  approximately  42 
seconds  but  achieved  its  maximum  value  at  motor  burnout.  Both  the 
output  from  the  exit  cone  radiometer  and  sink  temperature  and  purge 
pressure  results  did  not  exhibit  these  anomalies.  Before  discussing 
possible  reasons  for  the  variations  shown  in  the  output  from  the 
throat  radiometer,  thte  following  observations  should  be  noted: 

•  The  output  from  both  radiometers  delayed  in  a  similar 
manner  at  motor  burnout  and  lagged  the  combustion  chamber 
pressure.  This  lag  was  expected  as  the  temperature  of 
the  nozzle  surface  which  is  viewed  by  the  radiometer 
decreases  more  slowly  than  the  chamber  pressure. 

•  Post-test  observation  of  both  radiometers  showed  that 
no  effects  of  excessive  heating  were  visible  and  that 
the  sapphire  lens  was  clear. 

■  £ 

•  Post-test  observation  of  the  exit  cone  cavity  showed 
that  it  did  not  contain  any  debris.  However,  the  side- 
walls  of  the  cavity  at  the  nozzle  throat  had  some  slag 
which  extended  to  the  radiometer,  but,  at  motor  burnout, 
this  slag  was  not  restricting  the  aperture  or  the  view 
of  the  radiometer . 

Based  on  the  above  observations,  one  possible  explanation  for  the 
variations  shown  by  the  output  of  the  nozzle  throat  -radiometer  is 
that  during  the  firing  either  ablation  products  or  alumina  particles 
were  falling  into  the  cavity  and  adhering  to  the  relatively  cool 
cavity  walls.  However,  the  purge  flow  kept  this  debris  from  con- 
tamininating  the  radiometer.  If  the  debris  continued  to  build  up 
on  the  cavity  walls,  it  would  restrict  the  view  of  the  radiometer 
and  cause  the  dips  shown.  If  this  debris  broke  off  and  was  blown 
out  by  the  purge  gas,  then  output  from  the  radiometer  would  increase. 
The  post-test  observations  substantiate  the  validity  of  the  maximum 
value  measured  at  motor  burnout.  Some  interpretations  concerning 
this  data  point  plus  the  output  from  the  exit  cone  radiometer  are 
presented  in  the  next  section. 


7.4.3  Interpretation  of  the  Radiation  Results 


The  validity  of  converting  the  measured  radiative  heat  flux 
from  a  narrow  view  angle  radiometer  to  a  total  hemispherical  radia¬ 
tive  heat  flux  was  discussed  in  Section  5.  In  that  section,  it 
was  shown  that  the  total  hemispherical  radiative  flux  is  a  true 
measurement  of  the  flux  being  emitted  by  the  radiative  source  for 
the  following  situations: 
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•  The  particle  laden  combustion  products  are  a  non-absorbing 
and  non-scattering  media,  and,  thus,  the  radiometer  views 
only  the  opposite  wall. 

•  The  particle  laden  combustion  products  form  a  constant  tem¬ 
perature  optically  thick  media  which  has  a  negligible  scat¬ 
tering  coefficient  but  has  an  absorption  coefficient  which 
is  independent  of  wavelength.  In  this  situation,  the  com¬ 
bustion  products  emit  radiation  as  a  gray  body 


At  the  nozzle  exit  cone  station  corresponding  to  an  A/A*  =  3.0,  the 
emissive  power  of  the  opposite  wall  being  viewed  by  the  radiometer 
corresponds  to  a  hemispherical  radiative  flux  of  approximately  380 
Btu/ft^sec  if  the  surface  emissivity  is  unity.  If  the  combustion 
products  emit  as  a  black  body  at  the  local  static  temperature,  the 
average  hemispherical  radiative  flux  through  the  firing  would  be 
540  Btu/ft^sec.  The  similar  values  at  the  nozzle  throat  are  420 
Btu/ft^sec  from  the  opposite  wall  and  about  900  Btu/ft2sec  from 
black  body  combustion  products.  The  maximum  measured  values  at  the 
exit  cone  station  and  throat  were  460  and  540  Btu/ft2sec,  respec¬ 
tively,  which  fall  between  the  extremes  identified  above.  These 
comparisons  show  that  the  particle  laden  combustion  products  are 
emitting  a  significant  amount  of  energy  in  the  form  of  radiation. 
However,  these  comparisons  do  not  reveal  whether  the  radiative  flux 
absorbed  by  the  sensor  is  primarily  from  the  opposite  wall  with  a 
small  contribution  from  the  particle  laden  combustion  products  (low 
effective  emissivity)  or  primarily  from  the  particle  laden  combus¬ 
tion  products  alone  (high  effective  emissivity) .  An  attempt  to  de¬ 
fine  which  is  the  valid  possibility  is  made  in  the  following  para¬ 
graphs  . 


To  determine  the  relative  significance  of  the  radiation  from 
the  particle  laden  combustion  products,  one  must  account  for  the 
contributions  from  both  the  gases  and  the  alumina  particles.  Pre¬ 
liminary  calculations  showed  the  radiation  from  the  gas  was  negli¬ 
gible.  However,  these  calculations  showed  that  the  radiation  from 
the  particles  was  significant.  These  calculations  are  summarized 
be  low . 


The  first  step  in  predicting  the 
was  to  define  the  optical  depth  which 


radiation  from  the  particles 
is  given  by 


(7-4) 
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Using  the  procedure  outlined  in  Reference  7-2,  the  extinction  coef¬ 
ficient  can  be  defined  by 


6 


\ 


N  o. 
P  t 


(7-5) 


where  -  particle  density 

o  -  extinction  cross-section 

Substituting  Equation  (7-5)  into  Equation  (7-4)  gives 


\ 


(7-6) 


The  particle  density  can  be  defined  from  the  aluminum  loading  and  by 
assuming  the  mean  particle  radius  is  given  bv  the  relationship 
(Reference  7-2) 


r  =  0.862  (In  P  -  4.17)  (7-7) 

pm  c 

where  r  -  mean  particle  radius  in  microns 
pm 

P  -  combustion  chamber  pressure 
in  psia 


The  extinction  cross-section,  ot,  can  be  calculated  using  Mie  theory 
(Ref.  7-2)  provided  the  real  and  imaginary  parts  of  the  alumina  in¬ 
dex  of  refraction  are  known.  The  calculation  of  these  cross-sections 
for  the  particular  particle  concentrations  in  the  char  motor  firing 
was  beyond  the  scope  of  this  program.  However,  for  the  preliminary 
calculations  presented  here,  these  cross-sections  can  be  estimated 
from  the  information  presented  in  Reference  7-2.  In  this  reference, 
it  was  shown  that  the  extinction  efficiency  ractor  which  is  defined 
by 


(7-8) 
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is  approximately  equal  to  two  if  the  ratio  2irr  /A  is  also  greater 
than  two.  Since  this  condition  is  satisfied,  the  extinction  cross- 
section  was  approximated  from  Equation  (7-8)  by  using  Equation  (7-7) 
to  calculate  the  mean  particle  density  and  by  assuming  an  extinction 
efficiency  factor  of  two.  The  results  of  this  calculation  including 
the  effects  of  nozzle  surface  recession  and  chamber  pressure  decay 
are  tabulated  below  in  Table  7-9.  Also  shown  in  the  table  is  the 
blackbody  emissive  power  based  on  the  local  static  temperature  of 
the  combustion  products. 

TABLE  7-9 

OPTICAL  DEPTH  AND'  BLACKBODY  EMISSIVE  POWER 


Chamber 

Pressure 


(psia) 


Area  Ratio  of 
Radiometer  a 


Exit 

Cone 


Optical 

Depth 

Blackbody  Emissive 
Power 

(Btu/ft2sec) 

No2zle 

Throat 

Exit 

Cone 

Nozzle 

Throat 

Exit 

Cone 

21.4 

6.92 

951 

495 

20.0 

6.86 

924 

510 

16.5 

8.85 

805 

585 

aChanging  area  ratio  for  the  exit  cone  location  due  to  surface 
recession  that  occurred  at  the  throat. 

Referring  to  Equation  (5-8)  ,  the  radiation  from  the  opposite  wall 
is  attenuated  by  the  factor  e-T/^.  If  p  =  1  (a  good  approximation 
for  a  narrow  view  angle  radiometer) ,  the  values  for  the  optical  - 
depth  in  Table  7-9  indicate  that  the  radiometer  does  not  see  the 
opposite  wall  and  that  the  measured  radiation  is  due  to  the  alumina 
particles.  This  conclusion  was  further  substantiated  by  solving  for 
the  emissivity  of  a  one-dimensional  slab  of  scattering  and  emitting 
particles.  This  emissivity  is  given  by  (Reference  7-2) 


(i)=  2M  £ e_K ' T (k 1  +  m) 


m)-2K^J 


(7-9) 


where 


D  =  (k*  +  m)2  eK'T  -  ( K' 
K'=  |  M ^2  (1-M)  B  +  m]  j*5 


v,-\  2  -K  '  t 
M;  e 


a  +  a 
a  s 


•  5>- 


/ 


o&,  a  -  absorption  and  scattering 
^  cross-sections 

B  -  fraction  of  radiation 
scattered  backward 


This  equation  was  solved  for  various  values  of  B  and  M  which  would 
bound  the  conditions  of  the  test  firing,  and  the  results  are  shown 
in  Figure  7-29.  For  the  previously  mentioned  values  of  the  optical 
depth,  the  emissivity  is  approximately  constant  for  a  fixed  value 
of  B  and  M.  The  similarity  in  the  exit  cone  of  the  variations  in 
the  measured  hemispherical  radiative  flux  and  in  the  predicted  black- 
body  emissive  power  substantiated  the  conclusion, th£t  the  emissivity 
of  the  combustion  products  is  approximately  constant  for  this  varia¬ 
tion  in  the  optical  depth.  In  addition,  the  range  of  emissivity  in 
Figure  7-29  is  consistent  with  the  measured  results.  To  be  more 
definitive  requires  an  extensive  analysis  to  evaluate  B  and  H  and, 
as  previously  mentioned,  this  is  beyond  the  scope  of  the  present 
contract.  ...  t 


In  summary,  the  radiometer  at  the  A/A*  =  3.0  location  apparently 
provided  a  valid  indication  of  the  tota^  radiative  heat  flux  through¬ 
out  the  test  firing.  This  radiation  was  from  the  particle  stream 
with  an  effective  emissivity  near  unity.  The  radiometer ,  at  the 
throat  location  apparently  indicated  a  significantly  lower  than 
actual  radiative  flux  during  all  but  the  early  and  later  parts  of 
the  test  firing.  Apparently  an  unexplained  partial  restriction  of  . 
the  radiometer  view  occurred. 
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SECTION  8 

CONCLUSIONS  AND  RECOMMENDATIONS 


Total  heat  flux  and  radiative  heat  flux  measurement  systems 
for  use  in  ablative  solid  rocket  nozzles  were  developed  and  demon¬ 
strated.  The  total  flux  system  is  a  thermocouple  instrumented  plug 
of  the  same  material  and  layup  angle  as  the  material  in  which  it  is 
installed.  The  radiative  heat  flux  system  is  a  narrow  view  angle, 
Gardon-type  radiometer.  The  final  evaluation  and  demonstration 
tests  were  two  nozzle  firings  performed  on  the  AFRPL  40-inch  char 
motor.  The  conclusions  and  recommendations  related  to  the  develop¬ 
ment  and  demonstration  of  these  systems  are  as  follows: 

•  Total  heat  flux  measurement  system  provides  a  repeatable 
and  accurate  measurement  of  total  heat  flux 

•  Radiative  heat  flux  measurement  system  provides  an  accu¬ 
rate  measurement  of  total  radiative  heat  flux  provided 
the  assumption  of  a  uniform,  grey-body  source  is  satisfied 

•  These  systems  have  been  designed  for  use  in  a  wide  range 
of  materials  and  firing  conditions,  including  high  chamber 
pressure 

•  The  only  limitation  to  the  total  heat  flux  system  is  that 
the  material  be  well  characterized  in  terms  of  thermal  and 
ablative  properties 

•  The  only  limitation  for  the  radiative  heat  flux  system  is 
that  the  installation  design  does  not  result  in  a  large 
thermal  load  to  the  radiometer  body 

•  The  reason  for  an  apparent  partial  obstruction  of  the 
radiometer  cavity  at  the  throat  of  the  second  nozzle  should 
be  investigated  and  any  necessary  corrective  action  defined 

•  A  radiative  heat  flux  measurement  system  which  measures 
total  (hemispherical)  radiation  heat  flux  directly  should 
be  developed. 
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